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Introduction 


The  contract  period  for  the  University  Research  Initiative  was 
from  September  15,  1986  to  September  30,  1989.  The  Pennsylvania 

State  University  served  as  the  prime  contractor  with  A.  J.  Ferraro 
as  Principal  Investigator;  subcontracts  were  awarded  to  the 
University  of  California  at  Los  Angeles  (A.  Y.  Wong)  and  S.A.I.C. 
(D.  Papadopolous) .  Reports  of  these  subcontractors  are  attached  to 
this  final  report.  Support  from  the  Geophysical  Institute  (J.  V. 
Olson)  was  provided  by  a  separate  ONR  contract  and  a  final  report 
was  submitted  earlier.  A  broad  definition  of  the  tasks  of  the 
above  institutions  are: 

Penn  State  -  Develop  the  science  utilizing  the  HIPAS  facility. 

U.C.L.A.  -  Operation  of  _he  HIPAS  facility  for  cooperative 

experiments. 

S.A.I.C.  -  develop  theoretical  concepts  related  to  ionospheric 
modification. 

Geophysical  Institute  -  Provide  information  on  the  development 

of  electro jet  current  events. 

There  were  four  experimental  campaigns  established  during  this 
contract  covering  the  period 

June  15-25,  1987 

October  13-23,  1987 

July  1 1 2 2  ,  1988 

July  12-  Aug.  8,  1989 

The  performance  summary  of  the  HIPAS  facility  is  described  in 
detail  in  the  U.C.L.A.  attachment  to  this  final  report. 


Overview 


The  broad  research  objective  of  this  program  has  been  to 
investigate  unique  concepts  for  generating  ULF,  ELF  and  VLF  using 
the  ionosphere  as  the  '^wireless'*'  antenna.  The  "Wireless  Antenna" 
or  HF  demodulation  experiment  was  demonstrated  by  A.  J.  Ferraro^in 
1980  using  the  Arecibo  Observatory's  heating  facility.  Modulation 
of  the  dynamo  current  system  at  that  latitude  was  the  major  source 


for  the  "wireless"  antenna. 

Some  of  the  major  experimental  accomplishments  are 

-  Excellent  quality  of  generated  ELF  and  VLF  signals  . 

-  Determination  of  D-region  electron  density  profiles  at  high 
latitudes  for  large  sunspot  numbers. 

-  slow  beam  steering  to  examine  the  structure  of  the 
electro jet  current  system  mapped  down  into  the  D-region 

-  determination  of  electro jet  activity  from  ELF/VLF 

polarization  characteristics 

-  created  multiple  VLF  ionospheric  sources  (two-element  VLF 
arrays) 

-  Fast  beam  steering  -  the  painting  technique;  a  potential 
method  of  increasing  ELF  power. 

-  Short  range  field  measurements  in  Alaska  with  the  Penn  State 
Mobile  Van. 

-  Demonstrated  the  stability  of  the  current  system  to  support 
bi-phase  PSK  and  Quad-phase  PSK  to  increase  data  rates. 


Based  upon  the  July  1988  campaign  it  was  estimated  that  the  ELF 
power  produced  by  HIPAS  for  one  spot  exceeded  one  milliwatt  80%  of 


the  time.  This  figure  was  arrived  at  by  taking  2000  Hz  data  and 
scaling  to  100  Hz. 


Geomagnetic  and  Solar  Conditions  during  the  HIPAS  Campaigns 

Attached  are  summarized  data  for  the  solar  flux  and  earth's 
magnetic  field.  One  can  refer  to  these  to  estimate  the  ionospheric 
conditions  during  the  experiment  times  that  some  campaigns  were 
conducted  during  the  increasing  active  phase  of  the  solar  cycle. 

AGARD  Symposium 

Attached  is  a  preprint  of  a  paper  presented  at  the  Spring 
1990,  Electromagnetic  Wave  Propagation  Panel  Symposium  held  in 
Bergen,  Norway. 

List  of  Titles  and  Abstracts  of  Key  Scientific  Achievements 


October  1987 
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Ottawa  10.7  cm  solar  flux  for  October  1987.  The  measurements 
are  adjusted  to  1  AU,  measured  at  17:00  UTC  daily  and 

expressed  in  units  of  10‘22  watts/m  ‘/hertz.  Data  taken  from  the 
National  Geophysical  Data  Center. 


July  1988 
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Ottawa  10.7  cm  solar  flux  for  July  1988.  The  measurements  are 
adjusted  to  1  AU.  measured  at  17:00  UTC  daily  and  expressed 

in  units  of  10‘2  watts/m7hertz.  Data  taken  from  the  National 
Geophysical  Data  Center. 
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Planetary  K  indices  for  October  1987.  Each  shaded  area  shows 
3-hour  range.  Data  taken  fmra  the  National  Geophysical  Data 
Center. 
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Planetary  K  indices  for  juiv  1988.  Each  shaded  area  shows 
3-hour  range.  Data  taken  from  the  National  Geophysical  Data 
Center. 


An  Experimental  Method  to  Determine  the 
Phase  Height  of  an  Artificially  Created  Ionospheric 

ELF/VLF  Source 

This  report  describes  the  design  and  operation  of  a  system 
which  can  determine  the  phase  height  of  extremely  long  wavelength 
electromagnetic  emissions  from  natural  ionospheric  currents  which 
have  been  modulated  by  radio  wave  heating.  Previously, 
experimental  campaigns  at  Arecibo  have  utilized  collocated  heater 
transmitters  and  ionospheric  signal  receivers  which  greatly 
simplifies  the  task  of  determining  a  phase  relationship  between  the 
heater  modulation  wave  and  the  received  signal  wave.  However, 
future  experiments  envision  extensive  portable  receiver  operation 
at  sites  quite  remote  and  distant  from  the  heater  transmitter. 

A  low  cost  portable  system  is  developed.  The  system  consists 
of  two  stages  for  operation  at  the  HIPAS  heater  site  and  the 
receiver  site.  Prior  to  an  experimental  session  the  two  systems 
are  compared  to  determine  an  initial  phase  relationship.  One 
system  then  provides  the  modulation  to  the  heater  wave  while  the 
other  provides  a  reference  signal  at  the  receiver  site.  A 
programmable  frequency  synthesizer  is  also  incorporated  into  the 
system  which  will  coexist  with  the  existing  equipment  controllers. 

A  function  relating  signal  phase  shift  to  source  phase  height 
is  developed  for  the  general  case  of  an  oblique  incidence  heater 
beam.  This  function  is  reduced  for  the  considered  case  of  a  normal 
incidence  heater  beam.  Selected  source  phase  height  results  from 
an  experimental  campaign  at  HIPAS  are  presented  to  illustrate  the 
systems  performance  in  actual  use  (figures  attached) . 


[m  @  24:40  A.S.T. 


An  Investigation  of  the  Obliquely  Downgoing 
ELF/VLF  Generated  Waves  During  Ionospheric 
Modification  Heating  Experiments  at  High  Latitudes 

This  document  has  reported  on  a  series  of  ionospheric 
modification  heating  experiments  conducted  at  the  High  Power 
Auroral  Simulation  (HIPAS)  facility  near  Fairbanks,  Alaska  over  a 
period  of  2  years  from  June  1987  to  August  1989.  Particular 
emphasis  is  devoted  to  the  study  of  the  polarization  of  the 
electric  field  of  the  obliquely  downgoing  ELF/VLF  electromagnetic 
waves  generated  by  irradiation  of  the  heated  volume  from  a  high 
power  modulated  HF  wave. 

A  single  slab  model  of  the  ionosphere  has  been  chosen  to 
simplify  the  solution  of  the  Booker  quartic  which  is  used  to 
compute  the  wave  path  and  direction.  The  ELF/VLF  source  is 
simulated  by  two  horizontal  current  sources;  one  in  the  east-west 
direction,  and  one  in  the  north-south  direction  due  to  modulation 
of  the  electrojet  which  flows  at  an  altitude  of  100  km  to  110  km. 
The  waves  resulting  from  these  two  orthogonal  sources  are 
superimposed  on  the  ground  to  give  the  tilt  angle  and  ellipticity 
of  the  incoming  waves  which  are  the  parameters  being  measured. 

The  polarization  ellipse  gives  a  very  good  indication  of  the 
presence  of  an  electro jet  event.  The  polar imeter  is  sensitive  even 
with  weak  ELF/VLF  returns.  During  all  the  experiments,  it  is  able 
to  provide  instantaneous  visual  information  on  the  ground  and  can 
provide  guidance  as  to  the  scientific  measurement  schedule  during 
campaigns.  When  pulsations  were  first  encountered,  it  showed  a 
very  distinct  signature.  This  is  a  big  advantage  over  the 


magnetogram  records  which  are  not  immediately  available  during  the 
course  of  the  experiments. 

It  has  been  shown  that  the  polarimeter  is  a  viable  part  of  the 
ELF/VLF  scientific  program  in  furnishing  diagnostics.  It  can  also 
provide  a  rough  estimate  of  electron  density  and  collision 
frequency.  For  certain  electron  densities  and  collision 
frequencies,  the  tilt  angle  of  the  polarization  ellipse  is  shown  to 
track  the  angular  displacement  of  the  horizontal  source  current. 
This  allows  the  angular  direction  of  the  source  current  to  be 
determined. 

The  attached  figure  shows  the  tilt  angle  of  the  polarization 
ellipse  versus  source  angular  displacement.  For  a  reasonable  range 
of  electron  collision  frequencies  the  expected  observed  tilt  does 
track  the  source  angular  displacement.  For  those  cases,  the 
observed  ground  level  polarization  is  a  reasonable  indicator  of 
current  direction  flow  in  the  ionosphere. 


Source  angular  displacement  (degree) 


An  Investigation  of  the  Polar  Electro jet 
Current  System  Using  Radio  Wave  Heating  of  the 

Lower  Ionosphere 

The  High  Power  Auroral  Stimulation  (HIPAS)  heating  facility 
has  been  used  to  modulate  D-region  ionospheric  currents  at  high 
latitudes,  producing  extremely  low  frequency  (ELF)  radio  wave 
emissions.  The  behavior  of  these  ionospheric  currents  can  be 
deduced  from  a  comprehensive  study  of  the  ELF  signals  received  at 
a  local  field  site.  This  document  examines  the  relationship 
between  the  ELF  magnetic  field  strength  measured  on  the  ground  and 
the  intensity  of  an  overhead  electro jet  current.  The  mapping  of 
the  polar  electro jet  current  from  the  E  region  down  through  the  D 
region,  where  it  can  then  be  modulated  by  the  heater  beam,  is 
investigated.  A  finite  difference  solution  to  the  electrojet 
mapping  problem  is  presented  in  which  arbitrary  conductivity 
profiles  can  be  specified.  Results  have  been  obtained  using  a 
simple  Cowling  model  of  the  electrojet.  These  results  indicate 
that  for  an  electro  jet  flowing  at  an  altitude  of  100  km  with  a 
scale  size  in  excess  of  100  km,  the  ^mapping  of  the  horizontal 
current  density  can  be  completely  characterized  in  terms  of  the 
Pedersen  and  Hall  conductivities.  This  indicates  that  the  mapping 
becomes  independent  of  scale  sizes  which  exceed  100  km.  The 
downward  mapping  of  the  electro jet  current  and  associated  electric 
field  in  the  presence  of  ionospheric  conductivity  irregularities 
has  been  studied.  It  was  demonstrated  that  D  region  conductivity 
irregularities  do  not  significantly  change  the  downward  mapping  of 
ionospheric  electric  fields.  However,  these  conductivity 
irregularities  were  found  to  have  a  major  influence  upon  the 


downward  mapping  of  ionospheric  currents.  A  promising  new 
diagnostic  technique,  for  studying  ionospheric  D  region  currents, 
has  been  implemented  using  the  HIPAS  facility.  This  technique 
involves  HF  beam  steering  for  localized  ELF  generation  in  the 
mapped  region  below  electrojets.  Beam  steering  has  been  used  to 
deduce  the  strength  and  current  distribution  of  the  polar 
electrojet,  and  for  charting  the  movements  of  overhead  currents. 
The  attached  figure  shows  ELF  relative  field  strength  for  a  conical 
scan.  The  dashed  curve  is  expected  ELF  return  from  a  uniform 
current  sheet  while  the  solid  curve  is  due  to  the  actual  overhead 
current  distribution.  The  difference  between  these  two  curves  is 
a  measure  of  the  horizontal  variation  in  current  density.  The  next 
figure  is  a  deduced  ionospheric  current  map  showing  current  density 
strengths  over  the  region  of  conical  scan. 


ical  scan  magnitude  data.  I  lie  data  Tile  was  CNS(,AN2I.29,  the  start  time 
04:52.00,  and  the  date  recorded  was  7/21/88 


Determination  of  D-region  Electron  Densities 
Using  the  Stepped  Frequency  Technique  with 
—  Ionospheric  Heating 

This  document  describes  a  new  technique  for  deriving  D-region 
electron  densities  from  the  Extremely  Low  Frequency  (ELF)  signals 
produced  by  heating  of  the  D-region  of  the  ionosphere.  A  catalog 
approach  was  employed  in  this  technique.  A  class  of  exponential 
density  models  were  used  to  create  a  catalog  of  transmission 
coefficients  -deduced  from  the  full  wave  theory,  and  stored  in  a 
file.  The  electron  density  profiles  were  derived  by  searching  the 
catalog  with  the  transmission  coefficient  extracted  from  the 
locally  received  ELF  signals  and  selecting  the  best  fit  model  to 
the  observation.  In  particular,  the  degree  of  electron  density 
ionization  of  the  D-region  at  high  latitudes  can  be  studied.  The 
heating  facility  near  Fairbanks,  Alaska,  known  as  the  HIPAS 
facility  was  used  for  this  work.  In  this  document  the  theory  of 
electron  density  synthesis  is  described,  results  of  actual  ELF  data 
presented,  and  preliminary  electron  density  profiles  determined. 
A  new  future  ionospheric  measurement  technique  FDTD  is  discussed. 

Attached  are  figures  giving  examples  of  stepped  frequency  data 
and  derived  D-region  electron  densities. 
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Example  of  stepped  data.  A  sample  of  the  magnitude  of  received  high  band 
stepped  frequency  data  taken  at  MIPAS  on  10/14/87  at  17:08  to  17:12  A.S.T.. 
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Electron  density  profiles.  1.  Mitro-Rowe  density  profile;  2.  density  profile  for 
17:28:00  of  10/14/87  AKT;  3.  density  profile  for  early  morning  of  10/16/87  AKT; 
4.  density  model  used  by  Barr. 


Development  of  a  Diagnostic  System  for 
Ionospheric  Modification  Studies  at  High  Latitudes 

This  report  describes  the  design  and  operation  of  a  diagnostic 
system  for  the  study  of  extremely  long  wavelength  electromagnetic 
emissions  from  natural  ionospheric  currents  which  have  been 
modulated  by  radio  wave  heating.  The  ionosphere  at  high  polar 
latitudes  is  complex  in  structure  and  dynamic  in  behavior, 
particularly  during  periods  of  auroral  activity  and  geomagnetic 
disturbance.  It  also  contains  strong  current  systems  which  hold 
great  potential,  when  properly  modulated,  as  radiators  of  ULF,  ELF, 
and  VLF  signals.  Two  experimental  campaigns  have  been  conducted  to 
investigate  the  generation  and  reception  of  signals  produced 
through  modulated  ionospheric  heating  from  the  High  Power  Auroral 
Stimulation  transmitters  near  Fairbanks,  Alaska. 

The  diagnostic  system  includes  a  dual-channel  receiver  capable 
of  operation  at  frequencies  below  8  kHz,  and  companion 
microcomputers  for  real-time  data  acquisition  and  processing. 
Software  has  been  developed  to  control  data  collection,  analysis, 
and  storage  in  a  "user-friendly"  fashion?  several  configurations 
and  operational  modes  provide  maximum  flexibility  under  a  variety 
of  experimental  conditions.  Specific  provisions  are  made  for 
coherent  demodulation  and  high-speed  digital  sampling  of  incoming 
signals. 

Additional  computer  control  the  transmitter  modulating 
frequency  and  receiver  tuning.  These  permit  synchronized, 
operation  of  the  transmitter  and  receiver  sites  according  to  a 
preselected  schedule,  reducing  the  need  for  manual  intervention  and 


permitting  rapid  reselection  of  operating  frequency.  Selected 
results  from  the  experimental  campaigns  are  presented  to  illustrate 
the  system's  capabilities  and  performance  in  actual  use. 

A  more  detailed  discussion  is  allocated  to  this  report  since 
the  receiver  plays  an  important  role  in  developing  the  science. 

Loop  Antennas 

Loop  antennas  are  simple  to  construct  and  provide  stable 
operation  in  -  all  climates.  With  the  addition  of  capacitive 
loading,  the  resonant  frequency  of  the  antenna  system  may  be 
lowered  to  obtain  voltage  amplification  of  the  output  signal  and 
improved  selectivity.  A  center- tap  allows  connection  to  balanced 
transmission  line,  and  the  azimuthal  pattern  variation  lends  itself 
well  to  measurement  of  wave  direction  and  polarity.  Three  types  of 
loop  antennas  are  available  for  use  with  the  reconf  igurable 
receiving  system. 

The  first  set  of  loops  was  used  in  previous  experimental 
campaigns,  and  each  consists  of  2  00  turns  of  number  14  wire 
(center- tapped)  upon  a  1  meter  diameter  wooden  form.  The  coils 
have  a  measured  DC  resistance  of  approximately  7  ohms  and  an 
inductance  of  .075  henry  with  a  self-resonant  frequency  near  8000 
Hz.  For  calibration,  a  single-turn  loop  of  0.5  m  diameter  is 
mounted  coaxially  at  a  distance  of  0.5  m.  The  drive  signaj.  is 
applied  through  a  50  ohm  series  resistor,  and  the  current  in  the 
calibration  coil  is  obtained  by  measurement  of  the  voltage  drop 
across  the  resistor.  Mathematical  formula  provide  a  relationship 
for  the  mutual  inductance  between  two  coaxial  loops  used  for 
calibration. 


The  considerable  size  and  weight  of  these  units  necessitated 
an  attempt  to  simplify  the  logistics  and  expenses  involved  in 
transportation  of  the  receiver  system.  A  second  set  of  loop 
antennas  was  constructed  with  800  turns,  of  number  28  wire  upon  0.5 
m  diameter  wooden  forms.  The  product  of  the  number  of  turns  and 
loop  cross-sectional  area  was  not  altered,  in  order  to  achieve 
—  electrical  operation  comparable  with  the  original  units.  The 
redesigned  antenna  could  be  easily  carried  by  one  person, 
eliminating  the  need  for  special  packaging  and  handling. 

Following  construction  of  wooden  coil  forms  and  scramble-wound 
center-tapped  loops,  the  DC  resistance  was  measured  as  98  ohms  and 
the  inductance  as  0.3  henry.  In  order  to  simplify  the  calibration 
process,  a  single-turn  calibration  loop  was  wound  directly  on  top 
of  the  receiving  loop  itself;  instead  of  the  complicated  mutual- 
inductance  expression  referred  to  above,  the  new  design  could  now 
be  treated  as  an  ideal  transformer.  Frequency  response 
measurements  revealed  that  these  loops  resonated  at  approximately 
1200  Hz  due  to  the  increased  inductance  and  self-capacitance  of  the 
smaller  design.  As  expected,  the  increased  resistance  of  a  greater 
number  of  turns  would  with  smaller  wire  also  reduced  the  Q  factor. 

For  ULF  reception  below  10  Hz,  several  loop  antennas  of  20,000 
turns  upon  a  mu-metal  core  were  obtained  from  geophysical 
researchers  at  the  University  of  Texas.  At  these  very  low 
frequencies,  calibration  is  performed  in  an  incident  magnetic  field 
produced  by  the  rotation  of  a  motor-driven  precision  bar  magnet. 
Although  a  very  long  and  slender  core  having  a  large  ratio  of 
length  to  diameter  will  maximize  the  permeability,  several  factors 


render  this  approach  impractical.  These  include  the  mechanical 
integrity  of  the  resulting  structure,  the  effect  of  elastic  strains 
upon  the  core's  magnetic  properties,  and  possible  desentization  by 
the  earth's  geomagnetic  field.  The  latter  is  due  to  magnetic 
saturation  of  the  core  reducing  the  effective  permeability, 
yielding  an  antenna  sensitivity  which  varies  greatly  with 
geographic  orientation. 

In  order  to  minimize  these  effects,  the  length  to  diameter 
ratio  is  chosen  to  be  much  smaller  than  the  theoretical  optimum, 
forcing  the  effective  core  permeability  to  be  dependent  solely  upon 
this  factor.  The  ULF  loop  antennas  satisfy  the  criteria  for  a 
"geometry  limited"  design  which  performs  in  a  stable  and 
predictable  manner. 

Antenna  Tuning  and  Transmission  Line  Interconnection 

The  loop  antenna  may  be  treated  as  a  parallel  LC  tank  circuit, 
and  its  resonant  frequency  lower  by  an  increase  in  capacitance 
through  the  external  addition  of  parallel  capacitors.  The  mu-metal 
core  ULF  antenna  is  operated  with  a  single  capacitor  across  its 
terminals  at  all  times,  but  the  air-core  ELF/VLF  receiving  loops 
are  tuned  to  the  specific  resonant  frequency  of  interest. 

The  original  technique,  developed  for  earlier  experimental 
campaigns,  utilized  a  large  rotary  switch  at  the  antcr.na  cite  which 
sequentially  connected  preselected  sets  of  mylar  and  silver  mica 
fixed  capacitors.  Assembly  of  these  capacitor  combinations  was 
largely  trial-and-error ,  necessitating  measurement  of  frequency 
response  following  each  modification.  Operation  was  limited  to 


pretuned  channels,  and  manual  selection  of  the  proper  switch 


setting  was  necessary.  Since  the  tuning  assembly  was  remotely 
located  at  the  antenna  loop,  adjustment  required  the  dispatch  of  an 
operator. 

Although  originally  acceptable,  this  method  became 
increasingly  unwieldy  when  a  comprehensive  program  of  stepped 
frequency  measurements  was  begun.  The  ability  to  rapidly  retune 
the  antennas  was  required  in  order  to  permit  operation  on  multiple 
frequencies  within  a  short  span  of  time;  otherwise,  the  dynamic 
nature  of  the  ionosphere  would  preclude  valid  comparison  of  results 
obtained  at  different  frequencies.  The  rotary  switch  tuning  box 
was  replaced  by  a  relay  switching  unit  which  contains  a  binary 
sequence  of  sixteen  composite  capacitors  connected  in  parallel 
through  relay  contacts.  A  16-bit  binary  number  is  placed  upon  the 
relay  control  lines;  each  energized  relay  switches  in  a  capacitance 
twice  that  of  the  preceding  (less  significant)  bit  position.  The 
resultant  combination  can  synthesize  any  required  value,  limited  in 
range  by  a  maximum  of  all  capacitors  simultaneously  in-circuit  and 
limited  in  resolution  by  the  capacitance  of  the  initial  increment. 

The  total  capacitance  at  the  tuning  unit  terminals  is  equal  to 
the  product  of  the  binary  value  upon  the  control  inputs,  and  the 
capacitance  associated  with  the  least  significant  bit  position.  As 
built,  the  base  unit  is  approximately  20  pf,  and  the  maximum  value 
available  is  approximately  1.31  uf.  These  values  permit  tuning  of 
the  1  m  diameter  loops  over  an  approximate  range  of  800  Hz  to  8000 
Hz,  and  the  0.5m  diameter  loops  over  200  Hz  to  1200  Hz.  In  order 
to  ensure  a  monotonic  sequence  without  missing  values,  the 


capacitance  associated  with  each  bit  position  was  carefully  trimmed 
on  an  impedance  bridge  equipped  with  a  digital  readout. 


Equipment  Control  at  Transmitter  and 
Receiver  Sites 

The  primary  functions  of  the  main  data  processing  system  are 
real-time  data  acquisition,  analysis,  display,  and  storage;  in 
actual  operation,  a  minimum  amount  of  processor  time  remains  for 
the  performance  of  other  tasks.  The  video  display  and  keyboard  are 
dedicated  to  instantaneous  display  of  incoming  data  and  interactive 
operator  communications.  The  assignment  of  additional  functions 
would  then  prove  detrimental  to  the  overall  system  operation. 
Further,  32  bits  of  parallel  output  are  required  for  control  of 
each  synthesizer's  frequency,  and  16  bits  are  used  to  set  the 
antenna  loading  for  each  tuned  loop;  as  equipped,  the  8088-based 
system  lacked  this  extensive  parallel  input-output  capacity. 

For  these  reasons,  two  Apple  II  computers  used  in  earlier 
experimental  campaigns  were  restored  to  service  as  equipment 
controllers  for  the  transmitter  and  receiver  sites.  Serial  ports 
are  available  on  the  multifunction  board  in  the  data  processing 
computer  to  establish  communications  links  with  these  systems,  and 
an  adequate  parallel  interface  capability  was  provided  by  existing 
Apple  II  peripheral  circuit  cards  which  had  seen  use  in  prior  data 
acquisition  systems.  Although  the  operation  of  the  Apple  II 
computers  under  control  of  the  interpreted  Applesoft  BASIC  language 
is  not  particularly  fast,  it  proved  sufficient  to  meet  the  demands 
of  controlling  the  attached  devices  with  fractional  second 
accuracy,  marking  a  significant  improvement  over  purely  manual 


methods. 


The  Apple  II  computer  systems  utilize  the  Synertek  6502 
microprocessor  operating  at  a  clock  speed  of  1.023  MHz,  with  48 
kilobytes  of  random  access  memory.  Integral  video  display 
circuitry  provides  the  24  line  by  40  character  monochrome  text  mode 
used  in  this  application.  Peripheral  cards  attach  to  an  eight  slot 
expansion  bus;  unlike  the  IBM  PC-XT  design  which  permits  addressing 
of  expansion  cards  independently  of  slot  position,  the  Apple  II 
motherboard  contains  address  decoding  circuitry  which  usually 
requires  assignment  of  a  specific  edge  connector  to  each  interface 
card. 

Since  all  recording  of  experimental  data  is  performed  on  the 
8088-based  data  processor,  only  the  software  for  equipment  control 
need  be  resident  on  the  Apple  II  system.  Adequate  disk  storage  is 
provided  by  a  5.25  inch  single-sided  single-density  floppy  disk 
drive  capable  of  recording  143  kilobytes  per  volume.  The  disk 
controller  card  resides  in  slot  6  of  the  Apple  II  expansion  bus. 

Slot  5  contains  a  Prometheus  VERSAcard  multifunction  board, 
which  includes  an  RS-232  serial  port,  parallel  printer  port  and 
real-time  clock/calendar.  Although  it  physically  occupies  a  single 
position  on  the  bus,  the  timekeeping  hardware  is  logically 
addressed  as  slot  1,  and  the  printer  port  as  slot  7;  these 
connectors  become  unusable  for  other  peripheral  cards.  Further, 
slot  0  is  reserved  for  access  to  additional  read-only  memory 
containing  software  enhancements.  As  a  result,  only  slots  2,  3, 
and  4  remain  available  for  installation  of  parallel  interface 
cards,  which  proves  sufficient  for  the  existing  receiver 
configuration  but  severely  limits  future  expansion. 


Design  and  Operation  of  Data  Acquisition 
Interface  Unit 


Three  modes  of  operation  were  envisioned  for  the  diagnostic 
reception  of  ULF,  ELF,  and  VLF  emissions.  In  Mode  1,  the  incoming 
ELF  or  VLF  signal  is  demodulated  by  the  lock-in  analyzer;  the 
resulting  in-phase  and  quadrature  outputs  are  digitized  at  a  rate 
far  below  that  necessary  for  direct  sampling  of  the  baseband 
signal,  and  recorded  on  disk  storage.  Mode  2  provides  high-speed 
digitization  of  the  ELF  or  VLF  signal  prior  to  demodulation, 
requiring  DMA  operation  to  achieve  the  necessary  throughput  of 
values  between  the  data  acquisition  board  and  random  access  memory. 
Mode  3  is  intended  for  reception  at  ULF  frequencies,  which  allow 
digitization  of  both  the  baseband  signal  and  demodulator  outputs  at 
speeds  comparable  to  those  of  Mode  1. 

Data  Collection  and  Recording  Loop 

Digitized  data  is  placed  in  a  data  buffer  as  it  is  collected 
by  the  interrupt-driven  background  task.  Two  separate  areas  are 
allocated  for  alternate  use  to  prevent  lost  information  during  the 
time  required  for  recording.  The  buffers  serve  as  shared  memory 
for  interprocess  communication  between  multiple  tasks,  and  ASYST 
provides  predefined  words  which  permit  examination  of  the  current 
buffer  status  from  the  foreground  program. 

The  main  program  continuously  monitors  the  active  buffer,  in 
order  to  update  the  displayed  percentage  of  data  accumulated. 
Every  10th  set  of  I  and  Q  measurements  is  extracted  from  the  buffer 
as  available  for  calculation  of  the  corresponding  magnitudes  and 
phases;  these  are  plotted  on  the  graphic  display  to  form  a  dynamic 


representation  of  the  incoming  signal.  More  frequent  screen 
updates  are  difficult  to  attain  with  the  computing  power  found  in 
the  present  equipment  configuration.  Further,  a  graphics 
presentation  containing  more  detail  would  likely  be  of  little  value 
when  displayed  with  the  limited  resolution  video  hardware  currently 
in  use.  Plotting  every  tenth  point  has  been  found  to  give 
satisfactory  results  in  actual  experimental  operation. 

The  foreground  task  also  monitors  the  keyboard  for  activity 
during  every  iteration  of  the  acquisition  loop.  If  a  keypress  is 
detected,  then  the  addition  of  points  to  the  magnitude  and  phase 
plots  is  suspended  and  appropriate  action  is  taken.  Upon 
completion  of  the  interactive  command,  the  display  is  updated  as 
quickly  as  possible  in  order  to  correctly  reflect  the  data  resident 
in  the  active  buffer. 

Unfortunately,  ASYST  permits  only  data  acquisition  operations 
to  be  concurrently  multitasked;  all  other  functions  such  as 
operator  communication  through  the  keyboard/display  and  disk  file 
access  must  be  handled  in  a  preemptive  manner  due  to  the  non- 
reentrancy  of  the  disk  operating  system.  Should  too  much  time 
elapse  before  completion  of  command  entry,  it  is  possible  for  both 
buffers  to  become  overwritten,  destroying  data  which  has  been 
collected  but  not  yet  recorded  on  disk. 

When  a  buffer  is  completely  filled  with  data,  the  foreground 
routine  detects  that  a  switch  to  the  alternate  data  area  has 
occurred.  It  then  writes  the  collected  data  to  disk  storage  and 
clears  the  dynamic  display  areas  before  repeating  the  entire 
procedure  on  the  other  buffer.  A  fla"g— variable ,  modified  by  an 


interactive  user  command,  is  tested  after  each  data  file  is 
written;  if  set,  the  acquisition  loop  is  terminated  and  control  is 
returned  to  the  ASYST  interpreter. 

Although  the  transmitter  modulator  itself  contains  waveshaping 
circuitry  intended  to  limit  the  high-frequency  content  of  input 
signals,  care  must  be  taken  to  avoid  rapid  transients  which  can 
cause  damage  to  the  output  stages  of  the  power  amplifiers.  For 
this  reason,  modulation  is  typically  removed  from  the  transmitter 
when  the  synthesizer  frequency  is  manually  changed  from  its  front- 
panel  switches;  remote  synthesizer  control  in  binary  format, 
however,  ensures  amplitude  continuity  and  allows  the  frequency  to 
be  reselected  under  full-power  operating  conditions. 

Fairbanks  Receiver  Site 

The  receiver  system  was  located  in  a  wooded  area  immediately 
adjacent  to  the  Elvey  Building,  which  houses  the  Geophysical 
Institute  on  the  campus  of  the  University  of  Alaska  at  Fairbanks. 
The  electronic  and  computer  equipment  was  set  up  in  a  trailer  which 
housed  support  circuitry  for  the  Institute's  induction 
magnetometers  and  earth  current  sensors.  This  allowed  convenient 
access  to  the  magnetometer  outputs  for  digitization  and  storage  by 
the  receiver's  data  acquisition  system. 

All  loop  antennas  were  placed  approximately  20  m  from  the 
trailer  in  order  to  minimize  interference  pickup  by  the  sensitive 
receivers.  An  orthogonal  layout  was  used  to  eliminate  the  manual 
rotation  performed  in  prior  experimental  campaigns  for  determining 
the  orientation  of  the  modulated  current  system.  For  convenience, 
the  loop  antenna  positions  are  referred  to  as  "A"  and  "B" .  The 


plane  of  the  "A"  loops  is  along  the  north-south  axis  with  maximum 
sensitivity  to  signals  produced  by  currents  flowing  in  an  east-west 
direction.  Conversely,  the  "B"  loops  display  maximum  sensitivity 
to  signals  generated  by  currents  aligned  in  a  north-south 
direction.  These  designations  are  carried  throughout  the  remainder 
of  the  dual-channel  receiver  system.  Figure  23  is  a  photograph  of 
the  air-core  loop  antennas  as  deployed  for  the  October,  1987 
experimental  campaign. 

The  1  m  diameter  loop  antennas  were  used  for  reception  at 
frequencies  between  800  Hz  and  8000  Hz,  while  the  0.5  m  loops  were 
used  at  frequencies  between  200  Hz-  and  1200  Hz.  Extreme  care  was 
taken  to  prevent  inadvertent  phase  reversals  throughout  the 
receiver  system,  with  regard  to  the  winding  sense  of  the  pickup  and 
calibration  coils  as  well  as  the  balanced  transmission  line 
interconnections . 

During  the  October  campaign,  routine  testing  controlled  by  the 
FRQSWP.001  software  revealed  a  drop  in  the  frequency  response  of 
the  large  (1  meter  diameter)  loop  antennas  near  1200  Hz.  This 
anomaly  is  due  to  absorption  in  the  self -resonant  circuit  formed  by 
the  adjacent  small  (0.5  meter  diameter)  loops,  which  were 
disconnected  during  the  tests;  indeed,  the  small  loops  were  known 
to  resonate  at  approximately  1200  Hz  based  on  measurements  made  at 
the  time  of  construction.  Accordingly,  operational  procedure  now 
requires  that  the  small  loop  antennas  be  short-circuited  at  their 
terminal  posts  when  not  connected  to  the  remainder  of  the  system, 
in  order  to  prevent  unsatisfactory  receiver  performance  in  this 
frequency  range. 


Tuning  of  Loop  Antennas 


Once  the  receiver  system  has  been  set  up,  the  loop  antennas 
must  be  tuned  for  resonance  at  the  anticipated  operating 
frequencies.  The  ASYST  program  TUNING. 001  is  used  to  determine  the 
optimum  tuning  unit  capacitance  settings  for  each  loop.  A  strong 
signal  is  coupled  into  the  receiving  loop  from  the  associated 
calibration  coil,  and  the  capacitive  loading  is  progressively 
adjusted  by  computer  control  until  maximum  signal  level  is  detected 
on  the  lock-in  analyzer. 

Upon  completion  of  the  tuning  process,  FRQSWP.001  is  executed 
to  measure  the  bandwidth  and  Q  value  at  each  center  frequency. 
Unlike  the  adaptive  tuning  routine  which  varies  capacitance  while 
maintaining  a  constant  excitation  frequency,  the  response 
measurement  program  sweeps  the  calibration  coil  drive  frequency  in 
order  to  locate  the  -3  dB  points  of  the  amplitude  transfer  function 
at  a  specific  loading. 


Remote  Sensing  of  ELF  Signals  in  a 
Penn  State  Mobile  Receiving  Van  for 
— the  Purposes  of  Evaluating 
Communications  Applications 

One  of  the  problems  to  be  evaluated,  during  the  1990  summer 
heating  campaign  at  the  HIPAS  Facility  in  Alaska,  was  the  ability 
of  the  ELF  signals  created  by  ionospheric  heating  to  be  injected 
into  the  earth- ionosphere  waveguide.  If  distant  communications 
were  to  be  achieved,  it  was  necessary  to  know  if  the  ELF  signals 
were  ducting  through  this  waveguide  or  being  radiated  out  of  the 
heated  region  in  some  other  direction.  ELF  signals  had  to  be 
received  with  reasonable  strength  at  locations  farther  than  the  25 
miles  separating  the  transmitter  from  the  fixed  Penn  State 
receiving  station  located  near  the  Geophysical  Institute  at  the 
University  of  Alaska,  Fairbanks  (U  of  AF) . 

A  single  channel  Penn  State  VLF/ELF/ULF  receiving  system  was 
mounted  in  a  26  foot  mobile  van  and  driven  to  locations  ranging 
from  130  miles  to  387  miles  from  the  heated  volume  in  an  attempt  to 
determine  the  signal  strength  at  each  of  the  locations.  The 
receiver  located  near  the  U  of  AF  was  receiving  the  same  signals 
during  the  experiment  to  allow  for  a  correction  to  the  variation  of 
the  signal  amplitudes  received  at  the  mobile  van  in  order  to 
compensate  the  variations  in  ionospheric  activity.  As  shown  in  the 
attached  map,  all  of  the  sites  were  located  along  a  line  that  is 
nearly  magnetic  south  from  the  HIPAS  facility. 

The  transmitter  at  HIPAS  was  modulated  with  the  ELF  carrier 
frequency  as  in  normal  ELF  generation  experiments.  This  ELF 
modulation  was  further  modified  by  turning  the  modulation  on  and 


off  (Amplitude  Modulation,  AM) ,  Quadrature  Phase  Shift  Keying 
(QPSK)  modulation  _or  CW  (continuously  modulated  with  a  constant 
amplitude  of  the  ELF  frequency) .  Another  method  of  modulating  the 
radiated  ELF  signal  strength  and  phase  was  by  rotating  the  beam  of 
the  HF  signal  in  a  circular  path  tilted  30  degrees  off  vertical. 
This  method  was  coined  "conical  scan"  and  caused  a  change  in  signal 
path  and  the  region  of  the  ionosphere  that  was  being  heated;  thus, 
causing  a  variation  in  phase. and  amplitude  to  be  detected  at  the 
receiver.  Plots  attached  show  a  typical  signal  amplitude  and  phase 
received  during  QPSK  modulation  and  the  next  figure  shows  the 
amplitude  and  phase  variations  caused  by  conical  scanning  the  beam 
as  noted  above.  Both  of  these  plots  are  derived  from  received  data 
taken  at  the  Cantwell,  AK  location. 

The  results  of  the  remote  experiment  indicate  that  the  ELF 
signals  are  being  injected  into  the  earth-ionosphere  waveguide  as 
predicted.  Furthermore,  the  two  station  reception  (PSU  at  U  of  AF 
and  remote  receiver)  allows  for  correction  of  the  radiated  signal 
strength  due  to  ionospheric  variations.  This  correction  coupled 
with  the  measurements  taken  at  various "distances  from  HIPAS  permit 
the  calculation  of  the  attenuation  of  the  ELF  signal  as  it  travels 
along  the  waveguide. 
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The  High  Power  Auroral  Stiaulation  (HIPAS)  heating  facility  haa  been  used  to 
modulate  D  region  ionospheric  currents  at  high  latitudes,  producing  very  low 
frequency  (VLF)  radio  wave  eaisaiona.  The  behavior  of  these  ionospheric  currents  can 
be  deduced  froa  a  coaprehensive  study  of  the  VLF  signals  received  at  a  local  field 
site.  This  paper  examines  the  relationship  between  the  VLF  magnetic  field  strength 
measured  on  the  ground  and  the  intensity  of  an  overhead  electrojet  current  for  the 
purpose  of  enhancing  communications.  The  mapping  of  the  polar  electrojet  current 
from  the  E  region  down  through  the  D  region,  where  it  can  then  be  modulated  by  the 
heater  beam,  is  investigated.  A  finite  difference  solution  to  the  eiectrojet  mapping 
problem  is  presented  in  which  arbitrary  conductivity  profiles  can  be  specified. 
Results  have  been  ootained  using  a  simple  Cowling  model  of  the  electrojet.  These 
results  indicate  that  for  an  eiectrojet  flowing  at  an  altitude  of  110  km  with  a  scale 
size  in  excess  of  ICO  km,  the  capping  of  the  horizontal  current  density  can  be 
completely  characterized  in  terms  of  the  Pedersen  and  Hall  conductivities.  This 
indicates  that  the  mapping  becomes  independent  of  scale  sizes  which  exceed  100  km.  A 
promising  new  diagnostic  technique,  for  studying  ionospheric  D  region  currents,  has 
been  implemented  using  tne  HIPAS  facility.  This  technique  involves  nigh  frequency 
IHF)  beam  steering  for  iocaiized  VLF  generation  in  the  mapped  region  below 
electrojets.  Beam  steering  has  been  used  to  estimate  the  strength  and  current 
distribution  of  the  polar  eiectrojet,  and  for  charting  the  movements  of  overhead 
currents . 


I.  Introduction 

The  generation  of  ELF/VLF  signals  by  modulation  of  the  dynamo  current  system 
using  the  HF  heating  facility  located  near  Arecibo,  Puerto  Rico,  has  been  reported  by 
Ferraro  et  al.  [1]  and  Ferraro  et  al.  [2].  Radiation  from  a  heated  and  modulated 
equatorial  eiectrojet  current  system  was  detected  by  Lunnen  et  al.  (3], 
lonospherically  produced  signals  resulting  from  periodic  plasma  heating  of  the  polar 
eiectrojet  current  system  have  been  studied  by  Stubbe  et  al.  [4J,  Rietveld  et  al. 

(5),  Rietveld  et  al.  (61,  and  more  recently  by  Ferraro  et  al.  (7). 

A  series  of  ionospheric  heating  experiments  have  been  conducted  at  the  HIPAS 
facility  near  Fairbanks,  Alaska,  over  a  period  of  two  years  from  June  1987  to  August 
1989.  The  transmitting  system  consists  of  eight  individual  HF  heating  transmitters 
which  feed  the  eight  antennas  that  comprise  the  HIPAS  heating  array.  The  HIPAS  array 
has  a  total  of  eignt  crossed-dipoles ,  seven  of  which  are  equally  spaced  around  the 
circumference  of  a  340  foot  radius  circle,  with  an  eighth  one  situated  at  the  center. 
Ionospheric  heating  is  achieved  by  modulating  the  HF  continuous  wave  output  of  the 
transmitter . 

A  reconf iguraole  ELF/VLF  receiver  was  positioned  47  km  west  of  HIPAS  at  a 
University  of  Alaska  Geophysical  Institute  field  site  (8).  A  polarimeter  was  set  up 
at  the  receiving  site  in  order  to  monitor  the  polarization  characteristics  of  the 
incoming  ELF/VLF  signals  at  the  ground.  This  polarimeter  consisted  of  two  orthogonal 
loop  antennas,  one  with  its  plane  oriented  in  the  magnetic  north-south  direction  and 
the  other  with  its  plane  in  the  east-west  direction.  The  components  of  the  ELF/VLF 
signals  which  were  intercepted  by  each  of  the  tuned  loops  would  undergo 
preamplification  and  filtering  Lock-in  analyzers  then  extracted  the  in-phase  (I! 
and  quadrature  IQI  components  of  the  signals  using  narrowband  coherent  detection. 

The  outputs  of  the  lock-in  analyzers  were  digitized  by  a  data  acquisition  system  for 
analysis,  display,  ana  storage  via  a  microcomputer.  Highly  stable  crystal  frequency 
standards  were  used  at  both  the  transmitter  and  receiver  locations  to  insure  precise 
frequency  coherence  and  sampling  accuracy. 

Experiments  were  carried  out  in  which  the  HIPAS  heater  beam  was  conically 
scanned  at  an  angle  of  30"  with  respect  to  the  vertical.  The  total  time  to  complete 
one  of  these  scans  was  two  minutes.  Some  conical  scanning  VLF  data  is  presented  and 
interpreted  in  this  paper. 


IX.  The  Downward  Mapping  of  the  Polar  Electrojet  Current 

A  simple  Cowling  oodei  was  used  by  Werner  and  Ferraro  (9)  to  represent  the  polar 
electrojet  current.  The  potential  associated  with  an  electrojet  current  was  found  to 

be 
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The  parameters  1.  and  X,  are  the  spatial  wavelengths  in  the  geomagnetic  east-west  and 
north-south  directions,'  respectively ,  of  an  electrojet  flowing  at  an  altitude  of  z  . 
The  spatial  wavelengths  determine  scale  size  of  an  electrojet  and  are  related  by 
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The  function  Z(z) 
characterized  by 

appearing  in  Eq.  (1)  is  the  solution 
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is  the  effective  wave  number  of  the 
for  a  plasma  in  a  magnetic  field  is 
potential  function  given  in  Eq.  (1), 

electrojet.  The  generalized  form  of  Ohm’s  law 
used  to  relate  the  current  density  to  the 
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13  the  ionospheric  conductivity  tensor.  The  tern  o.  is  known  as  the  direct 
conductivity,  o.  the  Petersen  conductivity,  and  o,  tne  Hall  conductivity.  An 
expression  for  the  normalized  magnitude  of  the  horizontal  current  density  is  given  by 
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and 


is  the  Cowling  conductivity.  The  corresponding  angle  of  the  horizontal  current 
density,  measured  in  degrees  south  of  west,  is  then 
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The  boundary  value  problea  described  by  Eq.  (5),  Eq.  (6),  and  Eq,  <71  aust  be 
solved  numerically,  with  the  exception  of  a  few  special  cases  of  the  conductivities 
(9).  A  finite  difference  scheme  was  used  by  Werner  (10)  to  obtain  a  numerical 
solution  to  this  boundary  value  problem.  The  ability  to  specify  arbitrary 
conductivity  profiles  has  been  incorporated  into  the  nuaerical  mapping  model.  The 
differential  equation  Eq.  (S)  can  be  approximated  at  z*z;  by  the  difference  equation 
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with  the  boundary  conditions 
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Figure  1  shows  some  normalized  horizontal  current  density  profiles  for  the  altitude 
range  between  100  km  and  60  km.  Profiles  for  several  different  effective  wavelengths 
are  included  (dashed  curves),  while  the  solid  curve  represents  the  profile  to  wmch 
the  horizontal  current  density  converges  with  ”5 1 00  km.  This  curve  is 
characterized  by 
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which  only  depends  on  the  Hall  and  Petersen  conductivities.  The  current  density 
profiles  of  Figure  1  were  obtained  assuming  the  presence  of  an  intense  magnetic  storm 
during  sunspot  maximum  daytime  conditions.  An  electrojet  source  height  of  100  km  was 
assumed  in  accordance  with  the  results  of  Kamide  and  Brekke  (11).  The  ionospheric 
conductivity  model  corresponding  to  sunspot  maximum  daytime  conditions  was  adopted 
from  Hughes  and  Southwood  [12]. 

III.  The  Ambient  and  Modulated  D  Region  Current  Densities 

The  VLF  ionospheric  source  can  be  treated  as  an  incremental  volume  element 
A*AyA*  with  a  modulated  horizontal  current  density  of  |AJ^|  centered  at  the  altitude 
of  maximum  heating.  For  typical  experimental  conditions,  the  extent  of  the  effective 
radiating  layer  &z  is  approximately  1  km  and  the  layer  is  assumed  to  be  centered  at 
TO  km.  The  current  density  of  the  modulated  ionosphere  AJ  can  be  related  to  the 
current  density  of  the  ambient  ionosphere  J.  by 
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and 
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For  a  vertical  geomagnetic  field,  the  aabient  (unprimed)  and  the  heated  (priaed) 
conductivity  tensors  are  given  by 
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It  follows  from  Eq.  123),  Eq.  (26),  and  Eq.  (27)  that 


,  :  1  1 

9  9\  -9.9.' 

:  ;  1 

V  a.  +  9:  J 

■  <7*  +  9", 

,  n\9'.  ~  9:9  .  \ 

1  9-9\  +  9:9':\ 

'  -  *;  J 

{  ' 

0 

0 

If  we  let 
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then  the  modulated  horizontal  current  densities  can  be  expressed  in  terms  of  the 
ambient  horizontal  current  densities  in  the  following  way 
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Using  Eq .  (31),  it  i 
current  density  and 


an  be  shown  that  the 
the  total  horizontal 


magnitude  of  the  total 
modulated  current  dens 


horizontal  ambient 
ty  are  related  by 
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The  horizontal,  nodulated  current  denaity  nay  be  eatiaated  froa  the  intensity  of 
the  VLF  magnetic  field  aeaaured  on  the  ground.  The  valuea  of  Ao  are  determined  uaing 
ionoaphenc  heating  theory  [13  1,  (14  1.  The  aabient  horizontal  current  denaity  can 
then  be  calculated  from  Eq.  (32).  Finally,  the  atrength  of  the  electrojet  current  J( 
can  be  found  uaing  the  results  of  napping  theory  Eq.  (22).  Figure  2  shows  a  block- 
diagram  of  the  aodel  developed  to  study  the  relationship  between  the  VLF  magnetic 
field  intensity  measured  on  the  ground  and  the  strength  of  the  polar  electrojet. 

IV.  Estimation  of  Electrojet  Current  Density 

The  Alaska  meridian  chain  of  magnetometers  was  used  to  confirm  the  presence  of 
an  electrojet  over  the  HIPAS  facility  when  heating  experiments  were  in  progress.  The 
Alaska  meridian  chain  consists  of  several  flux-gate  magnetometers  located, 
approximately,  along  a  line  of  constant  geomagnetic  longitude.  Each  magnetometer 
measures  three  components  of  the  earth’s  vector  magnetic  field  in  a  coordinate  system 
in  which  the  H-component  represents  magnetic  north,  the  D-component  represents 
magnetic  east,  and  the  Z-component  points  down  towards  the  ground.  In  general, 
currents  flow  in  the  east-west  direction  in  the  high-latitude  ionosphere. 

Figures  3,  4,  and  S  show  the  H-,  D- ,  and  Z-tracea,  respectively,  for  July  22, 

1988  [15].  The  magnetometer  stations  of  greatest  interest  are  Talkeetna  (TLK)  which 
is  south  of  Fairbanks.  and  Fort  Yukon  1 FYL' )  which  is  north  of  Fairbanks.  At  0335 
local  time  (1135  UT1  on  July  22,  the  HIPAS  HF  heater  beam  was  being  modulated  at  a  5 
kHz  rate  in  the  vertical  position.  The  corresponding  B-field  intensity  measured  on 
the  ground  was  0.68  pT.  The  magnetograms  for  this  period  indicated  that  there  was  a 
negative  perturbation  in  the  H-component  measured  at  the  FYU  and  TLK  stations.  There 
was  no  perturbation  in  the  0-component  measured  at  FYU  and  TLK.  The  Z-component  was 
positive  to  the  north  (FYU)  and  negative  to  the  south  (TLK)  of  Fairbanks.  These 
measurements  indicate  that  there  was  a  westward  traveling  eiectrojet  present  over 
Fairbanks  during  this  time. 

A  VLF  B- f ield  . intenp l ty  of  0.68  pT  corresponds  to  a  horizontal  modulated  current 
density  of  9.25x10'  i./km',  assuming  that  the  altitudes  of  maximum  nesting  is  70  km. 
Papadopoulos  et  ai .  (14)  evaluate  the  level  of  conductivity  modulation  at  various 
ionospheric  heights  as  function  of  incident  HF  power  denaity.  The  values  of  Ao.  and 
Ao,  can  be  found  from  these  curves  assuming  an  altitude  of  70  km  and  a  maximum  HF 
ground  power  for  HIPAS  of  1.6  MW.  These  modulated  conductivity  values  used  in 
conjunction  with  Eq.  (32)  imply  £hat  the  horizontal  ambient  current  density  |J,|  at  a 
70  km  altitude  was  2.56x10*  A/km'.  If  it  is  assumed  that  the  electrojet  was  flowing 
at  an  altitude  of  110  km  and  had  an  effective  wavelength  which  exceeded  100  km, ,then 
Eq.  (22)  suggests  that  the  strength  of  the  electrojet  J,  must  have  been  5.1  A/km'. 

This  electrojet  current  density  is  within  the  range  of  values  reported  by  Kamide  and 
Brekke  [111. 

V.  The  Conical  Scanning  Experiment 

A  promising  new  diagnostic  technique,  for  studying  ionospheric  D  region 
currents,  has  recently  been  implemented  et  the  HIPAS  facility  [16).  This  technique 
permits  a  localized  cross  section  of  the  ionosphere  to  be  probed  by  steering  the  HF 
heater  beam.  The  VLF  signals  resulting  from  this  beam  steering  are  monitored  and 

recorded  at  a  field  site  on  the  ground  via  a  coherent  detection  scheme.  This  beam 

steering  diagnostic  technique  was  successfully  demonstrated  during  a  HIPAS 
ionospheric  heating  campaign  which  took  place  in  July  of  1988. 

Scanning  of  the  HIPAS  heater  beam  in  a  conical  fashion  was  investigated  during 
the  July  1988  campaign.  The  conical  scan  was  primarily  used  for  diagnostics  of  the 
ionospheric  current  system  in  the  vicinity  of  the  HIPAS  facility.  Under  ideal 
conditions  of  a  uniform  overhead  current,  the  signature  of  the  VLF  received  on  the 
ground  resulting  from  a  conical  scan  should  exhibit  a  certain  characteristic  shape. 
Any  deviation  from  tnis  shape  may  indicate  the  movement  of  a  current  into  the  path  of 

the  conically  scanning  heater  beam.  Since  the  position  of  the  heater  beam  is  known, 

it  is  possible  to  use  the  conical  scanning  technique  to  pinpoint  the  geograpnic  lor 
geomagnetic)  location  of  the  influxing  current. 

The  geometry  illustrating  the  conical  scanning  mode  of  operation  is  shown  in 
Figure  6.  The  location  of  the  receiver  site  is  47  km  due  west  of  the  HIPAS  facility, 
i.e.  d=47  km  and  A' =270".  The  coelevation  0  of  the  HIPAS  heater  beam  is  fixed  at  30' 
while  the  azimuth  A  of  the  beam  is  incremented  18*  every  6  seconds  with  geographic 
east  (A*901)  as  the  starting  point.  This  produces  a  conical  scanning  of  the  heater 
beam  which  traces  out  a  circle  in  the  ionosphere  above  HIPAS  once  every  two  minutes. 

An  incremental  volume  of  current  radiating  in  free  space  was  adopted  to  model 
the  VLF  ionospheric  source.  The  incremental  current  volume  was  chosen  because  it  can 
be  used  as  a  basic  unit  to  construct  more  complex  radiating  structures.  The  VLF 
source  model  relates  the  magnetic  field  strength  measured  on  the  ground  to  the 
strength  of  the  modulated  horizontal  current  density  at  the  altitude  of  maximum 
heating.  This  elementary  source  model  was  used  in  the  analysis  of  conical  scan  data. 

Typical  VLF  magnitude  and  phase  data  obtained  from  a  conical  scan  is  shown  in 
Figures  7  and  8,  respectively.  These  measurements  were  made  from  0558-0600  local 
time  (1358-1400  UT )  just  prior  to  an  eiectrojet  event  on  July  21,  1988.  The 
transition  of  the  heater  beam  from  a  vertical  position  to  the  scanning  mode  is 
visible  during  the  first  few  seconds  of  the  conical  scan  data  and  should  be 
disregarded.  The  frequency  of  VLF  signals  generated  during  the  conical  scanning 
experiment  was  5  kHz. 

The  dashed  curve  appearing  in  Figure  7  represents  the  normalized  field  strength 
that  would  result  from  a  uniform  ionospheric  current  flowing  at  an  altitude  of  70  km. 


The  magnitude  data  is  not  exactly  symmetrical  indicating  that  heating  by  sideiobee 
may  produce  VLF  radiation  which  interferes  with  the  radiation  generated  by  heating 
froa  the  aain  beaa.  Figure  9  is  a  current  map  showing  the  relative  intensity  of  the 
ionospheric  currents  associated  with  the  E-W  conical  scan  aagnitude  data  of  Figure  7. 
The  length  of  each  line  segment  in  the  current  map  is  proportional  to  the  ratio  of 
the  measured  cSTTIcal  scan  data  to  the  theoretical  values  for  a  uniform  current.  If 
'he  length  or'  the  line  segment  is  less  than  one  unit,  then  the  ionospnenc  current  is 
weaker  than  that  predicted  by  the  uniform  current  model.  And  if  the  length  of  the 
line  segment  is  greater  than  one  unit,  then  the  ionospheric  current  is  stronger  than 
that  predicted  by  the  uniform  current  model.  The  coordinates  used  for  the  current 
map  are  geomagnetic.  Hence  the  current  is  strongest  to  the  north-west  of  HIPAS, 
which  is  in  agreement  with  magnetometer  records  of  this  period. 

The  results  presented  above  suggest  that  conical  scanning  could  be  used  for 
communications  purposes  as  well  as  diagnostics.  The  heater  beam  could  be  adaptively 
steered  to  the  region  of  the  ionosphere  where  the  strongest  current  was  detected  by  a 
conical  scan.  This  procedure  could  be  repeated  periodically  to  account  for  any 
changes  in  the  location  of  the  most  intense  current  brought  about  by  the  dynamic 
polar  ionosphere. 

Figures  10  and  11  show  the  VLF  magnitude  and  phase  data  resulting  froa  a  conical 
scan  made  during  0434-0436  local  time  (1234-1236  UT)  on  July  21,  1988.  The  aagnitude 
data  for  this  conical  scan  reveals  that  the  strongest  current  was  located  to  the 
geographic  north-west  of  HIPAS.  If  this  were  an  adaptive  conical  scan  the  heater 
bean  would  subsequently  be  parsed  at  a  coelevation  of  30s  and  an  azimuth  of  324’. 

The  current  map  corresponding  to  the  E-W  magnitude  of  this  conical  scan  is  shown  in 
Figure  12. 

VI.  Conclusion 

The  relationship  between  the  VLF  magnetic  field  strength  measured  on  the  ground 
and  the  intensity  of  an  overheao  electrojet  current  has  been  treated  in  this  paper. 

The  currents  that  are  being  modulated  in  the  D  region  map  down  from  the  E  region 
where  the  electrojet  source  current  flows.  A  finite  difference  scheme  was 
implemented  to  find  a  numerical  solution  to  the  electrojet  mapping  boundary  value 
problem.  This  numerical  mapping  ooael  accommodates  arbitrary  conductivity  profiles. 

The  VLF  source  was  modeled  as  an  incremental  volume  of  current  radiating  in  freb 
space.  This  source  model  relates  the  magnetic  field  strength  measurea  by  a  ground- 
based  receiver  to  the  strength  of  the  modulated  horizontal  current  density  at  the 
altitude  of  maximum  heating.  The  associated  ambient  (DC)  current  was  then 
determined.  Finally,  using  this  information,  the  strength  of  the  polar  electrojet 
could  be  estimated.  The  results  of  this  analysis  agreed  favorably  with  values  of 
electrojet  current  densities  reported  elsewhere  in  the  literature. 

The  use  of  the  HF  beam  steering  VLF  generation  technique  as  a  diagnostic  tool 
was  introduced.  Results  were  presented  from  beam  steering  experiments  performed 
during  an  electrojet  event  which  occurred  on  July  21,  1988.  Conical  scanning  data 
was  used  to  infer  ionospheric  D  region  currents. 
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Figure  Captions 


Figure  1. 

Figure  2. 
Figure  3. 

Figure  4 . 

Figure  5. 

Figure  6. 
Figure  7. 

Figure  8. 

Figure  9. 


Normalized  horizontal  current  density  profiles  for  the 
altitude  range  between  100  km  and  60  km. 


Block-diagram  of  the  model. 


Magnetometer  records  for  July  22,  1988.  Alaska  RGON 
magnetometer  chain  H-traces. 


Magnetometer  records  for  July  22,  1988.  Alaska  RGON 
magnetometer  chain  D-traces. 


Magnetometer  records  for  July  22,  1988.  Alaska  RGON 
magnetometer  chain  Z-traces. 


Geometry  for  the  conical  scanning  experiment. 


Conical  scan  VLF  magnitude  data  for  the  period  1358-1400 
UT  on  July  21,  1988. 


Conical  scan  VLF  phase  data  for  the  period  1358-1400  UT 
July  21,  1988. 


The  relative  intensity  of  the  E-w  component  of  the  D 
region  ionospheric  currents  associated  with  the  E-W  VLF 
magnitude  data  of  Figure  7. 


Figure  10.  Conical  scan  VLF  magnitude  data  for  the  period  1234-1236 
UT  on  July  21,  1988. 


Figure  11.  Conical  scan  VLF  phase  data  for  the  period  1234-1236  UT 
on  July  21,  1988. 


Figure  12.  The  relative  intensity  of  the  E-W  component  of  the  D 
region  ionospheric  currents  associated  with  the  E-W  VLF 
magnitude  data  of  Figure  10. 
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PROJECT  REVIEW 


During  the  period  fom  9/15/86  to  9/30/90,  UCLA'?  Plasma  Physics 
Laboratory  was  under  subcontract  to  The  Pennsylvania  State  University,  to 
collaborate  with  Professor  A.J.  Ferraro  to  conduct  research  on  polar  VLF 
excitation  and  generation.  The  whole  program  was  under  the  sponsorship  of 
the  U.S.  Congress  Universities  Research  Initiative  <URI),  which  was 
administered  by  the  U.S. Navy's  Office  of  Research,  Washington  D.C.  UCLA’s 
contributions  to  the  UR I  initiative  were  its  very  high  power  radio  frequency 
array  3S  miles  east  of  Fairbanks,  Alaska;  called  HIFAS,  which  stands  for  High 
Power  Auroral  Stimulation  and  UCLA's  ionosonde  diagnostics. 

Four  separate  campaigns  were  conducted  over  the  course  of  the 
URI/PennStat e  subcontract;  namely, 


- 1 

June 

15- 

-25, 

1987 

II 

October 

13-23,  1987 

III 

July 

11 

-22, 

1988 

IV 

July 

12- 

-  Aug 

8,  1989 

The  four  campaigns  are  summarized  m  Table  I,  which  emphasizes  the 
operation  of  the  HIPAS  radio  frequency  heater;  naroeiv,  the  campaign  dates, 
the  durations  of  uninterrupted  operation,  the  total  power  radiated  <le 
8xl00KW  means  8  transmitters  at  iOOKW  each  or  a  total  radiated  power  of 
800KW) ,  and  comments  on  the  over  all  operation  during  the  specified  period  of 
operation. 

Prior  to  the  first  campaign,  the  total  radiated  power  was  brought  up  to 
800  KW  at  3.345  MHz  under  the  advice  of  John  Carroll  (the  orgmal  Platteville 
chief  engineer)  who  came  and  worked  at  the  site.  The  sytem  was  next  retuned 
to  2.85  MHz,  which  is  close  to  the  electron  cyclotron  frequency,  and  which 
should  efficiently  couple  to  the  electrojet  to  generate  VLF  and  ELF.  This 
meant  retuning  both  the  transmitters,  baluns,  arid  antennas.  The  antennas 
were  tuned  by  the  addition  of  "end-straps"  wThich  were  a  pair  of  cables 
attached  to  the  ends  of  each  antenna  and  spread  in  angle  and  anchored  to  the 
ground  to  control  both  the  impedance  and  frequency  of  a  dipole.  This  method 
of  antenna  tuning  u’as  the  contribution  of  Mr.  William  Harrison,  a  consultant 
to  HIPAS.  The  baluns  were  tuned  by  replacing  the  coaxial  cable  that  acted  as 
capacitive  ends  loads  by  high  voltage  variable  vacucr  capacitors  which  could 
be  continously  tuned.  The  tops  of  the  baluns  were  covered  bv  large  plastic 
boxes  to  protect  the  antenna  connections  from  the  elements,  principally  ice 
during  the  u’inter,  hut  rain  during  the  summer.  The  coaxial  feed  lines  were 
also  pressurized  to  a  few  inches  gauge  to  thwart  any  intake  of  mosture  laden 
air  into  the  coaxial  lines  and  baluns,  with  resultant  build  up  of  ice  during 
the  winter.  Initially,  this  air  was  dried  with  commercial  "Dri-  Rite". 

Later  a  refrigeration  type  of  drying  sytem  was  added.  The  pressurization  of 
the  transmission  lines  and  baluns  vTith  dry  air  soon  showed  that  the  system 
could  be  operated  during  the  winter  at  temperatures  as  low  as  -40°F. 


TABLE  I 

TRANSMITTER  PERFORMANCE  SUMMARY 


CAMPAIGN 


I 


II 


DATES 

TIMES 

DURATIONS 

POWERS 

REMARKS 

(hours) 

6/15/87 

0600-0900 

3 

8xl00KW 

10°C  ambient 

(Monday) 

1530-1830 

3 

ti 

rain 

6/16/87 

0700-1030 

3.5 

8xl00KW 

1800-2100 

3 

23°C 

6/17/87 

0400-1600 

12 

8xl00KW 

Coolant 

problems 

6/18/87 

0300-1100 

8 

8xlOOKW 

1700-2300 

6 

22°C 

6/19/87 

0600-1200 

6 

8xl00KW 

6/22/87 

1000-1230 

2.5 

8x10 OKU 

stopped  due  to 

(Monday) 

2030-2400 

3.5 

arc  in  #0 

6/23/87 

1230-1600 

3 .  S 

1900-2200 

3 

6/24/90 

0800-1424 

6.4 

8x10 OKU 

#6  out,  insulator 

1424-1500 

0.6 

6xl00KW 

failure,  #7  off  for 
balance 

1500-1524 

0.4 

8xl00KW 

repaired 

1S24-1600 

1.0 

4xl00KW 

&  4x  70KW  generator 

overheated 

6/25/90 

0630-1100 

4.5 

8x10 OKU 

2200-2400 

2 

8xl00KW 

trips  due  to  air 
in  water  interlocks 

TOTAL  hours 

73.1 

10/13/87 

1300-1900 

5 

CO 

X 

o 

o 

T—’ 

10/14/87 

1300-1500 

2 

ft 

1500-1600 

1 

8x  4 OKU 

fly  over 

1620-2000 

3.3 

8x10 OKU 

10/15/87 

0800-1400 

6 

8xioo:*:u 

1°C  ambient 

10/16/87 

0000-  600 

6 

8X100KW 

10/19/87 

1530-2200 

6.5 

SxlOO.KW 

10/20/87 

0830-1200 

- 

- 

interlock  problems 

1330-1730 

4 

8xl00KW 

10/21/87 

1530-1830 

3 

8x10 OKU 

2000-2400 

4 

7x10 OKU 

#6  off  line 

10/22/87 

1500-2300 

8 

8xl00KW 

10/23/87 

1900-2300 

4 

8X100KU 

TOTAL  hours 

50.3 

(continued) 
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TABLE  I 
(CONTINUED) 


CAMPAIGN 

II 


DATE 

TIMES 

DURATION 

POWERS 

( hours ) 

7/11/88 

1100-1600 

5 

8xl00KU 

( Monday 

)  1900-2200 

3 

If 

7/12/88 

1800-2230 

4.5 

If 

7/13/88 

1100-1200 

1 

7xl00KW 

1200-1400 

o 

i. 

8x10 0KW 

2000-2400 

4 

7/14/88 

1800-2400 

6 

ll 

7/15/88 

2040- 

6.3 

ii 

7/16/88 

0302 

7/17/88 

2100- 

5.2 

M 

7/18/88 

0220 

7/18/88 

2100- 

6 

»• 

7/19/88 

0300 

7/19/88 

1800-2400 

6 

ll 

(Tuesday ) 

7/20/88 

0900-1500 

6 

M 

7/21/88 

0300-0500 

2 

II 

0545-9000 

3.2 

1 1 

7/22/88 

0000-0600 

6 

II 

TOTAL  hours 

65.2 

REMARKS 

#5  off  IS  minutes 
80°F  ambient 
Power  out  at  UA 
#8  off  (meter 
shunt  bad) 

77°F  thunderstorm 

86°F  thunderstorm 
lightening 
beautiful  sunset 
at  0019 


85°F  ambient 
2  frequency  mod 
2. 85  8,2.8525  MHz 

oil  hose  broke  on 
diesel  generator 
52°F  ambient, 
rain  at.  0432 


< continued) 
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TABLE  I 
(CONTINUED) 


CAMPAIGN  DATES 

TIMES 

DURATIONS 

POWERS 

REMARKS 

( hours ) 

IV  7/12/89 

0000-0600 

6 

8xl00KW 

rotating  and 

7/13/89 

0000-0345 

3.75 

If 

scanning  beams 

0345-0400 

0.25 

6x10 OKU 

0  &  7  off  due  to 

0400-0630 

2.5 

8xl00KW 

interlock  problem 

7/14/89 

0000-0300 

3 

It 

64°F  ambient 

7/14/89 

2100- 

5 

ii 

90°F  ambient 

7/15/89 

0200 

7/15/89 

2000-2230 

2.5 

ii 

74°F  ambient 

2230-2300 

0.5 

6xl00KW 

0  8,  7  off  line 

7/16/89 

2300-0100 

2 

8xl00KW 

7/17/89 

2100-2400 

3 

II 

63°F 

0000-0245 

2.75 

8x1 2 OKU 

7/18/89 

0245-0300 

0.25 

8x10 OKU 

74°F 

7/18/89 

2140- 

4.75 

li 

delayed  by  #4 

7/19/88 

0230 

FCC  Anchorage 

called 

7/19/89 

2130- 

4.5 

M 

#6  interlock 

7/20/89 

0200 

delay 

7/20/89 

1900-2330 

4.5 

Ii 

55°F 

7/21/89 

2330-0300 

3.5 

8x1 2 OKU 

78°F 

7/21/89 

1900-0300 

8 

II 

7/31/89 

2100- 

6 

8x10 OKU 

65°F  ambient 

8/01/89 

0300 

8/01/89 

2100- 

6 

ii 

"painting" 

8/02/89 

0300 

8/02/89 

2100- 

6 

it 

66°F  ambient 

8/03/89 

0300 

8/03/89 

2100- 

7 

1 1 

8/04/89 

0400 

8/04/89 

2100- 

6 

1 1 

59°F  ambient 

8/05/89 

0300 

8/06/89 

2100- 

6 

li 

78°F  ambient 

8/07/89 

0300 

8/07/89 

2100- 

6 

1 1 

75°F  ambient 

8/08/89 

0300 

8/08/89 

2000-2130 

1.5 

7x1 0 OKU 

#1  off  line 

2130-0200 

4.5 

8xl00KU 

Finis 

TOTAL  hours 

104.5 

The  second  diesel 

generator 

was  brought  on- 

line  and 

the  two  diesels  were 
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actually  run  in  parallel  (ie  synchronized).  Although  interesting  as  an 
exercise  in  power  generation,  parallel  operation  was  abandoned  in  favor  of 
the  simpler  arrangment  in  which  one  diesel  drove  four  transmitters  while  the 
other  drove  the  other  four.  This  arrangment  avoided  the  problem  and  dangers 
of  loss  of  generator  synchronization.  With  the  two  diesels  operational,  all 
eight  transmitters  could  be  run  at  100KW  each  for  long  durations. 

The  original  Henry  Radio  radio  frequency  amplifiers  that  drove  the 
Platteville  transmitters  were  replaced  by  ENI  (Model  240L,  50  Watt)  broadband 
units,  which  were  much  more  reliable. 

During  the  first  campaign,  73  hours  of  operation  were  delivered.  On  one 
day,  (6/17/87)  the  system  suported  the  Penn  State  experiments  continously  for 
12  hours.  Occationally ,  transmitters  would  «-top  when  interlocks  over 
reacted,  but  generally  the  transmitters  could  put  were  back  on-line  within 
minutes. 

The  HIPAS  heater  is  described  in  a  paper  that  has  been  accepted  for 
publication  in  Radio  Science  m  1990.  It  is  attached  as  Appendix  I,  and  is  a 
very  complete  description  of  the  HIPAS  heater  during  all  of  the  campaigns. 
Noteworthy  of  the  array  is  the  fact  that  each  antenna  is  driven  from  a  single 
Platteville  transmitter  and  each  can  be  individually  controlled  in  phase, 
permitting  the  RF  beam  pattern  to  be  scanned  or  pointed  or  defocused. 

Another  feature  of  the  array  are  the  crossed  radiating  dipoles  which  produce 
either  right  or  left  circular  polarized  patterns  depending  on  how  the 
antennas  are  connected  to  their  balun  terminals.  The  circular  polarization 
of  the  radiating  beam  is  described  in  considerable  detail  the  Appendix  I 
paper . 


At  the  onset  of  the  first  campaign,  any  polarization  change  had  to  be  made 
Ip'  climbing  the  antenna  towers  and  changing  the  antenna  connections  by  hand, 
which  meant  that  half  a  day  was  typically  needed  to  change  the  polarization 
of  the  the  array.  Fy  the  time  of  the  third  campaign  switches  had  been 
installed  on  the  tops  of  each  antenna,  so  that  the  change  could  be  made  from 
“he  ground.  The  polarization  ran  now  be  changed  in  less  than  five  minutes. 

The  polarization  switches  were  of  HIPAS  design  and  were  fabricated  both  at 
-he  site  arid  with  the  aid  of  the*  I'CLA  Physics  Department  machine  shop. 

During  the  first  two  campaigns  the  array  was  phased  or  pointed  manually, 
via  delay  boxes  with  switches  that  set  the  delays  on  the  inputs  to  the 
individual  Platteville  amplifiers  and  associated  antennas.  The  beam  was 
pointed  by  pre  calculation  of  the  array  pattern.  Two  operators,  working 
furiously,  could  reset  the  direction  of  the  beam  in  about  2  minutes.  By  the 
time  of  the  third  campaign,  voltage  variable  phase  shifters  had  been 
installed  and  the  pointing  oc  dephasing  was  controlled  by  a  personel  computer 
i PC) .  A  beam  forming  time  of  10  milliseconds  was  achieved  with  this  first  PC 

controlled  svstem.  The  svstem  was  as  shown  in  Figure  4  of  the  Appendix  I 
paper.  At  the  onset  of  the  last  campaign  an  even  faster  system  had  been 

installed  with  15  microscond  beamforming  times.  The  new  svtem  permitted 
"painting"  experiments. 


To  phase  the  array  properly  each  antenna  lias  to  be  sampled  and  compared 
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with  the  others.  The  base  of  each  antenna  (see  also  Figure  4,  Appendix  I) 
has  a  sampling  loop  which  is  now  connected  to  a  1000  foot  long  piece  of  RG-9U 
coaxial  cable.  The  cables  ends  are  all  brought  to  a  distribution  panel  at 
the  control  console;  there  the  phases  are  compared  with  either  an 
oscilloscope  or  vector  voltmeter.  The  excess  cable  on  each  run  are  left 
wound  on  the  cable  spools  to  serve  as  chokes  to  keep  RF  induced  on  the 
outerbraid  from  interfering  with  the  phase  measurments  (on  the  center 
conductors).  The  longer  cables  to  the  most  distant  antennas  <ie,  #’s  1,  2, 

6,  &  7)are  finished  by  by  winding  ends  on  magnetic  cores.  The  phasing  cables 
were  orginally  RG-58U  coax  which  was  too  light  weight.  Many  were  broken,  RF 
burned,  and  poorly  spliced.  They  were  all  replaced  with  the  forementioned 
RG-9V  cables  with  end  chokes. 

During  the  first  campaign,  harmonics  of  the  radiating  HF  were  found  to 
interfere  with  television  reception  of  channel  2  within  the  local  community 
of  Two  Rivers,  Alaska.  The  problem  was  solved  by  series  resonant  "traps" 
across  the  main  tank  circuit  of  the  Platteville  transmitters.  The  traps  were 
tuned  while  monitoring  televison  reception  at  the  bunk  house,  one  half  mile 
away.  There  was  no  problem  in  subsequent  campaigns 


Ue  also  procured  a  broadband  isotropic  RF  radiation  safety  monitor  ( ie 
NARDA  Model  8616  with  8662R  Efield  ( 0 . 3-1000MHz )  and  8652  H  field  (0.3-30 
MHz)  probes).  The  monitor  was  used  to  carefully  survey  the  site  for  unsafe 
levels  of  radio  freuqency  radiation,  especially  when  the  array  was  radiating 
a  the  megawatt  level.  At  3  MHz  the  US  standard  for  whole  body  radiation, 
averaged  over  0.1  hours,  is  100  mW/cm^.  This  corresponds  to  approximately  500 
volts/meter  amplitude.  This  standard  was  not  exceeded  beyond  100  feet  of  any 
antenna  when  the  array  was  radiating  at  a  megawatt.  During  operation,  access 
roads  into  the  antenna  field  are  barricaded,  a  flashing  light  on  the  top  of 
the  generator  building  is  activated,  and  the  antenna  field  surveyed  from 
control  room  for  inadvertent  intrusion. 

Between  the  third  and  fourth  campaign,  Dr.  Michael  McCarrick  of  this 
laboratory  moved  with  his  family  to  the  HIPAS  site,  for  the  purposes  of 
bette’-  cooordinating  the  scientific  work  at  the  facility.  It  was  during  his 
tenure  at  HIPAS  that  diagnostics  were  set  up  at  the  NOAA,  Gilmore  Creek, 
approximately  35  kilometers  North  West  of  HIPAS,  the  fast  beam  switching 
system!  was  installed,  and  the  array  was  completely  controlled  bv  computers. 
With  an  lonosonde  at  NOAA,  lonograms  could  be  recorded  while  the  heater  was 
operating.  During  this  same  period  special  ELF  detectors,  developed  by  Dr. 
Dave  bent. man  of  ICLA  Space  Sciences  Group  under  a  contract  with  the  US  Air 
Force,  were  installed  at  the  NOAA  site.  With  these  special  detectors,  we 
were  able  to  detect  heater  induced  ULF  excitation  at  11  Hz,  between  the  first 
two  (7  &  14  Hz)  Schuman  earth  E&M  resonances. 


Data  links  were  also  set  up  between  the  PSU  receiver  site  at  the 
lniversitv  of  Alaska  arid  the  HIPAS  transmitter  control  ,  which  allowed 
instantaneous  monitoring  of  t hie  VLF  resonance  at  the  transmitter  control 
center,  and  allowed  an  experimenter  to  see  cause  and  effect  in  an  on-line 
fashion.  This  was  very  important  during  the  beam  steering  experiments.  The 
same  data  line  via  IBM  PC  was  also  established  between  NOAA  Gilmore  Creek  and 
the  HIPAS  control  room  allowing  us  to  monitor  the  ionospheric  conditions 
while  the  heater  was  running  continuously. 
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The  electric  starting  motors  on  the  diesel  generators  were  replaced  by 
air  starters,  which  are  much  more  reliable;  particularly  since  replacement 
electric  motors  are  no  longer  available. 

The  reliablity  of  HIPAS  depends  on  regular  maintenance,  especially  during 
severe  winter  periods.  The  diesel  generators  are  checked  on  a  weekly  basis 
and  all  the  heating  and  cooling  systems  now  have  back  ups.  HIPAS  weathered 
two  severe  winters  <one  cohere  site  temperatures  dropped  to  -70°F)  and  was 
able  to  run  campaigns  after  each.  Ue  attributed  this  performance  to  the 
excellant  on-site  staff  mebers,  who  were  able  to  foresee  potential  problems 
and  take  action. 

Reference  to  Table  I,  Campaign  IV  includes  long-duration  runs  at  total 
radiated  power  of  one  megawatt.  The  transmitters  ran  for  a  total  of  104. S 
hours  with  less  than  1~27  down  times. 

Appendix  II  contains  a  copy  of  a  paper  written  by  R.G. Brandt  for 
presentation  at  a  special  conference  (May  1990)  in  Tromso,  Norway  on 
lonosphereic  modification.  It  is  an  apt  summary  of  the  HIPAS's  contribution 
to  the  four  Fenn  State  VRI  campaigns.  D; .  R.G. Brandt  is  the  Program  Manager 
of  the  Office  of  Naval  Research's  which  supported  the  entire  project. 

SUMMARY 

In  summary,  the  PSU  and  HIPAS  staffs  worked  very  well  together  during  the 
four  URI  campaigns.  HIPAS  scientists  would  have  liked  to  have  had  greater 
access  to  the  data  gathered  during  the  campaigns,  as  well  as  participated  in 
more  post  campaign  discussions,  as  was  orginally  agreed  upon  the  onset  of 
the  program.  In  reptrospect,  the  "opportunity"  mode  (in  which  experiments 
are  conducted  when  conditions  are  favorable)  would  have  been  more  preferable 
to  the  "campaign"  mode  (in  which  experiments  are  conducted  during  a  block  of 
time  regardless  of  conditions),  particularly  due  to  the  fact,  that  the 
ionosphere  is  so  variable. 

During  the  URI  program  the  HIPAS  facility  has  been  brought  to  5  high 
degree  of  operational  reliablity.  It  can  now  be  operated  any  time  of  the 
vear.  In  addition  the  facility  has  been  instrumented  to  a  stage  that  it  can 
be  controlled  locally  or  remotely  by  computers.  The  invention  of  the  "phase- 
modulation"  technique  by  HIPAS  sceintists  made  it  possible  to  modulate  the  HF 
radiation  at  any  low  frequency  without  affecting  the  diesel  power  sources. 
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APPENDIX  I 


Paper  submitted  for  publication  in  Radio  Science  which  descrobes  the 
HIPAS  radio  frquency  heater  essentually  as  it  was  during  the  four 
PennState  URI  campaigns: 

A.Y.Wong.et  all.  High  Power  Radiating  Facility  at  the  HIPAS  observatory, 
Radio  Science,  accepted  for  publication. 
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ABSTRACT 

UCLA's  radio  frequency  ionospheric  heater,  25  miles  east  of  Fairbanks, 
Alaska,  is  described.  The  heater  consists  of  eight  crossed  dipole  antennas 
arranged  in  a  circular  pattern  to  give  a  gain  of  18.4  db  over  isotropic  at 
2.85  MHz  (-  2nd  electron  cyclotron  harmonic).  At  1.2  MW  total  radiated  power 
the  array  has  a  calculated  ERP  of  84  MW.  Each  antenna  is  driven  by  a  150  KW 
transmitter,  originally  from  the  Platteville  heater.  The  eight  transmitter- 
antennas  are  managed  by  a  personal  computer  which  controls  power,  modulation 
and  beam  steering.  Methods  of  tuning  the  antennas,  to  achieve  either  right 
(O-mode)  or  left  (X-mode)  circular  polarized  radiated  beams,  are  described. 
The  heater  is  powered  by  two  1500  HP  diesel  electric  generators.  It  can  be 
operated  throughout  the  year  over  -30°C  to  40°C  ambient  temperature  extremes. 
Future  improvements  include  the  construction  of  an  even  larger  840  x  840 
meter  (24  x  24  antenna)  array  with  a  gain  of  37.5  dbi. 
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1.  INTRODUCTION 


The  HIPAS  Observatory  was  developed  as  an  outdoor  laboratory  for 
research  on  nonlinear  space  plasma  physics  using  high  power  radio  frequency 
(RF)  electromagnetic  waves  to  actively  stimulate  the  auroral  ionosphere, 
hence  the  acronym  HIPAS  (High  Power  Auroral  Stimulation). 

The  Observatory  was  begun  in  1981,  [A.  Wong,  19811.  The  site  has  the 
feature  that  (like  the  one  in  Tromso,  Norway)  the  RF  beam  can  be  launched 
along  the  earth's  magnetic  field  [Stubbe,  19821.  It  also  passes  under  the 
auroral  oval,  resulting  in  spectacular  auroral  displays  throughout  the  late 
fall  to  early  spring  months  [Elvey,  19571.  It  is  accessable,  over  good 
roads,  to  Fairbanks  (2S  miles  to  the  west),  the  University  of  Alaska  (UA), 
the  UA's  Geophysical  Institute,  UA's  Poker  Flat  rocket  launching  site,  the 
NOAA  Gilmore  Creek  Satellite  Tracking  Station,  Eielson  AFB,  and  Fort 
Wainwright . 

The  heart  of  the  system  is  a  circular  eight  antenna  radio  frequency 
array  of  84  MW  Equivalent  Radiating  Power  (ERP),  when  the  transmitter 
outputs  are  1.2  MW  total  (i_e.  8  x  150  KW).  The  array  is  summarized  in 
Table  I. 

The  first  column  of  Table  I  lists  the  frequencies  at  which  the  system 
has  been  operated.  The  system  is  not  frequency  agile.  The  antennas  have 
to  be  retuned  whenever  the  operating  frequency  is  changed  more  than  10  KHz. 
The  present  operating  frequency  is  2.85  MHz  which  is  close  to  the  second 
harmonic  of  the  electron  cyclotron  frequency  about  the  earth's  (~  1/2 
Gauss)  magnetic  field,  and  also  one  of  our  licensed  frequencies. 

A  top  view  schematic  of  the  array  is  given  in  Figure  1.  The  array 
consists  of  seven  antennas  equally  spaced  around  a  circle  of  208  meters 
diameter  with  the  eighth  at  the  center.  Each  antenna  is  connected  to  its  own 
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ISO  KW  (max)  amplifier  through  seven  Inch  outside  diameter  (2-1/2  Inch  center 
conductor),  dry  air  pressurised  <1-2  psi  gauge)  coax  made  from  agricultural 
Irrigation  pipe.  The  Individual  final  amplifiers  and  coaxial  lines  are  all 
modifications  from  the  now  legendary  Platteville  heater,  which  ceased 
operation  in  the  late  1970's  [Carroll,  19741.  This  equipment  was  loaned  to 
UCLA,  which  reconfigured  each  transmitter  with  new  individual  power  supplies. 
The  eight  transmitters  are  now  housed  in  a  building  on  the  edge  of  the 
antenna  field,  as  indicated  in  Figure  1.  Close  by  is  a  transformer  yard  and 
a  large  barn-like  building  that  houses  the  two  1500  horsepower  (or  1.2  MW  each) 
diesel  electric  generators  which  power  the  system. 

As  shown  in  Table  I,  at  2. 85  Hz,  the  array  has  a  calculated  gain  of  18.4 
dbi.  At  a  total  radiated  power  of  1.2  MW,  this  corresponds  to  an  ERP  of  84 
MW,  which  means  that  at  a  range  (R)  of  100  Km,  a  RF  flux  (F  =  ERP/4tiR2)  of 
0.67  mW/m2,  or  a  peak  vacuum  RF  electric  field  amplitude  of  0.7  V/m 
[calculated  from  the  flux  and  the  chacteristic  impedance  of  free  space  <ZD  = 
377  ohms);  namely,  EpeakRF  =  ^2Z0F  ]  - 

A  recent  (1988)  photograph  of  one  of  the  antennas  (i.e.,  #4  in  Figure  1) 
is  shown  in  Figure  2.  The  crossed  dipoles  are  46  feet  above  the  ground. 

Seen  also  in  this  picture  are  the  balun,  the  coaxial  feed  line,  and  the  "end¬ 
loading-straps"  that  are  used  to  retune  the  antenna(s)  to  frequencies  in  the 
range  of  2.8  -  4.9  MHz.  A  second  antenna  is  in  the  background.  The  picture 
also  gives  a  sense  of  the  environment  at  the  HIPAS  site  in  the  spring  (April). 
During  the  winter,  the  antenna  field  is  completely  snow-covered,  but  usually 
not  more  than  a  few  feet  deep. 

By  phasing  the  north-south  dipole  currents  90  degrees  with  respect  to  the 
east-west  dipole  currents,  the  array  can  be  made  to  radiate  either  right  (O 
mode)  or  left  (X  mode)  circularly  polarized  patterns,  depending  on  the 
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connection  of  the  east-vest  dipoles  to  the  balun,  as  shown  schematically  In 
Figure  3.  When  the  antennas  are  connected  to  radiate  a  left  hand  circularly 
polarized  beam,  the  radiated  electric  field  rotates  in  the  same  direction  as 
free  electrons  in  the  ionosphere  about  the  earth's  magnetic  field.  This  is 
called  the  X-mode.  The  other,  right  hand  circular,  polarization  is  called 
the  O-mode.  It  propagates  through  the  ionosphere  with  less  loss. 

The  whole  system  now  <1988-89)  operates  rather  routinely  at  a  total 
radiated  power  of  one  megawatt  CW  throughout  the  year,  over  -30°C  (winter)  to 
+40°C  (summer)  ambient  temperature  extremes.  The  transmitters  can  be 
amplitude  modulated  from  0.15  -  10  KHz.  The  antennas  can  also  be  rephased 
for  the  purpose  of  modulating  the  far  field  pattern  at  ULF  frequencies,  or 
for  steering  the  beam  30°  off  the  vertical  at  rates  up  to  100  Hz. 

Details  about  the  contructlon  of  the  present  RF  ionospheric  heater  system 
with  its  eight  Platteville  transmitters,  new  antennas,  diesel  power  sources, 
arctic  water  cooling  system,  and  original  pulsed  transmitter  are  included 
below. 

II.  RADIO  FREQUENCY  SYSTEM 

A.  CW  TRANSMITTERS 

Figure  4  shows  the  whole  HIPAS  RF  system  schematically.  One  can  start 
at  the  triaxial  balun  at  the  right  side  of  the  page.  The  antennas,  as  seen 
earlier  in  Figure  3,  are  connected  to  the  inner  and  mid  balun  conductors. 

The  7  inch  agricultural  coaxial  lines  connect  to  the  bottom  of  the  balun 
through  a  coaxial  right  angle  reducer.  The  coaxial  lines  (see  Figure  1)  run 
to  the  transmitter  building  where  they  connect  to  coaxial  power  meters  (Bird 
Company)  and  then  to  the  Platteville  cabinets,  containing  the  4CV100,000 
(class  C)  final  and  3-1000Z  (class  B)  intermediate  amplifiers.  The 
Platteville  intermediate  amplifier  is  driven  by  a  linear  ENI  240L  driver,  a 
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linear  Mini  Circuit  ZFL  500  preamplifer,  a  variable  attenuator,  a  manual 
phase  shifter,  an  electronic  phase  shifter  (Merrimac  PSE-4,  0-360  degree 
voltage  variable  phase  shifter),  and  finally  the  output  port  of  a  10  way  RF 
power  splitter,  the  input  of  which  is  connected  to  a  frequency  synthesizer 
which  drives  the  system  at  low  level-  The  other  seven  antennas  are 
identically  connected  to  the  other  output  ports  of  the  10  way  power  splitter. 

The  array  can  be  amplitude  modulated  at  audio  frequencies  through  a  TO. 
circuit,  driven  by  another  audio  frequency  synthesizer.  The  lowest  AM 
frequency  is  presently  1S6  Hz,  due  to  one  Henry  chokes  in  each  4CV100.000 
plate  power  supply.  For  even  lower  frequency  modulations,  half  of  the  antennas 
are  shifted  in  phase  by  180°  at  the  modulation  frequency.  This  change  in  the 
phasing  of  the  array  destroys  the  central  lobe  of  the  radio  beam.  Such  a 
method  of  antenna  phase  modulation  has  been  used  to  modulate  the  array  at 
very  low  frequencies  (0.1-160  Hz)  without  stressing  the  final  plate  power 
supplies  or  harming  the  diesel  generators  through  ar.y  sympathetic  mechanical 
resonances. 

The  most  recent  improvement  in  the  whole  RF  system  has  been  the 
introduction  of  a  Personal  Computer  (PC)  to  contol  and  manage  the  operation 
of  the  eight  transmitters.  The  PC  now  controls  the  second  frequency 
synthesizer  which  amplitude  modulates  the  transmitters  at  audio  frequencies. 

It  also  controls  the  phasing  of  the  antennas  for  the  purposes  of  modulating 
the  antenna's  far  field  pattern  at  very  low  (ULF)  frequencies  and  for 
pointing  the  RF  beam  as  much  as  30°  off  the  vertical.  Its  description  will 
be  the  subject  of  a  subsequent  paper. 

A  Platteville  amplifier  is  shown  schematically  in  Figure  5.  The  final 
Eimac  4CV100,000  power  tetrode  has  its  screen  grid  grounded,  and  as  a  result 
requires  two  power  supplies.  The  plate  supply  is  it  KV.  The  screen  supply 


is  contlnously  variable  <0-1200  V)  and  is  used  to  control  the  output  RF 
power.  It  is  connected  between  the  ground  and  the  tube's  cathode.  At  low 
cathode  voltage  the  tube  is  cut  off  and  the  RF  output  is  zero.  With 
Increasing  cathode  voltage  the  tube  conducts  periodically  in  the  class  C 
mode. 

The  Eimac  3-1000Z  drives  the  control  grid  of  the  final  amplifier  through 
its  own  resonant  plate  circuit.  It  operates  as  a  Class  B  grounded  control 
grid  triode.  Its  cathode  is  connected  to  ground  by  a  SO  ohm  resistor  that 
also  terminates  the  cable  from  the  remote  ENI  240L  amplifier.  Both 
amplifiers  have  to  be  retuned  whenever  the  operating  frequency  is  changed. 

They  are  described  in  more  detail  in  the  original  Platteville  paper  [Ibldl. 

The  ENI  and  Minicircuit  amplifiers  are  linear  types  which  require  no  tuning. 

The  Bird  power  meters  at  the  output  of  the  Platteville  amplifers  monitor 
the  forward  and  reflected  powers  in  the  coaxial  lines  to  the  individual 
antennas.  Whenever  high  standing  waves  are  present  in  a  line,  the  Bird  meter- 
circuitry  trips  off  the  final  amplifier's  high  voltage  supply.  Additional 
trips  are  included  in  the  plate  and  screen  current  meters  on  each  4CV100.000 
tube.  The  <UCLA  constructed)  4CV100.1000  power  supplies  are  only  shown 
symbolically  in  Figure  5.  In  reality,  both  are  unregulated  three  phase 
rectified  units,  with  1  H  filter  chokes  and  10  uF  output  filter  capacitors. 

The  Platteville  transmitters, along  with  their  individual  power  supplies, 
were  mounted  in  pairs  within  surplus  highway  trailers.  This  reassembly  was 
all  done  at  UCLA.  On  completion,  the  trailers  were  shipped  to  the  site  and 
attached  to  the  side  of  a  specially  constructed  building  on  the  edge  of  the 
antenna  field,  as  seen  in  Figure  1.  The  central  building  now  functions  as 
the  transmitter  operation  center.  Figure  6  is  a  photograph  of  the  resulting 
transmitter  complex.  The  southwest  corner  of  the  generator  building  is  seen 
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In  the  lower  left  edge  of  the  picture.  The  transformer  yard  is  between  the 
generator  building  and  transmitter  complex.  The  antenna  field  is  in  the 
background.  The  picture  was  taken  on  the  same  day  as  Figure  2. 

B.  TRANSMITTER  COOLING  SYSTEM 

The  4CV100,000  tubes  were  designed  to  be  cooled  by  the  evaporation  of 
water.  Each  (massive  copper)  anode  fits  "upside-down"  in  a  EIMAC  boiler, 
filled  with  distilled  water,  as  shown  schematically  in  Figure  7.  A  side  arm 
connects  the  boiler  to  an  overhead  heat  exchanger  (also  of  EIMAC  manufacture) 
via  an  electrically  insulating  4  inch  diameter  pyrex  glass  tube.  Condensed 
water  is  brought  back  to  a  make  up  reservoir  on  the  side  of  the  amplifier 
cabinet.  Water  is  then  returned  to  the  boiler  via  a  12  foot  long 
(electrically  insulating)  p)ece  of  1  inch  ID  silicone  rubber  hose.  Since  the 
boiler  and  tube  anode  are  both  at  a  combined  RF  and  DC  voltage  of  32  KV  peak, 
the  pyrex  vapor  line  and  long  rubber  return  have  to  be  carefully  maintained. 
The  boiler  and  tube  are  supported  (off  ground)  by  6  inch  long  ceramic 
insulators. 

In  addition,  the  4CV100,000's  tank  circuit  inductor  (made  from  2  inch 
diameter  copper  tubing)  and  RFC  choke  are  both  cooled  with  distilled  water 
from  the  make  up  reservoir.  These  cooling  schemes  were  all  part  of  the 
original  Platteville  design  [Ibid]. 

To  run  the  Platteville  amplifiers  in  Alaska,  four  separate  secondary 
(glycol)  cooling  loops  were  built.  Each  has  its  own  radiator  and  fan,  which 
can  be  seen  in  Figure  6,  on  the  ends  of  the  trailer  roofs.  Each  glycol  loop 
circulates  water  through  two  4CV100,000  steam  heat  exchangers,  which  are 
mounted  in  the  small  sheds  on  the  roof  of  each  trailer  (seen  also  in  Figure 
6).  In  addition,  one  of  each  pair  of  glycol  loops  is  used  to  cool  (through 
another  heat  exchanger)  the  distilled  water  system  for  four  transmitters. 


This  secondary  loop  can  be  completely  bypassed,  to  avoid  freezing  of  the 
distilled  water  during  the  winter  by  glycol  cooled  to  sub  zero  temperatures 
by  the  ambient  air. 

Air  is  also  used  to  cool  the  interior  of  the  Platteville  cabinets.  It  is 
drawn  through  a  floor  duct.  During  winter,  the  duct  is  closed  to  prevent 
freezing  of  components  by  the  cold  air. 

To  further  prevent  freezing  of  distilled  water  and  harm  to  the  electrical 
components,  the  whole  transmitter  complex  is  heated  by  forced  air  and  backup 
electrical  heating  systems.  Winter  temperatures  in  Alaska  can  even  exceed 
the  low  temperature  ratings  on  many  oil  filled  capacitors.  An  alarm  system 
has  also  been  installed  to  monitor  both  the  coolant  and  air  temperatures 
throughout  the  transmitter  complex. 

The  water  cooling  system  is  completely  interlocked.  The  transmitters 
cannot  operate  without  proper  water  flow.  Failure  of  any  loop  shuts  down  its 
4CV100,000  plate  supply.  For  a  transmitter  to  operate,  water  flows  have  to 
be  correct,  plate  and  screen  currents  below  maximum  rated  values  of  12  and  1 
amperes  respectively,  and -reflected  power  from  the  antennas  below  S  KW.  The 
complete  water  cooling  system  is  shown  schematically  in  Figure  7. 

C.  DIESEL  ELECTRIC  GENERATORS 

The  HIPAS  site  is  at  the  very  end  of  the  Golden  Valley  Electric  Company's  line 
from  Fairbanks.  At  the  site,  the  commercial  line  has  only  50  KW  capacity. 

The  problem  of  powering  the  eight  Platteville  transmitters  was  solved  by 
shipping  three  1500  horsepower  <1.2  KVA  each)  4800  volt-3  phase  diesel 
electric  generators  from  surplus  sites  in  California.  Two  of  the  generators 
are  now  used  to  operate  (4  each)  transmitters,  via  the  transformer  yard, 
which  also  powers  the  water  cooling  systems.  The  third  generator  is  a  source 
of  spare  parts. 
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The  transmitter  and  generator  control  circuity,  building  lighting,  and 
building  heating  are  all  run  off  the  Golden  Valley  line. 

After  the  generators  had  been  installed  on  concrete  foundations,  the 
generator  building  was  erected.  The  third  generator  bay  now  functions  as  a 
combined  meeting  and  recreation  hall.  The  second  floor  portion  of  the 
building  also  functions  as  the  site  office,  a  library,  a  conference  room,  and 
storage  space.  The  location  of  this  building  was  shown  in  Figure  1;  an  edge 
of  it  was  seen  in  Figure  6. 

Figure  8  is  a  photograph  of  one  of  the  diesel-electric  generators,  which 
had  recently  been  upgraded  with  pneumatic  starters,  one  of  which  can  even  be 
seen  in  the  picture. 

When  the  system  is  radiating  1  MW  RF,  the  two  diesels  together  consume 
120  gallons  of  fuel  per  hour. 

D.  PULSED  TRANSMITTER 

Prior  to  the  installation  of  the  Platteville  transmitters  at  the  HIPAS 
site,  a  1.8  MW  pulsed  sytem  was  built  and  used.  It  was  lost  (March  1984)  in 
a  building  fire.  Though  this  pulsed  system  is  no  longer  available,  it  is 
reviewed  because  it  may  be  the  key  to  achieving  very  high  (100  GW  ERP)  peak 
powers  with  a  new  array.  It  is  a  resource  since  all  the  drawings  are  still 
avialable,  the  engineering  is  completed,  and  the  system  successfully  operated 
at  the  site. 

The  original  pulsed  system  was  a  master  oscillatoi — amplifier  chain,  not 
unlike  the  present  (Figure  4)  CW  system.  It  started  with  commercial  ENI 
403LA  (37  db)  and  ENI  240L  (50  db)  linear  preamplifiers,  an  Intermediate 
Power  Amplifier  (IPA;  with  21  db  gain  and  EIMAC  4CX250BC  tetrode),  a  Power 
Amplifier  (PA,  with  17  db  gain  and  EIMAC  4CX5000A  tetrode),  and  finished  with 
a  Final  Power  Amplifier  CFPA,  with  17  db  gain  and  a  Machlett  ML  7560  triodel. 
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The  FPA  had  a  26  KV  supply  with  a  166  microfarad  <45  KJ)  capacitor  bank. 
Pulsed  RF  output  of  the  system  was  1.8  MW  for  one  millisecond.  The  whole 
amplifier  is  described  in  more  detail  in  other  publications  [Wong,  1983; 
Santoru,  19841. 

We  have  recently  considered  reproducing  12  of  these  pulsed  high  power 
amplifiers  for  the  purpose  of  achieving  an  ERP  of  100  GW  with  a  new  array. 

III.  ANTENNA  SYSTEM 

The  present  antenna  array  was  shown  earlier  in  Figures  1,  2,  and  6. 

When  all  the  antennas  are  driven  in  phase,  the  radiated  beam  is  concentrated 
vertically  overhead.  The  calculated  radiation  pattern  for  this  case  is 
included  in  Figure  9.  The  center  of  this  pattern  (see  Table  I)  has  a 
theoretical  gain  at  2.85  MHz  of  18.4  dbi  and  a  half  power  width  of  21.5 
degrees,  which  corresponds  to  a  diameter  of  30  Km  at  100  KM  range.  The  sun 
and  moon,  in  comparison,  are  only  one-half  of  a  degree  wide  <10"^  radians). 

By  adjusting  the  phases  of  the  antennas  relative  to  one  another  the  beam 
pattern  can  be  changed.  An  example  is  shown  in  Figure  9,  which  corresponds 
to  the  case  where  the  beam  is  tipped  30°  off  the  vertical.  At  this  angle, 
the  peak  intensity  is  a  factor  of  two  lower  than  the  vertical  case  (Figure 
6),  while  the  side  lobe  power  has  increased. 

The  two  dipoles  of  each  antenna  are  physically  supported  46  feet  above 
the  ground  by  a  central  (original)  steel  tower  and  by  four  wooden  telephone 
poles  topped  with  G-10  fiber  glass  standoffs,  as  seen  in  Figure  2.  The 
central  tower  also  supports  the  single  "bazooka  balun"  which  converts  the 
unbalanced  coax  cable  input  to  a  balanced  drive. 

Figure  3  showed  that  there  are  two  ways  of  connecting  the  two  antennas  on 
each  tower  together.  One  connection  produces  a  right  circularly  polarized 
radiation  pattern  (called  the  O  or  ordinary  mode).  The  other  connection 
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produces  a  left  circular  pattern  (called  the  X  or  extraordinary  node).  For 
circularly  polarized  emission,  the  currents  in  the  two  orthogonal  dipoles 
have  to  be  90°  out  of  phase  with  respect  to  one  another.  For  this  to  happen, 
both  antennas  have  to  be  specially  tuned;  first,  the  east-west  dipoles  are 
disconnected  from  the  (already  tuned)  triaxial  balun.  The  north-south 
dipoles  are  then  tuned  to  have  an  impedance, 

Zns  =  SO  +  jSO  ohms  (1) 

at  the  transmitter  operating  frequency.  Next,  the  tuned  north-south  dipole 
is  disconnected  from  the  balun  and  the  east-west  dipole  reconnected.  This 
dipole  is  then  tuned  to  an  impedance  of 

Zew  =  SO  -  jSO  ohms.  (2) 

When  the  two  antennas  are  reconnected  together  in  parallel  (on  the  balun)  the 
resulting  impedance  Zr  of  the  combination  is, 

1/Zr  =  <1/Zns)  +  <1/Zew)  =  (50+j50)-1  +  (50-J50)-1  =  1/50,  (3) 

or  SO  ohms  of  real  resistive  impedance,  which  properly  terminates  the  50  Ohm 
coaxial  cable  that  drives  the  antenna.  When  the  parallel  combination  is 
driven  by  an  RF  voltage  of  V,  the  currents  to  the  dipoles  are  the  applied  RF 
voltage  divided  by  the  individual  dipole  impedances;  namely, 

Ins  =  v/zns  =  V<1  ‘  J>/100  ohES  <4) 


and 


lew  =  V/Zew  =  V(1  +  J)/100  ohms 

As  can  be  seen,  the  two  current  are  phased  90°  with  respect  to  one  another. 
The  current  into  the  common  connection  is  the  sum  of  the  two;  namely, 

Iin  =  Jns  +  lew  =  v/s0  ohas  <S> 

When  the  antenna  currents  are  in  quadrature,  the  radiation  from  the  antenna 
is  circularly  polarized.  Right  or  left  circular  polarization  is  now  only 
determined  by  the  parallel  connection  of  the  east-west  dipole  to  the  north- 
south  dipole;  namely,  whether  the  east  half  dipole  is  connected  to  the  south 
half  dipole  or  to  the  north  half  dipole,  as  diagrammed  earlier  in  Figure  3. 
This  rather  elegant  scheme  for  producing  a  right  or  left  circularly  polarized 
radio  beam  was  suggested  to  us  by  the  Platteville  design. 

Originally,  the  HIPAS  array  was  changed  between  X  and  0  modes  by  climbing 
the  antenna  towers  and  manually  Interchanging  the  east-west  dipole 
connections.  We  recently  11989)  developed  an  antenna  switch,  which  now 
let  us  quickly  change  the  antenna  polarizations  manually  from  the  ground. 

Soon  ,  electric  motors  will  be  added  and  the  polarization  changed  from  the 
control  room. 

Each  balun  also  has  to  be  tuned  to  the  transmitter’s  operating  frequency 
before  the  antennas  can  be  tuned.  The  balun  is  basically  a  quarter  wave 
resonator,  orginally  resonant  at  4.9  MHz.  A  0-1800  pF  vacuum  variable 
capacitor  was  placed  between  the  outer  coaxial  balun  sleeve  and  the  outer 
conductor  of  the  coxial  feed  line,  so  that  the  balun  circuit  could  also  be 
resonated  at  lower  operating  frequencies.  To  tune  a  balun,  the  antennas  are 
disconnected.  Next,  the  center  balun  conductor  is  connected  to  the  outer 
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most  (triax)  conductor  by  a  low  inductance  copper  shorting  strap.  The 
variable  vacuum  capacitor  Is  then  tuned  until  the  transmission  line  appears 
to  have  an  open  circuit  at  the  transmitter's  operating  frequency. 

Afterwards,  the  short  is  removed  and  the  dipoles  individually  tuned  as 
described  above.  The  tuning  of  a  balun  and  its  antennas  is  expedited  by  the 
use  of  a  commercial  network  analyzer  such  as  the  Hewlett  Packard  87S3A. 

For  frequencies  below  4.9  HHz,  the  antennas  have  to  be  lengthened  and 
retuned.  This  is  now  done  by  adding  a  pair  of  conductors  to  the  ends  of  each 
dipole.  The  ends  of  these  conductors  are  each  terminated  with  a  20 
centimeter  diameter  stainless  steel  kitchen  bowl,  v+iich  functions  as  a  corona 
shield.  Polypropolene  rope  is  attached  to  a  "hook  eye"  within  the  bowl.  By 
pulling  the  ropes  from  the  ground,  the  added  conductors  or  "end-straps”  are 
extended  in  a  straight  line  away  from  the  ends  of  the  dipoles.  A  dipole  is 
then  retuned  by  spreading  the  end-straps  away  from  one  another  and  by 
adjusting  their  angle  with  respect  to  the  ground.  The  spread  between  the  two 
straps  primarily  changes  the  reactance  of  the  antenna,  while  the  angle  of 
approach  to  the  ground  principally  varies  the  antenna's  resistance.  The  end- 
straps  on  either  end  of  a  dipole  are  manipulated  symmetrically  until  the 
dipole  fulfills  either  equation  1  or  2  at  the  new  operating  frequency. 

Tuning  a  dipole  is  a  five  man  operation,  with  one  person  on  each  end-strap 
rope,  and  the  fifth  manning  the  network  analyzer,  connected  at  the  output  of 
the  power  amplifier  (in  the  transmitter  building).  Once  the  end-straps  are 
properly  positioned,  stakes  are  driven  In  the  ground  and  the  ropes  "tied 
off."  Each  new  operating  freqeuncy  requires  a  new  set  of  end-straps  and  stake 
locations.  The  antenna  end-straps  and  corona  bowls  are  seen  clearly  in 


Figure  2.  The  end-straps  are  actually  lengths  of  RG-9U  coaxial  cables  which 
use  the  outer  braid  as  the  RF  conductor.  By  this  technique  we  have  achieved 
VSW's  of  l.i. 

The  "dog  house"  on  the  top  of  the  antenna,  and  seen  In  Figure  2,  protects 
the  balun  and  antenna  connections  from  rain,  snow,  and  ice.  With  these 
recent  protective  houses,  the  system  can  be  operated  during  rain  storms,  as 
well  as  in  the  winter.  In  the  past,  the  open  antenna  connections  became 
encased  in  ice  and  the  system  could  not  be  operated. 

As  noted  earlier,  the  irrigation  coax  lines  and  baluns  are  pressurized  to 
about  1  psi  gauge  with  dry  air  to  prevent  the  accumulation  of  water  or  ice 
inside.  This  pressurization  scheme  is  another  reason  that  the  system  can  be 
operated  in  the  winter  (-30°).  A  commercial  refrigeration  system  is  now  used 
to  dry  the  air. 

The  present  125  KW  power  capability  of  each  antenna  is  due  to  a  three 
inch  diameter  right  angle  coax  bend  between  the  Plattevllle  coax  and  the 
vertical  balun.  At  Plattevllle,  the  individual  transmitters  could  all  be  run 
at  160  KW  outputs.  To  operate  again  safely  at  this  power,  the  baluns  and 
right  angle  bends  would  all  have  to  be  enlarged  in  diameter.  The 
transmitters  have  however  been  tested  at  160  KW  for  short  durations. 

Loop  couplers  (shown  in  Figure  A)  have  been  installed  in  each  antenna 
coax,  close  to  the  transition  to  its  balun,  for  the  purpose  of  phasing  the 
antennas.  The  individual  monitoring  loops  are  connected  to  a  1000  foot  long 
piece  of  RG-8  cable,  which  ends  in  the  transmitter  control  room.  With  equal 
length  cables  on  each  monitoring  loop,  the  system  can  easily  be  phased  with 
either  a  high  frequency  X-Y  oscilloscope  or  a  Vector  voltmeter.  Excess  cable 
is  left  on  the  individual  shipping  spools  so  that  the  remaining  cable 
functions  as  chokes  to  surpress  RF,  induced  on  the  outer  conductor,  from 
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getting  Into  the  control  room  and  effecting  the  phase  measurements.  The 
monitoring  cables  to  the  more  distant  antennas  are,  for  the  same  reason, 
wrapped  around  large  ferrite  cores  before  entering  the  control  room. 

IV.  PERFORMANCE  REVIEW 

The  HIPAS  system,  as  described  above,  operates  regularly  at  1  MW  total 
output.  It  has  produced  some  remarkable  ionospheric  effects,  particularly 
when  the  ionosphere  critical  frequencies  match  the  transmitter  frequency  and 
the  array  operates  in  the  O-mode.  This  usually  means  waiting  for  the 
ionosphere's  plasma  frequency  to  either  fall  in  the  evening  through  the 
transmitter  frequency  or  rise  during  the  morning.  Unlike  at  Platteville,  the 
antennas  are  not  broadband.  Also  the  frequency  agility  of  the  amplifiers 
(originally  possible  at  Platteville)  has  not  yet  been  implemented. 

The  heater  has  also  excited  ELF  waves  by  modulating  the  ionospheric 
currents,  with  the  array  radiating  mainly  in  the  X-mode  continuously  for  up 
to  24  hours.  These  effects  are  the  subjects  of  companion  papers. 

V.  FUTURE 

The  HIPAS  site  now  represents  one  MW  of  installed  RF  power  at  auroral 
lattitudes.  The  total  radiated  power  could  be  increased  to  1.3  MW  by 
operating  each  of  the  Platteville  amplifiers  again  at  160  KW.  Such  an 
improvement  will  require  new  baluns,  and  would  increase  the  ERP  to  88  MW. 

We  could  even  increase  the  radiated  power  to  200  KW.  the  limit  of  the 
4CV100,000  tubes.  This  improvement  would  double  the  ERP  of  the  present 
array  to  140  MW  and  would  require  higher  voltage  plate  supplies,  as  well  as 
the  new  baluns. 

We  are  now  considering  other  ways  of  increasing  the  system  ERP  by 
dramatically  Increasing  both  the  size  of  the  array  and  also  increasing  the 
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operating  frequency.  One  idea  envisions  a  new  24  x  24  or  S76  dipole  array, 
830  x  830  meters  on  a  side  (i.e.,  2700  x  2700  ft),  operating  at  7.S  MHZ  (40 
meter  wavelength).  Such  an  array  would  have  a  gain  of  37  dbi  and  a  beam  width 
of  only  2.S0.  With  twelve  (12)  new  pulsed  1.8  MW  peak  transmitters,  identical 
to  the  one  originally  built  for  HIPAS  (see  II. D),  a  peak  ERP  of  110  GW  would 
be  realized.  At  100  Km  range,  this  system  would  achieve  a  flux  of  0.9  W/m^ 
or  a  peak  RF  (vacuum)  amplitude  of  26  volts/meter.  Flux  values  for  other 
ranges  are  given  in  Table  II. 

With  twelve  of  the  present  Platteville  transmitters,  operating 
continuously  at  160  KW  each,  an  ERP  of  10  GW  would  be  realized.  These  upgrade 
studies  are  on  going  and  will  be  the  subject  of  future  papers. 

VI.  SUMMARY 

A  "new"  1  MW  (70  MW  ERP)  continous  wave  radio  frequency  ionospheric 
heater  is  now  operational  in  Alaska.  The  radiated  beam  can  be  either  O  or  X 
mode  of  circular  polarization.  The  heater  uses  eight  RF  power  amplifiers 
from  the  Platteville  heater.  It  has  been  operated  at  frequencies  of  2.85, 
3.349,  4. S3,  and  4.905  MHz.  It  presently  (April  1988-  December  1989) 
operates  at  2.85  MHz,  which  corresponds  to  the  second  harmonic  of  the 
electron  cyclotron  frequency.  The  heater  can  be  amplitude  modulated  160  Hz  - 
20  KHz,  and  effectively  amplitude  modulated  at  lower  (ULF)  frequencies  by 
dephasing  half  of  the  antennas  in  the  array.  Beam  pointing  (+  30  )  is  now 
controlled  by  a  personal  computer.  The  transmitters  are  powered  by  two  1500 
HP  diesel  electric  generators.  The  system  is  capable  of  almost  year  round 
operation  and  has  been  operated  without  any  interruption  for  as  long  as  12 
hours. 

Future  improvements  include  the  construction  of  a  larger  array  operating 
at  higher  frequencies  to  minimize  the  array  size.  Arrays  with  37  dbi  gains 
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are  being  considered,  along  with  pulsed  amplifiers  to  raise  the  ERP  to  -100  GW. 
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FIGURE  CAPTIONS 


Figure  1:  Schematic  of  the  present  HIPAS  ionospheric  radio  frequency  heater 
showing  the  generator  building,  transformer  yard,  transmitter  complex,  and 
antenna  field. 

Figure  2:  Photograph  of  one  of  the  HIPAS  antennas  taken  April  1988.  The 
antennas  were  orginally  constructed  for  4.9MHz  operation.  Their  operating 
frequencies  were  lowered  by  the  addition  of  the  enc-straps  visible  in  the 
picture.  When  the  picture  was  taken  the  antenna  was  tuned  for  2. 85  MHz. 

Also  seen  in  the  picture  is  the  balun  and  rigid  coax  to  the  antenna’s 
nominal  12S  KW  amplifier.  The  dipoles  are  48  feet  above  ground. 

Figure  3:  Schematic  view,  from  above,  of  a  HIPAS  antenna,  showing  the 
connections  to  the  balun  to  produce  either  right  ( O-mode )  or  left  (X-mode> 
circularly  polarized  radiation  patterns.  Also  included  in  this  picture  is 
the  direction  of  the  earth's  magnetic  field  (into  page)  and  direction  of 
rotation  of  a  free  electron  in  the  field. 

Figure  4:  Schematic  of  the  HIPAS  radio  frequency  system. 

Figure  5:  Schematic  of  one  of  the  Platteville  transmitters. 

Figure  6:  Photograph  (April  1988)  of  the  transmitter  building,  with  the  array 
and  transmission  lines  in  the  background.  The  edge  of  the  generator  building 
is  seen  in  the  lower  left  corner  of  the  picture.  The  transformer  yard  is 
also  in  part  visible.  The  four  radiators  for  the  glvchol  cooling  loops  are 
on  the  far  ends  of  the  building. 


Figure  7:  Schematic  of  the  cooling  system 

Figure  8:  Photograph  of  one  of  the  (two)  1500  HP  diesel  electric  generators, 
used  to  drive  the  RF  system. 

Figure  9:  Top:  Calculated  HIPAS  beam  pattern  at  2.85MHz  when  all  eight 
antennas  are  in  phase;  namely,  18. S  dbi  gain  over  isotropic,  with  a 
divergence  21.5  degrees  (FWHM);  bottom:  Calculated  HIPAS  beam  pattern  for 
2.85  MHz,  when  the  antennas  are  phased  to  tip  the  RF  beam  30°  off  the 
vertical . 
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TABLE  I 


HIP AS  OPERATING  PARAMETERS 


Frequencies 

Calculated 

Antenna 

ERP 

<@1.2MW> 

FLUX 

(@100km) 

Peak  Amplitude 
(Vacuum  §  100km) 

MHZ 

Gain  (dbl) 

(MW) 

rnW/rn^ 

Volts/m 

4.905 

17.7 

71 

0.56 

0.64 

4.530 

3.349 

2.85* 

18.4 

84 

0.67 

0.71 

•1987,  1988,  and  1989  operating  frequency 
Number  of  crossed  half  wave  dipole  antennas:  8 
Array  diameter:  207  meters  <680  feet) 

Polarization:  Left  (X  mode)  or  Right  (0  mode)  Circular 
AM  Modulation:  160  Hz  -  20  KHz 
Antenna  Phase  Modulation:  >  0.1  Hz  -  156  Hz 
Beam  Direction:  <30°  off  Zenith,  at  >0.01  seconds 
Geographic  Location:  62°  52'  22"  N,  146°  50 ’  06"  W. 

Magnetic  Location:  28°  02'  E  declination,  76°  30'  dip  angle 
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TABLE  II 


PERFORMANCE 

OF  THE  PROPOSED 

110  GW  ERP  ARRAY 

RANGE 

FLUX 

RF  AMPLITUDE 

(Kilometers ) 

<W/m2)  i 

:EpeakRF»volts/m> 

25 

14.4 

104 

50 

3.6 

52 

100 

0.9 

26 
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APPENDIX  II 


Paper  prepared  by  R.G. Brandt,  ONR  Program  Manager,  for  presentation  at 
the  May  1990  AGARD  meeting  in  Tromso,  Norway.  This  paper  summarises  both 
the  URI  results  as  well  as  later  experiments  on  ULF  generation  with  the 
HIPAS  facility; 


R.G. Brandt,  Ionospheric  Modification  Research  at  HIPAS, 


IONOSPHERIC  MODIFICATION  RESEARCH  AT  HIPAS 

by 

Richard  G.  Brandt 
Office  of  Naval  Research 
800  N.  Quincy  Street 
Arlington,  VA  22217 
U.S.A. 


SUMMARY 

The  HIPAS  ionospheric  heating  facility  radiates  a  total  power 
of  1.2  MW  with  an  ERP  of  84  MW.  It  presently  operates  at  an  HF 
frequency  of  2.85  MHz  but  is  tunable  to  about  5  MHz.  Electro jet 
modulation  experiments  have  been  conducted  at  frequencies  from  5  Hz 
to  5  kHz.  The  magnetic  field  amplitudes,  measured  close  to  the 
heater,  can  be  1  pT  or  larger  under  very  strong  electro jet 
conditions.  Even  under  much  weaker  conditions  when  the  amplitudes 
are  highly  variable,  the  phase  of  the  ELF  signal  is  relatively 
stable.  The  efficiency  of  converting  HF  to  ELF  is  presently  too 
low  for  a  practical  communication  system.  Beam  painting  has  been 
proposed  as  a  method  for  improving  the  conversion  efficiency  in  D 
region  heating  by  causing  a  much  larger  area  of  the  ionosphere  to 
radiate  coherently;  this  concept  will  be  tested  using  microsecond 
beam  steering.  Use  of  shorter  heating  pulses  (lower  duty  cycle) 
already  seems  promising.  Even  larger  gains  are  expected  for  E 
region  heating  as  compared  to  D  region  heating. 

PREFACE 

This  review  of  ionospheric  modification  research  at  HIPAS 
starts  with  a  brief  description  of  the  facility  which  began 
operating  at  the  megawatt  average  power  level  in  1986.  The 
research  conducted  at  HIPAS  since  that  time  is  described.  This 
research  has  focused  on  the  generation  of  ELF  radiation.  The 
review  concludes  with  an  indication  of  future  research  directions. 

1.  THE  HIPAS  FACILITY 

The  HIPAS  ionospheric  heater  is  located  in  the  auroral  zone  40 
km  east  of  Fairbanks,  Alaska.  The  acronym  HIPAS  stands  for  High 
Power  Auroral  Stimulation.  It  was  originally  conceived  as  a  high 
pulse  power  heater,  rather  than  high  average  power,  because  of  the 
recognition  that  many  plasma  nonlinearities  are  excited  by  large 
instantaneous  electric  fields  (Wong  and  Santoru,  1981].  When  HIPAS 
became  operational  in  October  1981,  it  delivered  1  ms  HF  pulses 
with  a  peak  power  of  2  MW  at  a  pulse  repetition  rate  of  1  Hz.  The 
average  power,  therefore,  was  only  2  kW,  and  the  antenna  gain  was 
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17.7  dB  at  a  radiating  frequency  of  4.905  MHz.  This  facility  was 
initially  used  for  probing  short-term  nonlinear  plasma  phenomena, 
especially  enhancement  of  the  incident  electric  fields  at  the 
resonant  altitude  and  creation  of  electromagnetic  cavitons  via 
ponderomotive  forces.  These  cavitons  are  localized  depletions  in 
electron  density  with  scale  sizes  from  1  to  100  m. 

In  1982  several  events  occurred  which  altered  the  original 
research  plan.  First,  in  March  1982  the  incoherent  scatter  radar 
at  Chatanika  was  moved  to  Greenland.  Much  of  the  detailed  work  on 
caviton  structure  and  dynamics  was  subsequently  carried  out  at  the 
Arecibo  facility  with  its  excellent  incoherent  scatter  radar. 
Second,  in  late  1982  the  HF  transmitters  from  the  Platteville 
heater  became  available,  and  a  decision  was  made  to  incorporate 
these  transmitters  in  the  HIPAS  facility,  transforming  it  into  a 
high  average  power  heating  facility. 

Eight  of  the  Platteville  transmitters  are  now  operating  at 
HIPAS,  each  connected  to  an  individual  element  in  the  transmitting 
antenna  array  [Wong  et  al,  1990].  This  modular  design  allows  the 
phase  of  the  high  power  radiation  emitted  by  each  element  in  the 
array  to  be  controlled  at  the  low-level  input  to  each  transmitter 
and  thus  permits  rapid  beam  steering  and  related  modulation 
techniques.  Each  transmitter  is  normally  operated  continuously  at 
100  to  150  kW,  although  the  maximum  rating  is  200  kW.  Instead  of 
relying  on  commercial  prime  power,  two  megawatt  diesel  electric 
generators  were  installed  at  HIPAS  so  that  operational  costs  would 
be  minimized. 

The  antenna  consists  of  a  circular  (208  m  diameter)  array  of 
seven  crossed  half-wave  dipoles  approximately  one  quarter  of  a 
wavelength  above  the  ground  plus  one  similar  crossed  dipole  in  the 
center.  Although  the  antenna  elements  were  designed  to  be  resonant 
at  4.905  MHz,  a  novel  scheme  was  developed  [2]  for  retuning  the 
antenna  to  lower  frequencies  by  making  use  of  a  pair  of  wires 
attached  to  the  end  of  each  dipole  element;  the  angle  between  these 
two  wires  and  the  angle  these  wires  make  with  the  ground  determine 
the  new  resonant  frequency.  Using  this  technique,  the  array  has 
been  successfully  retuned  and  operated  at  the  additional 
frequencies  of  4.503,  3.349  and  2.805  MHz.  The  antenna  gain 
remains  essentially  constant  at  about  18  dB.  At  2.85  MHz  the 
calculated  effective  radiated  power  (ERP)  is  84  MW  at  a  total 
radiated  power  of  1.2  MW. 

With  the  upgrade  of  HIPAS  to  the  megawatt  average  power  level, 
the  research  directions  changed.  These  new  research  activities 
capitalized  on  the  existing  diagnostic  capabilities  which  included 
ionosondes,  an  HF  imaging  radar,  magnetometer  chain  data,  riometer 
data,  and  ELF  receivers.  More  recently,  a  50  MHz  coherent  radar 
has  become  operational.  In  addition,  satellite  receivers  and 
beacons  have  been  utilized  at  various  times.  The  initial 
scientific  studies  were  concerned  with  optimal  approaches  for 
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coupling  energy  into  the  ionospheric  plasma  [Wong  et  al,  1989]. 
Beginning  in  1987  the  research  has  focused  on  the  subject  of 
generation  of  ELF  radiation  in  the  ionosphere  by  modulation  of  the 
powerful  electro jet  current  flowing  in  the  auroral  region. 

2 .  ELF/VLF  GENERATION 

The  ELF  generation  experiments  at  HIPAS  began  in  1987.  This 
team  effort  involved  A.  Ferraro  (Pennsylvania  State  University) ,  A. 
Wong  (University  of  California,  Los  Angeles) ,  D.  Papadopoulos 
(Science  Applications  International  Corporation),  J.  Olson 
(University  of  Alaska)  and  coworkers.  In  these  experiments  the  ELF 
receivers  were  located  about  50  km  west  of  HIPAS.  The  receivers 
[Baker  et  al,  1990]  used  standard  coherent  detection  technigues, 
and  the  heater  and  the  receivers  were  synchronized  using  separate 
freguency  standards.  Most  of  the  measurements  were  made  at 
modulation  freguencies  between  1  and  5  kHz  where  the  detection 
sensitivity  was  best.  The  heater  was  almost  always  amplitude 
modulated  with  100  percent  modulation  depth. 

These  experiments  continued  through  1989  during  which  time 
progressively  more  sophisticated  beam  steering  capability  was 
implemented.  Initially,  the  heater  beam  could  be  steered  by 
manually  adjusting  the  eight  phaseshifters  in  the  lines  leading  to 
the  individual  antenna  elements.  In  1988  electronic  phaseshifters 
were  introduced  and  placed  under  computer  control,  and  a 
beamforming  time  of  10  ms  was  achieved.  In  1989  this  capability 
was  upgraded  to  provide  15  jjs  beamforming  time. 

The  initial  experiments  [Ferraro  et  al,  1989]  served  to 
demonstrate  the  ability  to  generate  ELF  and  confirm  earlier  results 
obtained  at  Tromso.  As  an  illustration.  Figure  1  shows  data 
obtained  at  a  modulation  rate  of  2.5  kHz  during  a  strong  electro jet 
event  which  occurred  overhead  of  the  HIPAS  facility  as  indicated  by 
magnetometer  data.  We  see  that  the  amplitude  and  phase  of  the  ELF 
signal  remain  constant  during  the  2-minute  period  of  modulation, 
after  which  the  amplitude  falls  to  zero  and  the  phase  becomes 
random.  The  magnetic  signal  strength  is  0.6  pT.  Even  when  the 
electrojet  is  very  weak  and  the  amplitude  of  the  ELF  signal  is 
small  and  highly  variable,  the  phase  remains  remarkably  constant  as 
illustrated  in  Figure  2. 

This  phase  stability  of  the  ELF  signal  suggests  that 
phaseshift  keying  might  be  a  viable  signal  modulation  technigue  for 
an  ELF  communication  system  based  on  ionospheric  generation.  To 
explore  this  hypothesis,  several  phase  shifting  experiments  were 
performed.  The  tests  performed  in  1987  were  relatively  simple, 
consisting  of  periodically  changing  the  phase  of  the  ELF  signal  by 
180  degrees.  Some  of  these  data  [5]  are  illustrated  in  Figure  3 
and  indicate  a  very  stable  phase  for  periods  of  0.5  minute. 
Following  these  initial  biphase  stability  tests,  various  guadphase 
tests  were  conducted  in  1989  in  which  the  phase  was  periodically 
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switched  from  180  to  90  to  0  to  -90  to  -180  degrees. 

Representative  data  [5]  are  shown  in  Figure  4  and  again  indicate 
phase  stability  for  periods  of  0.5  minute.  All  these  measurements 
were  made  at  the  receiving  site  50  km  from  HIPAS  and  at  a  frequency 
of  2.5  kHz.  In  1989  several  attempts  were  made  to  receive  data 
using  a  mobile  receiving  station  located  at  more  remote  sites  in 
Alaska.  Data  obtained  from  Cantwell,  Alaska,  about  210  km  from 
HIPAS  [5]  are  shown  in  Figure  5.  Switching  between  the  various 
phase  states  is  clearly  evident,  although  the  transitions  between 
these  states  are  less  sharp. 

Other  experiments  were  performed  to  determine  whether  these 
ELF  signals  could  be  generated  at  less  than  full  power  levels. 

Some  of  these  data  [5]  are  shown  in  Figure  6.  The  curious  spikes 
in  the  data  were  caused  by  a  local  thunderstorm.  In  this 
experiment  the  modulation  frequency  was  2.5  kHz,  and  the  total 
radiated  power  was  varied  from  a  maximum  of  800  kW  to  a  minimum  of 
80  kW.  At  maximum  power  the  measured  magnetic  field  amplitude  was 
less  than  0.1  pT  which  indicates  that  there  was  no  strong  overhead 
electro jet.  Nonetheless,  ELF  signals  were  generated  at  all  HF 
transmitter  power  levels  including  at  80  kW  when  the  measured 
magnetic  field  amplitude  was  less  than  0.02  pT.  There  is  some 
indication  that  the  effect  is  beginning  to  saturate  at  the  higher 
power  levels,  suggesting  that  there  would  be  little  advantage  in 
using  a  heater  with  an  even  higher  ERP. 

3.  BEAM  PAINTING/LOW  DUTY  CYCLE  HEATING 

The  experimental  results  described  above  clearly  demonstrate 
the  ability  to  generate  ELF  signals  under  a  variety  of  ionospheric 
conditions.  Although  the  amplitude  of  these  signals  may  be 
variable,  the  phase  is  reasonably  stable.  The  efficiency  of  the 
generation  process,  however,  is  quite  low.  The  total  radiated  HF 
power  is  approximately  1  MW.  Except  under  strong  electrojet 
conditions,  the  ELF  power  is  very  modest.  When  there  is  no  strong 
electro jet  directly  overhead,  the  ELF  power  is  no  more  than  10  mW, 
corresponding  to  a  conversion  efficiency  of  10'8.  This  ELF  power 
is  calculated  from  the  measured  near-field  amplitudes  using  a  far- 
field  formula  for  dipole  radiation  in  free  space. 

Papadopoulos  and  coworkers  [Papadopoulos  et  el,  1990]  have 
developed  a  theory  which  suggests  an  approach  for  achieving 
dramatic  improvements  in  HF-to-ELF  conversion  efficiency,  namely, 
beam  painting.  For  low  altitude  heating  they  calculate  that  the 
ERP  of  the  current  HIPAS  heater  is  sufficient  to  achieve  the  full 
saturated  value  of  Hall  conductivity  modulation.  Any  further 
increase  in  power  density  would  be  wasted.  Papadopoulos  et  al  [6] 
conclude  that  the  only  way  to  achieve  an  increase  in  conversion 
efficiency  in  low  altitude  heating  is  to  cause  a  much  larger  area 
of  the  ionosphere  to  radiate  coherently.  The  predicted  ELF  power 
is  then  proportional  to  the  square  of  the  radiating  area.  In 
principle,  this  can  be  accomplished  by  beam  painting,  that  is, 
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using  the  present  beam  spot  size  but  rapidly  steering  the  beam  such 
that  many  different  spots  are  heated  and  not  allowed  to  cool  off 
before  they  radiate.  At  75  km  altitude  the  time  required  to 
increase  the  electron  temperature  by  a  factor  of  2  (a  sufficient 
increase  to  produce  most  of  the  conductivity  change)  is  typically 
about  10  ps;  this  heating  time  decreases  somewhat  at  higher  power 
densities.  At  the  same  altitude  the  cooling  time  is  of  the  order 
of  100  ps,  at  least  for  temperatures  below  2000  K.  The  concept  of 
beam  painting,  therefore,  requires  that  many  spots  be  visited  and 
that  the  temperature  rise  produced  during  each  visit  is  not 
dissipated  before  that  same  spot  is  revisited.  The  total  time 
available  for  heating  is  a  half  cycle  of  the  ELF  modulation  or 
about  5  ms.  During  that  time  the  temperature  of  each  spot  must  be 
increased,  in  cummulative  fashion,  by  at  least  a  factor  of  2. 

Papadopoulos  and  coworkerc  [6]  have  also  analyzed  the  case  of 
heating  at  100  km,  the  altitude  at  which  the  electro jet  flows. 

Here  the  situation  is  more  complex.  At  these  altitudes  the 
Pederson  conductivity  modulation  dominates  over  the  Hall 
conductivity  modulation  and  does  not  saturate  until  power  densities 
of  at  least  10  mW/m2  are  reached,  a  factor  of  10  higher  than  that 
available  at  HIPAS.  When  saturation  occurs,  the  predicted  increase 
in  conductivity  modulation  is  by  a  factor  of  100  which  leads  to  an 
increase  in  the  efficiency  of  ELF  generation  bv  a  factor  of  104. 
This  increase  would  be  obtained  without  beam  painting.  The 
practical  difficulty  in  realizing  the  high  altitude  scheme, 
especially  under  day  time  conditions,  is  the  fact  that  self 
absorption  at  lower  heights  could  deplete  the  power  density  before 
reaching  the  higher  altitudes.  Another  potential  disadvantage  of 
high  altitude  heating  is  that  the  coupling  of  the  ELF  fields  into 
the  earth-ionosphere  waveguide  may  be  more  inefficient. 

Nonetheless,  the  potential  gain  in  generation  efficiency  by  a 
factor  of  104  is  a  significant  prize  worth  pursuing,  and  several 
schemes  for  defeating  the  self  absorption  problem  have  been 
suggested.  However,  none  of  these  ideas  for  increasing  the 
efficiency  of  high  altitude  generation  of  ELF  has  yet  been  tested 
experimentally. 

The  recent  experiments  at  HIPAS  were  designed  to  maximize 
heating  in  the  low  altitude  regime  because  of  the  choice  of 
operating  frequency  (2.85  MHz  which  is  near  the  second  harmonic  of 
the  electron  cyclotron  frequency)  and  the  use  of  X  mode 
polarization  (same  sense  of  rotation  as  the  electron  cyclotron 
motion) .  An  initial  effort  was  made  to  demonstrate  that  beam 
painting  could  produce  an  efficiency  enhancement  in  generation  of 
ELF  at  low  altitudes.  At  present  the  beam  forming  time  at  HIPAS  is 
15  jis.  Each  spot  must  be  heated  almost  10  ps  during  each  visit. 
Therefore,  in  order  to  heat  10  spots,  for  example,  about  250  jjs  is 
required  which  is  longer  than  the  cooling  time  of  each  spot.  This 
simple  arithmetic  clearly  illustrates  that  a  significantly  faster 
beam  forming  time  is  required  before  an  efficiency  enhancement  can 
be  demonstrated.  The  initial  tests  at  HIPAS,  therefore,  did  not 
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yield  an  efficiency  enhancement,  but  the  dependence  of  generated 
signal  strength  on  beam  painting  parameters  was  in  accordance  with 
theory. 

Other  experiments  were  performed  [Papadopoulos  and  Ferraro, 
1989]  which  directly  relate  to  the  beam  painting  concept.  A  series 
of  measurements  were  made  in  which  the  duty  cycle  of  the  ELF 
modulation  of  the  HF  waves  was  varied.  The  ELF  modulation 
frequency  was  833  Hz,  and  the  duty  cycles  selected  were  50  percent 
(600  JL1S  on,  600  jas  off),  37.5  percent  (450  jis  on,  750  jas  off),  25 
percent  (300  jis  on,  900  ^is  off),  12.5  percent  (150  }is  on,  1050  ps 
off),  6.25  percent  (75  ^as  on,  1125  jis  off),  2.5  percent  (30  jas  on, 
1170  >is  off),  and  1.25  percent  (15  jjs  on,  1185  ps  off).  The 
amplitudes  of  the  generated  ELF  signals  were  measured  for  all  duty 
cycles,  and  the  data  are  shown  in  Figure  7.  The  error  bars 
represent  the  scatter  in  the  results  for  repeated  measurements. 

The  results  clearly  show  that  the  same  ELF  signal  strengths  are 
generated  at  25  percent  duty  cycle  as  at  50  percent  duty  cycle. 
Below  25  percent  duty  cycle  the  amplitudes  fall  off  approximately 
linearly.  The  conclusion  is  that  the  efficiency  of  generating 
magnetic  field  signals  can  be  increased  by  a  factor  of  2  and  thus 
the  efficiency  of  generating  ELF  power  can  be  increased  by  a  factor 
of  4  merely  by  decreasing  the  duty  cycle  from  50  percent  to  25 
percent,  that  is,  by  reducing  the  average  HF  power  by  a  factor  of 
2.  For  single  spot  heating  only  about  the  first  300  jis  of  the 
heating  pulse  is  effective;  the  rest  is  wasted.  Pulses  shorter 
than  300  ps  do  not  saturate  the  conductivity  modulation  although 
they  do  produce  a  measurable  effect.  The  time  between  heating 
pulses  is  longer  chan  the  cooling  time,  and  thus  energy  cannot  be 
accumulated. 

4.  ULF/ELF  GENERATION 

ELF  generation  was  also  investigated  in  the  frequency  range 
below  100  Hz  [McCarrick  et  al,  1990],  In  these  experiments  it  was 
initially  not  possible  to  use  full  amplitude  modulation  as  in  the 
higher  frequency  experiments  because  of  various  transmitter  system 
resonances.  Instead,  an  alternate  procedure  was  developed  known  as 
beam-dephasing  modulation.  In  this  procedure  all  eight 
transmitters  were  operated  continuously,  while  the  phase  of  four  of 
these  transmitters  was  switched  in  a  square-wave  fashion  between  0 
and  180  degrees  with  respect  to  the  phase  of  the  other 
transmitters.  In  this  way  the  beam  was  alternately  focused  and 
defocused  at  the  desired  modulation  rate.  When  the  beam  is 
defocused,  there  still  is  some  residual  heating  of  the  ionosphere, 
both  overhead  and  in  the  sidelobes,  and  thus  the  effective 
modulation  depth  is  3 ess  than  in  100-percent  amplitude  modulation. 
It  was  calculated  that  the  power  density  in  the  main  lobe  is 
reduced  by  more  than  75  percent  during  the  out-of-phase  half¬ 
period. 

The  receiver  used  in  these  experiments  is  located  c. t  the  NOAA 
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Gilmore  Creek  Facility,  35  km  from  HIPAS.  The  system  was  developed 
by  D.  Sentman  for  the  purpose  of  studying  the  natural  Schumann 
resonances  of  the  earth-ionosphere  waveguide.  The  receiver 
consists  of  a  pair  of  orthogonal  magnetic  coils,  a  vertical 
electric  field  sensor,  associated  amplifiers  and  signal 
conditioning  electronics,  and  a  dedicated  data  acquisition  and 
storage  system.  To  reduce  the  sensitivity  to  60  Hz  background 
radiation,  a  low-pass  filter  is  placed  between  the  coils  and  the 
amplifiers.  The  attenuation  is  4.4  dB  at  76  Hz  and  becomes  very 
large  above  100  Hz.  The  data  obtained  at  this  site  are  transmitted 
in  real  time  to  the  HIPAS  control  room  where  they  are  displayed 
together  with  environmental  data  from  the  ionosondes,  the 
magnetometer  chain,  and,  more  recently,  from  the  coherent  radars. 

Using  the  dephasing  modulation  technique  and  the  receiving 
system  just  described,  experiments  were  performed  at  the 
frequencies  of  5,  11,  21,  23,  41,  and  76  Hz.  In  Figure  8  is  shown 
a  typical  background  noise  spectrum  obtained  from  a  4 -minute 
incoherent  average  of  the  east-west  coil  output  [8].  The  first 
three  Schumann  resonance  peaks  at  7.5,  14  and  21  Hz  are  well 
defined.  In  order  to  observe  the  generated  ELF  signals  which  are 
often  weaker  than  the  natural  atmospheric  noise  background, 
coherent  detection  techniques  were  used,  including  conventional 
lock-in  amplifiers  and  digital  signal  processing  routines.  A 
typical  coherent  ELF  spectrum  from  the  east-west  coil  [8]  is  shown 
in  Figure  9.  Here  we  see  the  11  Hz  signal  generated  by  HIPAS  as 
well  as  the  third  harmonic  of  this  signal  at  33  Hz.  A  4-minute 
averaging  period  was  used.  Similar  quality  data  were  obtained  from 
the  north-south  magnetic  sensor  as  well  as  the  vertical  electric 
field  sensor.  Note  that  the  measured  magnetic  field  strength  at  11 
Hz  is  in  excess  of  1  pT. 

These  experiments  were  performed  under  a  variety  of 
ionospheric  conditions  in  order  to  correlate  the  strength  of  the 
observed  signals  with  measures  of  geomagnetic  activity.  Basically, 
whenever  the  magnetic  activity  is  high  enough,  as  determined  by  the 
Kp  index  or  by  magnetometer  chain  data,  large  ELF  signals  can  be 
generated.  In  terms  of  Kp  index,  which  is  an  indication  of  how  far 
south  in  geomagnetic  latitude  the  auroral  oval  has  expanded,  it  was 
found  that  signals  in  excess  of  1  pT  could  be  generated  whenever 
the  Kp  index  exceeded  3.  Similarly,  ELF  signal  strength  correlates 
strongly  with  strength  of  the  electro jet  as  determined  from 
magnetometer  chain  data.  Figure  10  shows  data  for  the  north-south 
magnetometer  at  College,  Alaska  during  a  time  period  corresponding 
to  an  ELF  generation  experiment  [8].  Data  from  a  single  ground- 
based  magnetometer  are  somewhat  difficult  to  interpret  since  the 
ionospheric  currents  move  in  geographic  location  as  well  as 
increase  and  decrease  in  amplitude  throughout  the  day.  However,  in 
Figure  11  we  see  the  measured  ELF  amplitudes  during  the  same  time 
period.  Both  signals  begin  at  high  levels  and  then  fall  to  minimum 
values  at  0700  UT,  followed  by  a  significant  rise.  Finally,  as  the 
magnetometer  output  decreases  and  then  changes  sign  corresponding 
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to  a  reversal  of  electrojet  current  direction,  there  is  a  decrease 
and  then  a  very  rapid  increase  in  ELF  amplitude.  At  the  same  time 
as  this  current  reversal  is  occurring,  the  otherwise  steady  phase 
of  the  ELF  signal  experiences  a  discontinuous  phase  shift. 

Other  favorable  situations  for  ELF  generation  include  the 
occurrence  of  visible  aurora  and  the  existence  of  a  sporadic  E 
layer  as  determined  by  ionosondes.  An  unfavorable  time  for  ELF 
generation  is  under  high  absorption  conditions,  as  indicated  by 
riometer  readings.  Such  high  absorption  conditions  are  caused  by 
high-energy  proton  and  electron  precipitation  during  major  solar 
events.  Under  these  conditions  most  of  the  HF  power  appears  to  be 
absorbed  at  a  very  low  altitude  and  cannot  effectively  modulate  the 
conductivity  where  ionospheric  currents  are  flowing. 

These  measurements  prove  that  the  beam-dephasing  modulation 
technique  works  successfully,  but  with  an  effective  modulation 
depth  of  less  than  100  percent.  Very  recently  [Wong  et  al,  1990] 
it  has  been  demonstrated  that  the  transmitters  can  be  amplitude 
modulated  at  frequencies  down  to  50  Hz  by  increasing  the  filter 
capacitance  in  the  high  voltage  power  supplies.  At  the  present 
time  four  of  the  eight  transmitters  have  been  so  modified.  Initial 
tests  have  been  performed  to  compare  the  performance  of  amplitude 
modulation  with  dephasing  modulation  at  76  Hz.  The  resultant 
electric  field  signal  amplitudes  are  shown  in  Figure  12.  In  all 
cases  the  total  radiated  power  was  kept  constant  at  400  kW.  The 
data  labeled  8 DM  indicates  that  all  eight  transmitters  were 
operated  at  the  50  kW  level  and  modulated  by  dephasing.  The  label 
4 DM  denotes  that  four  transmitters  were  modulated  by  dephasing, 
each  operated  at  the  100  kW  level.  When  four  transmitters  were 
full  amplitude  modulated  and  operated  at  the  100  kW  level,  the 
resultant  signal  amplitudes  are  labeled  4AM.  Comparing  the  4AM 
results  with  the  4 DM  results,  we  see  that  the  4AM  amplitudes  are 
larger  than  the  4 DM  amplitudes  by  a  factor  of  more  than  3, 
corresponding  to  an  increase  in  ELF  power  by  at  least  a  factor  of 
10.  This  increase  is  achieved  at  the  same  total  radiated  power  and 
the  same  ERP.  The  only  change  is  an  improvement  in  effective 
modulation  depth.  In  the  8 DM  case  the  ERP  is  much  higher  and  even 
larger  amplitudes  are  obtained,  indicating  that  the  conductivity 
modulation  effect  is  not  saturated  at  the  lower  power  density.  It 
is  anticipated  that  similar  gains  can  be  realized  when  all  eight 
transmitters  are  converted  to  full  amplitude  modulation. 

5.  FUTURE  RESEARCH  DIRECTIONS 

Additional  ELF  generation  experiments  will  be  performed  at 
HIPAS  later  in  1990  in  order  to  investigate  more  fully  various 
approaches  for  enhancing  the  HF-to-ELF  conversion  efficiency.  The 
beam  steering  capability  will  be  upgraded  to  a  2  )is  beam  forming 
time,  allowing  more  definitive  tests  of  the  beam  painting  concept. 
Other  experiments  will  evaluate  the  use  of  short  pulses  (low  duty 
cycle  heating)  to  improve  efficiency.  All  transmitters  will  be 


8 


modified  to  allow  full  amplitude  modulation  down  to  50  Hz.  The 
receivers  will  be  optimized  for  detection  in  the  75  to  150  Hz 
range,  and  measurements  will  be  made  at  distances  of  hundreds  of 
kilometers  from  the  heater  using  a  mobile  receiving  station.  It  is 
anticipated  that  a  single  HF  heating  frequency  will  be  used,  chosen 
to  maximize  absorption  in  the  D  region,  although  dual  frequency 
operation  would  permit  exploration  of  some  of  the  concepts  for 
improved  efficiency  via  E  region  heating.  The  measured  ELF 
amplitudes  will  be  correlated  with  real-time  environmental  data 
from  ionosondes,  the  riometer  and  the  magnetometer  chain  plus  data 
from  50 -MHz  coherent  radars  operating  near  Anchorage  and  at  HIPAS. 
The  understanding  gained  from  these  new  measurements  and  their 
interpretation  will  be  valuable  in  future  research  planning  and  in 
the  design  of  new  facilities. 
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Figure  1:  Magnitude  and  phase  of  strong  2.5  kHz  signal 
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Figure  2:  Magnitude  and  phase  of  weak  2.5  kHz  signal. 
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Figure  3:  Detected  phase  of 
2.5  kHz  signal  in  phase  keying 
test.  Phase  inverted  by  180 
degrees  at  30  second  intervals. 
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Figure  4:  Detected  phase  of  2.5  kHz  signal 
in  phase  shift  keying  test.  The  phase  is 
shifted  from  180  to  90  to  0  to  -90  to  -180 
degrees. 
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Figure  5:  Detected  phase  of  2  kHz  signal  in 
phase  shift  keying  test.  The  phase  is  shifted 
from  180  to  90  to  0  to  -90  to  -180  degrees. 

The  receiver  is  located  210  km  from  the  heater. 
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Figure  6:  Magnitude  of  received  2.5  kHz 
signals  in  power  reduction  test.  Total 
HF  power  is  reduced  from  800  to  640  to 
480  to  320  to  160  to  80  kW. 
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Figure  7:  Relative  signal  amplitudes  at  various 
heater  duty  cycles 


Figure  8:  Background  noise  spectrum 
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Figure  9:  Signal  spectrum  resulting  from  modulation 
at  11  Hz. 
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Figure  12:  Signal  amplitude  at  76  Hz  for  various  modulation 
schemes.  AM  denotes  amplitude  modulation;  DM  denotes  dephasing 
modulation;  4  and  8  indicate  number  of  transmitters  being 
modulated.  The  total  radiated  power  in  all  cases  is  400  kW. 
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ABSTRACT 

The  scaling  laws  that  control  the  efficiency  of  convening  ground  based  HF  power 
to  ELF  power  by  using  modulation  of  the  polar  electrojet  current  is  discussed.  The 
analysis  is  based  on  kinetic  calculations  of  the  modification  of  the  ionospheric  conductivity 
by  HF  waves  in  conjunction  with  the  experimental  results  reponed  from  the  Tromso  Max 
Planck  and  the  Alaska  HIP  AS  facilities.  It  is  shown  that  the  efficiency  can  be  increased  by 
more  than  a  factor  of  104  by 

(i)  Using  phasing  to  sweep  the  antenna  beam  over  an  area  spanned  by  a 
maximum  tilt  of  35°.  on  a  timescale  faster  than  the  cooling  rate  at  the  heating 
altitude  ( <  100  (isec). 

(ii)  Localization  of  the  heating  at  the  E-region  (90- 100km)  where  electron 
runaway  can  be  induced  resulting  in  substantial  modification  of  the  Pederson 
conductivity. 


1.  INTRODUCTION 

Generating  ULF/VLF/ELF  waves  by  utilizing  the  ionosphere  as  an  active  medium 
is  an  exciting  prospect.  The  technique  provides  frequency  agility  and  avoids  many  of  the 
geopolitical  and  economic  difficulties  associated  with  large  and  inefficient  ground  based 
facilities.  By  adjusting  the  HF  beam  geometry,  the  low  frequency  waves  can  be  made  to 
propagate  upwards  in  the  magnetosphere  for  use  in  active  magnetospheric  stimulation 
experiments  or  downwards  into  the  earth  ionosphere  waveguide  for  communication  and 
geological  probing  purposes.  One  of  the  mechanisms  discussed  for  downconverting  HF 
power  to  UUF/ELF/VLF  power  by  interactioning  with  the  ionosphere  requires  the  presence 
of  ambient  ionospheric  currents,  such  as  the  auroral  or  equatorial  electrojets  (Geimantsev  et 
al.  1974,  Kotik  and  Trakhtengerts  1975,  Stubbe  and  Kopka  1977,  Chang  et  al.  1981, 
Ferraro  et  aT.  1982),  while  another  relies  on  pure  plasma  non-linearities  and  does  not 
require  the  presence  of  ionospheric  currents  (Papadopoulos  et  al.  1982,  1983. 
Papadopoulos  and  Chang  1985,  Ko  et  ai  1986,  Ganguli  et  al.  1986).  Although  both 
mechanisms  have  been  verified  experimentally,  the  predominance  of  the  experiments  dealt 
with  the  modulation  of  ionospheric  currents  because  this  technique  demands  lower  HF 
transmitter  power  and  antenna  directionality. 

Generation  of  ULF/ELF/VLF  waves  by  modulating  ionospheric  currents  has  been 
confirmed  for  both  the  equatorial  and  the  auroral  electrojet  in  experiments  conducted  in  the 
USSR  (Migulin  and  Gurevich  1985,  Belyaev  et  al.  1987),  the  Max-Planck  Tromso  facility 
(Stubbe  et  al.  1982,  Barr  and  Stubbe  1984a,  b,  James  et  al.  1984,  Barr  et  al.  1985)  and  the 
U.S.  (Ferraro  et  al.  1982,  Ferraro  1988,  Ferraro  et  al.  1988).  Most  of  the  experiments 
were  performed  using  HF  frequencies  in  the  range  of  2-5  MHz,  while  the  power  density  at 
the  interaction  region  varied  between  Kh4  -  10*3  W/m2.  The  most  exhaustive  studies  were 
performed  at  the  Tromso  Max  Planck  facility  and  the- typical  resuits  arc  summarized  in  Fig. 
1.  These  results  arc  in  general  terms  consistent  with  results  produced  at  other  facilities, 
(Ferraro  et  al.  1988,  Belyaev  et  al.  1987)  although  the  precise  values  of  the  detected  field 
amplitudes  depend  on  local  conditions,  characteristics  of  the  HF  facility  and  other  specific 
factors.  The  purpose  of  this  paper  is  to  explore  the  potential  and  the  limitations  of  the 
ionospheric  generation  of  ULF/ELF/VLF  waves  using  as  a  guide  the  current  experimental 
results  in  conjunction  with  theoretical  extrapolations.  Namely,  starting  from  the  current 
experimental  and  theoretical  status,  we  develop  scaling  laws  and  use  them  to  explore  the 
efficiency  with  which  HF  power  can  be  convened  to  the  desired  low  frequency  range.  It  is 
expected  that  the  conclusions  will  be  useful  in  guiding  future  facilities  and  experiments. 
For  specificity  we  selected  for  our  study  the  ELF  frequency  range  (i.e.  50  Hz  -  500  Hz), 
since  in  this  range  the  interpretation  of  the  resultslis  not  complicated  by  waveguide 


resonances,  high  order  modes  and  skin  layer  thickness.  The  plan  of  the  paper  is  as 
follows.  The  next  section  summarizes  the  status  of  understanding  of  ELF  generation  by 
current  modulation  and  the  resulting  scaling  laws.  Section  3  discusses  the  level  of 
conductivity  modification  as  a  function  of  the  incident  power  density  and  altitude.  Contrary 
to  previous  studies  which  used  fluid  analysis  to  calculate  the  electron  heating,  a  completely 
kinetic  treatment  is  used  here.  Such  a  treatment  is  required  for  power  densities  exceeding 
10'3  W/m2.  Based  on  the  analysis  of  Sections  2  and  3,  Section  4  discusses  techniques  for 
increasing  the  HF  to  ELF  conversion  efficiency.  The  final  section  summarizes  the  results 
and  pinpoints  high  leverage  research  issues. 

2.  IONOSPHERIC  ELF  GENERATION-SCALING  CONSIDERATIONS 

Relating  the  amplitude  of  the  ELF  waves  on  the  ground,  to  the  extent  and 
characteristics  of  the  modified  region  in  the  ionosphere  and  to  the  design  characteristics  of 
the  HF  transmitter  is  a  complex  problem.  Aspects  of  the  problem  have  been  examined  by 
many  authors  (Kotik  and  Trakhtengerts  1975,  Belustin  et  al.  1977,  Tripathi  et  al  1982, 
Fejer  and  Krenzien  1982,  Barr  and  Stubbe  1984a,b).  It  involves  a  series  of  sequential 
steps.  First,  based  on  the  transmitter  characteristics  (power,  gain,  H.F.  frequency  and 
modulation  frequency)  and  the  applicable  model  of  the  ambient  ionosphere  the  conductivity 
modulation  is  computed  as  a  function  of  altitude.  This  allows  for  the  calculation  of  the 
modulated  current  altitude  profile,  for  a  specified  ambient  ionospheric  electric  field 
(assumed  or  measured).  This  calculation  involves  computation  of  the  primary  current,  the 
polarization  currents  required  to  set  up  quasineutraiity,  and  the  induction  currents  caused  by 
the  time  dependence  of  the  magnetic  field.  Third,  from  the  modulated  current  profile  the 
excitation  and  propagation  of.  the  ULF/ELF/VLF  waves  in  the  earth’s  ionospheric 
waveguide  can  be  computed  by  using  either  the  extended  source  technique  (Tripathi  et  al. 
1982)  or  with  an  equivalent  ULF/ELF/VLF  moment  in  the  ionosphere  the  reciprocity 
principle  (Galejs  1968,  1971;  Barr  and  Stubbe  1984a,b).  A  comprehensive  analysis  of 
these  factors  lies  beyond  the  scope  of  this  paper.  Within  the  experimental  and  theoretical 
uncertainties  our  purpose  can  be  accomplished  by  starting  from  the  current  experimental 
results  in  the  frequency  range  of  50-500  Hz  and  analyzing  them  under  the  assumption  that 
the  power  generated  is  proportional  to  the  square  of  the  dipole  moment 

MsIL  (1) 

where  I  is  the  total  modulated  current  contributing  to  the  ELF  field  on  the  ground  and  L  is 
the  linear  size  of  the  modulated  region.  All  other  factors  entering  the  efficiency  calculation 
will  be  taken  from  the  current  experimental  data  basis.  For  concreteness  we  use  as  a 


baseiine  the  measurements  and  analysis  of  the  Tromso  results  (Barr  and  Stubbe  1984a). 
These  results  air  similar  to  the  ones  reported  from  HIPAS  (Ferraro  et  al  1988).  The  range 
of  50-500  Hz  has  been  selected  in  order  to  avoid  effects  associated  with  waveguide 
resonances  which  arise  for  frequencies  above  1  kHz. 

The  experimental  results  indicate  that  when  the  heating  transmitter  operated  at  a 
power  level  of  150  MW  ERP  the  amplitudes  of  the  ELF  field  measured  on  the  ground  are 
-100  mV/m  or  1  pT.  These  correspond  to  about  10-100  mW  of  radiated  ELF  power  in  the 
200-500  Hz  range  (Barr  and  Stubbe  1984a)  The  equivalent  radiating  horizontal  dipole  at 
75-80  km  altitude  in  the  ionosphere  is  IL  »  3-5  xlO4  A-m  which  corresponds  to  an 
equivalent  ground  based  electric  vertical  dipole  with  IL  =  2-4X103  A-m.  The  polarization  is 
consistent  with  predominance  of  Hall  current  modulation.  For  the  frequency  range, 
computations  based  on  a  fluid  model  indicate  that  for  an  assumed  electric  field  Eo  =  25 
mV/m  the  peak  values  of  the  modulated  current  density  are  in  the  range  of  10*8  A/m2  and 
are  located  between  75-80km  in  altitude.  The  effective  horizontal  radiated  current  moment 
M  can  be  estimated  by  height  integration  as  . 

•  M  =  EL  =AjAzL2  =  AoAzEqL2  (2) 

where  Aj  is  the  modulated  current  density,  Az  is  the  extent  of  the  effective  radiating  layer,  in 
altitude  and  Aa  is  the  modulated  conductivity.  For  the  typical  experimental  conditions  L  » 
20km  (i.e.  at  a  heater  beam  width  of  7.5°),  Az  =*  10km,  Eo  =  25  mV/m .  the  value  of  IL  =» 
3-5X104  A-m  corresponds  to  a  peak  value  of  Aa  =  3-4  x  103  sec*1  This  is  achieved  with 
an  incident  HF  power  density  of  2  mW/m2  at  75-80  km  height.  Finally  it  should  be  noted 
that  the  value  of  IL  is  independent  of  frequency  and  the  frequency  dependence  in  the  ELF 
power  is  attributed  to  the  scaling  of  the  excitation  efficiency  (Galejs  1971). 

Since  the  quantities  that  can  be  controlled  from  the  ground  are  the  value  (and 
possibly  the  altitude  location)  of  the  conductivity  modulation  Aa  and  the  size  L,  it  is 
instructive  to  cast  the  experimental  data  as  interpreted  by  Barr  and  Stubbe  (1984a)  in  the 
form  of  Fig.  2.  Notice  the  important  scaling  valid  for  each  ELF  frequency,  i.e. 

Pelf^M2  -  (Aa)2L4  (3) 

For  instance,  if  we  increased  the  peak  conductivity  modulation  by  a  factor  of  10  while 
keeping  L=»  20km,  the  radiated  power  at  500  Hz  would  increase  from  100  mW  to  about 
10W.  Similarly,  if  we  increased  the  size  L  by  a  factor  of  5  while  keeping  A  a  =»  3-4  x  I03 
sec*1,  the  radiated  power  would  increase  by  a  factor  of  (5)4  -  6  x  102. 


A  critical  scaling  needed  for  the  determination  of  the  factors  that  optimize  Pelf  is 
the  dependence  of  the  conductivity  modification  Ac  on  the  HF  power  density  S  at  the 
modified  height. 


e  _  pHF 


(4) 


where  Phf  is  the  ground  HF  power  and  L  is  the  spot  size  at  the  appropriate  height.  Of 
course  this  neglects  absorption  at  lower  heights,  a  point  that  we  will  return  later  on. 
Assuming  that 
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(5) 


we  find  from  (3)  -  (5)  that 

PHLF-PfgL4^)  (6) 

For  a  =  1 ,  we  find  that  Pelf  is  independent  of  the  antenna  gain  and  scales  as 

Pelf -Phf  (7) 

Namely  the  HF  to  ELF  conversion  increases  as  the  square  of  the  HF  power.  The 
conversion  efficiency  depend  predominately  on  Phf  for  a  »  1,  while  for  a  «  1  more 

efficient  conversion  requires  large  spot  sizes.  A  computation  of  the  value  of  a  for  various 
altitudes  and  HF  power  densities  is  performed  in  the  next  section. 

3.  DEPENDENCE  OF  CONDUCTIVITY  MODIFICATION  ON  ALTITUDE 

AND  POWER  DENSITY 

In  this  section  we  evaluate  the  level  of  conductivity  modulation  at  various 
ionospheric  heights  as  a  function  of  the  incident  HF  power  density  S.  The  calculation  is 
"local”  and  neglects  transport.  From  the  modulated  conductivity  we  can  easily  evaluate  the 
modulated  cuirent  for  a  given  ionospheric  electric  field.  As  opposed  to  previous  studies 
(Stubbe  and  Kopka  1977,  Tomko  1981,  Chang  et  al  1981,  James  1985)  which  used  fluid 
equations  to  compute  the  variation  of  the  electron  temperature  Te  as  a  linear  approximation, 
the  present  study  uses  the  complete  time  dependent  kinetic  equation  for  the  election  energy 
distribution  function  f(e).  For  a  HF  electric  field  of  amplitude  Eo  and  frequency  coo  at  the 
chosen  height  f(e)  is  given  by  (Gurevich  1978) 


(8) 


■—  f  (e)  =  ~=.^-(e3/2D  -  L(e) 

de  OE 


D()  =  I  £§2l  v(e) 
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where  v(e)  is  the  energy  dependent  electron-neutral  collision  frequency  at  the  chosen  height 
and  ft  is  the  electron  cyclotron  frequency.  The  ±  signs  correspond  to  o  (+)  and  x  (-)  mode 
heating  correspondingly.  L(£)  is  an  operator  that  represents  the  energy  loss  due  to  various 
inelastic  processes.  It’s  form  is  discussed  in  the  Appendix  which  includes  excitation  of 
rotational,  vibrational  and  optical  levels  as  well  as  ionization  and  attachment  for  N2  and  O2. 
Note  that  the  latter  process  is  not  important  for  the  power  densities  analyzed  here.  Equation 
(8)  was  solved  numerically  for  f(£,t)  at  various  altitudes,  HF  power  density  values  S  and 
modulation  frequencies  co.  The  time  dependence  of  the  conductivity  is  found  from 
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where  aPi  ah,  az  are  the  Pederson,  Hall  and  paraiM-conducriviries  and  n  is  the  electron 
density.  As  noted  in  the  introduction  a  kinetic  analysis  is  an  absolute  requirement  for 
exploring  high  power  densities.  The  initial  f(E,  t  =  0)  was  taken  as  Maxwellian  at  Te  * 
0.025  eV.  For  the  studies  reported  below  o0  =  1.8  x  107  rad/sec  which  corresponds  to  a 
heater  frequency  of  2.8  MHz. 

We  report  below  results  for  daytime  ionospheric  conditions  corresponding  to 
altitudes  between  70-120  km.  The  ionospheric  model  used  is  shown  in  Fig.  3. 

The  need  for  a  kinetic  description  at  high  power  and  high  altitude  is  obvious  from 
the  distribution  function  shown  in  Fig.  4.  Figure  4a-c  show  the  initial  and  steady  state 
distribution  functions  for  ionospheric  heating  at  75, 90  and  100  km  altitude  and  at  values  of 
S  =  10*3,  10*2  W/m2,  correspondingly.  The  time  required  to  reach  steady  state  was  in  all 


cases  much  shorter  than  10‘3  sec.  This  implies  that  steady  state  is  reached  at  times  much 
shorter  than  the  relevant  modulation  frequencies.  It  is  seen  that  at  low  altitude  and  low 
power  density  (i.e.  75  km,  10*3  W/m2)  f(e)  does  not  deviate  much  from  Maxwellian  and 
the  fluid  description  is  a  good  approximation.  This,  however,  is  not  true  for  the  other  two 
cases  where  the  high  energy  tails  of  the  distribution  functions  become  the  dominant  part.. 

Figures  5-6  show  the  time  dependence  of  the  Hall  and  Pedersen  modulation  for  the 
above  three  cases.  At  75  km  and  10'3W/m2  the  Hall  conductivity  modulation  is  about 
7xl03sec*1  which  is  consistent  with  the  value  quoted  in  Section  2  for  the  Tromso 
experiments.  The  Pederson  conductivity  modulation  is  significantly  smaller.  This  is 
reversed  for  the  90  and  100km  cases  at  10‘2  W/m2.  The  Hall  conductivity  modulation 
becomes  progressively  smaller  and  is  negligible  at  100  km.  Furthermore  the  level  of  the 
Pederson  conductivity  modulation  is  about  3  x  104  sec*1  at  90  km  and  6  x  104  sec*1  at  100 
km.  It  is  clear  that  if  the  size  L,  which  is  controlled  by  the  antenna  gain  was  the  same  for 
all  three  cases  and  the  increase  in  the  power  density  was  entirely  due  to  an  increase  in  Phf 
by  a  factor  of  10,  the  radiated  Ppt  p  would  increase  by  a  factor  of  20  and  102  for  the  90  and 
100  km  cases  over  the  75  km  case.  Notice  that  as  noted  before  steady  state  is  established 
for  all  cases  much  earlier  than  the  low  frequency  oscillation  period,the  values  of  Act  are 
independent  of  ELF  frequency.  In  order  to  determine  the  scaling  of  conductivity 
modification  with  power  density  and  altitude,  we  conducted  a  survey  of  the  level  of  steady 
conductivity  modification  for  three  altitudes  (70, 100, 120  km)  and  values  of  S  in  the  range 
of  10-4  - 1  W/m2.  The  results  for  the  Hall  and  Pederson  conductivities  are  shown  in  Figs. 
7  and  8  The  boundary  altitudes  were  chosen  in  a  way  that  they  reflect  the  range  of 
variation  of  the  Pedersen  and  Hall  conductivities  with  altitude.  A  more  detailed  analysis  of 
the  consequences  of  the  resuits  shown  in  Figs.  7  and  8  and  of  the  efficiency  optimization  at 
intermediate  altitudes  will  be  presented  elsewhere.  We  restrict  our  discussion  here  to  the 
consequences  of  the  general  trends  derived  from  Figs.7  and  8.  These  are: 

(i)  For  low  altitudes  (~  70-75  km)  the  Hall  conductivity  provides  the  dominant 
contribution.  The  value  of  Acth  increases  very  weakly  with  power  density  (a  <  1/2)  and 
saturates  at  a  value  of  S  about  2-10*3  \s  /m2.  Further  increase  in  the  power  density  does  not 
produce  any  increase  in  the  modulated  current  density.  According  to  the  discussion  after 
Eq.(7),  optimization  of  the  conversion  efficiency  requires  increase  in  L  under  constant  S. 

(ii)  For  high  altitude  (>  90  km)  heating  the  modification  of  the  Pederson 
conductivity  dominates.  The  value  of  Aop  increases  as  S^i.e.  a  =  2),  up  to  power  density 
of  10*2  W/m2  and  saturates  slowly  afterward .  There  is  an  obvious  premium  in  increasing 
Phf>  since  in  this  case  Pelf  -  P*hf»  while  keeping  S  =  10'2W/m2. 

(iii)  Maximum  value  of  Aa  is  achieved  for  high  aidtude  heating. 


4.  EFFICIENCY  OF  HF  TO  ELF  CONVERSION 

In  this  section  we  combine  the  results  of  Barr  and  Stubbe  (1984)  as  shown  in  Fig. 

2  with  the  results  reported  in  Section  3,  and  use  them  to  determine  the  HF  to  ELF 
conversion  efficiency  and  techniques  by  which  it  can  be  optimized.  Based  on  Hall 
conductivity  modulation  of  the  polar  electrojet  at  70*75  km  altitude  Barr  and  Stubbe  (1984) 
estimate  a  power  conversion  efficiency  of  5-10  mW  per  MW  of  HF  at  200  Hz.  This 
implies  an  overall  conversion  efficiency  of  about  IQ*8,  compared  with  10'6  conversion 
achieved  by  the  Wisconsin  Test  Facility.  For  low  altitude  heating  the  results  of  Section  3 
indicate  that  increasing  the  power  does  not  have  any  significant  effect  in  the  ELF  power. 
However,  as  shown  in  Eq.  (3)  for  constant  power  density  and  therefore  constant  Ac,  Pelf 
-  L4.  An  increase  in  area  while  maintaining  the  same  approximately  power  density  can  be 
achieved  by  using  phasing  to  sweep  the  antenna  beam  over  an  area  spanned  by  a  maximum 
tilt  of  8m  degree  in  each  direction  at  a  rate  faster  than  the  cooling  rate,  which  for  75km  is 
approximately  50-100  (isec.  The  Ppt  f  then  increases  according  to  Eq.  (3)  by  a  factor  of 
(sin0m/sin0o)4  where  0O  is  the  antenna  half-width.  Taking  Qn  g  7,5°  and  9m  »  35°.  we 
find  an  increase  in  Pelf  by  a  factor  of  400.  Namely  modified  the  to  sweep  angle  to  a  30  0 
cone  at  the  Tromso  or  HIPAS  facility  over  50-100  nsec  is  expected  to  produce  2-4W  and 
have  an  efficiency  of  10^,  comparable  to  the  Wisconsin  Test  Facility. 

Alternately,  higher  efficiency  as  well  as  higher  Pelf  can  be  produced  by  high 
altitude  heating.  If  HF  heating  localized  in  the  95-100  km  region  can  be  achieved  with 
power  density  of  the  order  10  mW/m2,  the  resultant  conductivity  modification  according  to 
Figs.  7  and  8  increases  by  a  factor  of  100  over  that  for  low  altitude  heating.  According  to 
Eq.  (3)  this  will  produce  a  factor  of  104  more  power  than  the  5-10  mW ,  which  results  in 
50-100  W  in  ELF  even  in  the  absence  of  sweeping.  The  system  efficiency  increases  by  a 
factor  of  103  giving  an  overall  conversion  efficiency  of  10*5.  Incorporating  a  similar  sweep 
as  for  the  low  altitude  heating  will  result  in  further  increase  of  the  efficiency  by  a  factor  of 
400  giving  an  overall  efficiency  better  than  4  x  10-3.  The  practical  difficulty  in  realizing  the 
high  altitude  scheme,  especially  under  day  time  conditions,  is  the  fact  that  self  absorption  at 
lower  heights  could  prevent  the  achievement  of  power  densities  of  10  mW/ra2  at  95-100 
km  altitude  for  frequencies  of  2.8  MHz  used  incur  calculations.  Such  heating  can  be 
achieved  by  using  higher  frequencies,  beating  two  HF  waves  at  the  local  plasma  frequency, 
or  using  short  pulses  that  allow  the  power  to  "sneak  through"  to  high  altitudes.  These 
possibilities  are  currently  under  study.  The  enormous  increase  in  efficiency  achieved  in 
high  altitude  heating  places  a  large  premium  in  realizing  them  even  if  their  efficiency  is  by 
an  order  of  magnitude  lower. 


5.  SUMMARY  AND  CONCLUSIONS 

We  presented  a  detailed  kinetic  study  of  local  heating  of  the  ionospheric  plasma  by 
modulated  radiowaves,  of  the  resultant  conductivity  modulation  and  of  the  associated 
modulated  current  density  in  the  presence  of  an  ionospheric  electric  field  in  the  ELF  region 
(50  -  500  Hz).  Regimes  where  the  conductivity  modulation  is  a  strong  or  weak  function  of 
the  incident  HF  power  density  were  identified.  Combining  these  results  with  the  recorded 
observations  in  this  ELF  region  and  assuming  that  to  zero  order  the  ELF  power  density  is 
proportional  to  (A a  )2  L4  we  found  that  the  HF  to  ELF  power  conversion  efficiency  can  be 
increased  by  more  than  two  orders  of  magnitude  if  the  heater  could  be  swept  over  a  30° 
cone  on  timescales  faster  than  the  plasma  cooling  rate  (~  100  ^isec)  for  low  altitude  heating. 
Heating  techniques  that  can  preferentially  deposit  their  energy  at  higher  altitude  (-  90-100 
km)  where  the  dominant  modulated  current  is  the  Pederson  current,  can  further  increase  the 
efficiency  by  a  factor  of  104,  although  part  of  this  increase  can  be  negated  by  the  potential 
inefficiency  of  high  altitude  heating  as  well  as  by  a  more  inefficient  coupling  to  the 
waveguide.  We  are  currently  examining  these  issues  theoretically  and  expect  to  resolve 
them  in  combination  with  the  Penn  State  HIP  AS  experimental  campaigns 
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FIGURE  CAPTIONS 

Fig.  1  Typicai  experimental  results  in  the  ELF  region  generated  by  the  Max  Planck 
Tromso  facility  (Stubbe  et  ai.,  1982). 

Fig.  2  ELF  paper  vs  frequency  for  the  Tromso  facility  as  determined  by  the  analysis  of 
Barr  and  Stubbe  ( 1984a).  The  ELF  power  scale  has  been  multiplied  by  (3xlQ3/Ao) 
x  (20  km/L)4  to  emphasize  the  scaling  with  size  L  and  conductivity  A  a.  For  the 
Tromso  results  this  factor  is  unity  (i.e.  A a.»  3  x  103,  L  =■  20  km). 

Fig.  3  Ionospheric  profile  used  in  the  Fokker-PIanck  runs. 

Fig.  4  Initial  and  steady  state  electron  distribution  for 

(a)  75  km,  S  =  10-3  N/m* 

(b)  90  km,  S  =  10*2  W/m2 
(3)  100  km,  S  =  10'2  W/m2 

Fig.  5  Hall  conductivity  modification  for  the  cases  of  Fig.  4. 

Fig.  6  Same  as  Fig.  5  for  the  Pederson  conductivity. 

Fig.  7  Hall  conductivity  modulation  vs  S  for  70, 100  and  120  km. 

Fig.  8  Pederson  conductivity  modulation  vs  S  for  70,  100,  and  120  km. 


APPENDIX 

The  inelastic  term  of  the  Fokker- Planck  equation  (Eq.  (8))  is  composed  from  the 
following  terms  (Gurevich  1978) 

(i)  The  excitadon  of  rotadonal  levels  is  taken  into  account  by  a  term 


s^-Li^XfL 


vfo)] 


(Al) 


where  Rr(v)  =  8{50Co  Ntf/mv,  a0  =  8nao2/15,  ao  =  h2/me2  is  the  Bohr  radius,  B0  =  2.5  x 
IQ4  eV  is  the  rotational  constant  and  Nn  is  the  N2  number  density. 

(ii)  The  excitation  of  vibrational  and  optical  levels 

Sv,o(f)  =  2/mv  X  Nn[(£  +  £k)  fo  (e  +  £k)  Ck(e  +  £k)-  effete)]  (A2) 
k 


where  e*  is  the  energy  of  the  k-state  and  c*  is  the  cross  section  for  that  state, 

(iii)  Attachment  . 


Sa  =  -  N0vCTaf 

where  aa  is  the  attachment  cross  section  and  N0  is  the  O2  number  density. 

Although  the  code  includes  ionization  processes  the  energy  range  for  the  powers  under 
consideration  in  this  study  did  not  energize  the  electrons  sufficiently  to  initiate  it 
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Abstract 


Low  frequency  electromagnetic  waves  have  been  repeatedly  generated  when 
the  flow  of  the  polar  electrojet  current  is  modified  by  amplitude  modulated 
heating  from  ground  based  HF  transmitters.  The  modified  current  flow  is 
equivalent  to  a  low  frequency  plasma  antenna  known  as  Polar  Electrojet  (PEJ) 
antenna.  Current  results  indicate  that  the  radiation  efficiency  of  the  PEJ  in 
the  ELF  range  is  by  almost  two  orders  of  magnitude  lower  than  the  ground  based 
ELF  facility.  It  is  shown  that  by  matching  the  operation  of  the  HF  heater  to 
the  nonlinear  response  of  the  ionospheric  plasma  the  efficiency  of  the  PEJ  can 
be  increased  by  almost  three  orders  of  magnitude,  thus  making  the  PEJ  a  unique 
plasma  device  for  studying  many  geonhysical  problems. 
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1.  Introduction 


The  ionosphere  is  an  active,  strongly  non-linear  plasraa  medium  containing 
free  energy  sources  in  the  form  of  ambient  currents  and  density  gradients 
maintained  by  the  input  of  the  solar  energy  and  the  earth's  rotation.  Tapping 
these  free  energy  sources  is  a  long  time  goal  of  space  plasma  physics.  The 
best  known  non-linear  effect  caused  by  the  ionospheric  plasma  is  the  so  called 
"Luxemburg  effect"  in  which  an  amplitude  modulated  high  power  HF  wave 
prooagating  through  the  ionosonere  causes  modulation  transfer  (cross 
modulation)  to  other  waves  passing  through  the  same  ionospheric  medium.  The 
ionospheric  plasraa,  thus,  acts  as  a  non-linear  frequency  transformer. 

Manifestations  of  the  frequency  transforming  properties  of  the  ionosphere 

have  been  confirmed  in  a  series  of  experiments  conducted  during  the  past 

decade  in  the  U.S.*”\  the  and  Norway^  The  majority  of  the 

experiments  demonstrated  downconversion  of  HF  power  in  the  1.5-10  MHz  region 

to  the  VLF-ELF-ULF  range  (i.e.  10  kHz  -  1  mHz).  The  most  comprehensive 

experiments  were  performed  with  the  "Heater"  facility  at  Tromso,  "orway 
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ooerated  by  the  Max  Planck  Institute  .  A  summary  of  the  ground 

measurements  of  the  magnetic  field  amplitude  vs.  frequency  are  shown  in  Fig. 

,,  1 6  17 

An  analysis  of  the  ’’eater  results’  '  gave  conversion  erficiencies  that 
varied  from  10  at  a  frequencv  of  2  kHz  to  10  at  200  Hz.  Due  to  mechanical 
resonances  of  their  diesel  generators  data  ware  not  taken  in  the  20  -  200  Hz 
regime.  As  seen  from  Fig.  I  the  field  amplitude  of  the  1-20  Hz  waves  is 
comparable  to  the  200  Hz  amplitude.  The  technique  of  generating  low  frequency 
'  .1)  waves  from  HF  waves  has  been  called  "Combination  Frequency  Signal"  (CFS) 
in  the  'J.S.S.R.  literature  because  the  low  frequency  ..  could  be  selected  by 
ising  two  carrier  HF  waves  with  frequencies  ^  such  that  ^  =  u 

instead  of  amplitude  modulating  the  HF  carrier..  The  "Heater"  experimenters 


In  this  note  we 


have  called  the  technique  "Polar  Eleccrojet"  (PEJ)  antenna, 
will  use  the  name  PEJ. 

Electromagnetic  (era)  waves  in  the  ELF-ULF  (1-500  Hz)  range  are  excellent 
diagnostic  tools  for  geophysical  probing  due  to  their  long  underground 
penecration  range.  In  fact  low  frequency  eo  techniques  are  often  used  in  oil 
exploration  instead  of  seismic  methods.  Their  large  skin  depth  in  sea  water 
combined  with  their  low  attenuation  rate  (1  db/Mm)  made  ELF  the  major 
techniaue  for  communicating  with  submarines.  Finally,  waves  in  the  ELF/ULF 
range  are  excellent  diagnostic  probes  of  the  magnetosphere. 

In  attempting  co  generate  ELF/ULF  waves  from  the  ground,  one  is  faced 

with  the  problem  of  the  long  antenna.  The  free  space  wavelength  of  these 

waves  is  well  in  excess  of  1000  km,  requiring  antennae  of  totally  unrealistic 

lengths  and  accompanied  by  major  ohmic  losses.  In  the  ELF  range  a  Navy 

transmitter  named  ELF  located  in  Michigan  uses  a  frequency  modulation  centered 

at  16  Hz  for  almost  worldwide  communications  with  submarines.  The  transmitter 

has  very  low  efficiency  (<  10  f>)  using  over  2  MW  of  ground  power  to  produce 

o p 

less  than  2  W  of  radiated  ELF  power"  .  The  PEJ  antenna  offers  a  simple 

_  a 

solution  to  the  size  issue.  Due,  however,  to  the  low  efficiency  s'i.e.  10  ) 

nuoted  above,  it  has  been  discarded  as  a  potential  ionospheric  plasma  device 
and  is  only  thought  of  as  a  diagnostic  tool  of  D-region  ionospheric 
conditions.  It  is  the  purpose  of  this  comment  to  demonstrate  that,  while  the 
theoretical  analysis  of  the  PEJ  efficiency  is  correct,  the  conclusion  that  PEJ 
constitutes  an  inefficient  plasma  device  is  incorrect.  As  will  be  shown  below 
all  tbe  current  PEJ  antennae  failed  to  match  their  operating  procedures  to  the 
non-lir.ear  ionospheric  plasma  response.  Adjusting  their  operation  to  match 
the  ionopsheric  plasma  response  results  in  an  efficiency  increase  of  almost 


three  orders  of  magnitude.  As  a  resuLt  PEJ  can  serve  as  a  unique  ionospheric 
plasma  device  for  studying  a  broad  class  of  geoohysical  problems. 

2.  Physics  of  the  PEJ  Antenna 

PEJ  antennae  have  been  studied  systematically  over  the  last  ten  years  and 
a  quantitative  understanding  of  its  physical  principles  well  validated  by 
experiments  has  emerged.  The  free  energy  source  for  its  operation  is  the 
auroral  current  system.  This  system  is  powered  by  the  flow  of  the  solar  wind 
plasma  intercepting  the  polar  magnetic  field  lines.  The  resultant  emf  drives 
currents  along  the  magnetic  field  lines  of  the  collisionless  raagnetospheric 
plasma  (Birkeland  currents)  which  close  across  the  magnetic  field  lines  at  an 
altitude  between  70-100  km  in  the  collision  dominated  ionospheric  plasma. 

This  DC  cross  field  current  is  called  the  auroral  electrojet  current  and  is 
confined  around  the  auroral  oval.  The  auroral  oval  is  a  ring  around  the 
geomagnetic  pole,  where  the  ionospheric  conductivity  is  enhanced  by  the 
continuous  precipitation  of  keV  electrons.  The  total  power  of  the  electrojet 
current  system  is  10  MW  and  its  ohmic  dissipation  results  in  an  increase  of 
the  ionospheric  temperature  by  100°  K.  The  electrojet  current  I  is  driven  bv 
an  ionospheric  electric  field  EQ  which  results  from  mapping  along  the  magnetic 
field  lines  of  the  emf  induced  by  the  solar  wind.  The  average  observed  value 
of  ZQ  is  about  25  mV/m,  although  at  times  it  can  reach  150  mV/m  or  be  as  low 
as  5-10  mV/m. 

The  PEJ  antenna  operates  on  the  following  principle.  An  HF  transmitter 
araolitude  modulated  at  a  frequency  or,  equivalently,  two  HF  transmitters 
with  frequencies  =*  .1  are  incident  on  the  lower  ionosphere.  Absorption 

of  the  HF  power  by  the  electrons  that  carry  the  electrojet  currents  results  in 
local  heating  and  modifies  the  temperature  dependent  value  of  the  local 


conductivity.  roc  values  of  .1  such  that  .. r  <<  1,  where  r,  is  the  electron 
■  c 

cooling  tine,  the  conductivity  modulation  reflects  the  shape  of  the 

envelope  imposed  by  the  heater  modulation  (for  the  lower  ionosphere 

t  o.l  ns).  The  associated  variation  of  the  current  density  uj  is  be  given 
c 

by 


ij.  *  Ao  •  E  (1) 

This  constitutes  a  primary  source  of  modified  current.  For  low  frequencies 

the  total  current  is  divergence  free.  As  a  result  a  secondary  current  is  set 

up  consistent  with  tr.e  divergence  free  condition.  The  total  current  consists 

of  the  current  I ^  flowing  through  the  heated  volume,  a  current  I2  flowing 

around  the  modified  volume,  and  a  field  aligned  current  1 3  flowing  through  the 

magnetosphere  to  the  conjugate  ionosphere.  The  equivalent  circuit  shown  in 
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Fig.  2  was  proposed  by  Stuhbe,  et  al.  to  describe  the  radiative  loop.  It  Is 
easy  to  see  for  the  ELF  range  (10-500  MHz)  the  radiating  current  is 
essentially  due  to  I,  and  only  the  R^,  ^  ^  part  of  the  circuit  is  active.  For 
the  !JLF  part  (1-10  Hz)  It  becomes  relevant.  For  our  frequency  range,  the 
field  aligned  currents  do  not  respond  ( ^.^  =  50  H)  and  the  radiating  currents 
are  essentially  due  to  the  locally  modified  currents.  We  refer  the  reader 
interested  in  more  details  to  Refs.  1,  12-20,  23. 

For  our  purposes  it  is  sufficient  to  consider  that  the  ELF  source  in  the 
lonopshere  is  an  horizontal  magnetic  dipole  source  with  moment 


M  =  IL 


(2) 


where  I  is  the  value  of  the  nodulated  current  and  L  Is  the  horizontal  linear 


size  of  Che  modified  region,  "or  Che  high  latitude  case 

I  *  AjAzL  3  E  AaAzL  (3) 

o 

where  Aj  is  the  modified  current  density  and  az  the  vertical  extent  of  the 
effective  radiating  layer.  Combining  (2)  and  (3),  we  find  that 

M  3  (E  Az)AoL2.  (4) 

o 

The  value  of  the  ELF  field  on  the  ground  depends  on  the  value  of  M  and  the 
efficiency  at  which  the  field  couples  to  the  earth  ionosphere  waveguide 
(excitation  factor).  For  the  ELF  range,  the  excitation  factors  varies  as  £7. 
f72.  The  ELF  power  on  the  ground  would  then  be 

Pr._  -  (E  Az)2  e2U)Ao2L4  (5)  _ 

LLr  O 


where  e(fl)  is  the  waveguide  excitation  factor^-  The  HF  to  ELF 

conversion  efficiency  n  is  then  given  by 


where  is  the  average  ground  based  HF  power.  For  the  case  analyzed  by  Barr 
and  Stubbe ^ ,  a  square  heating  pulse  was  used.  The  value  of  was  1  MW 
and  at  the  ELF  generation  region  the  value  of  L  -  10  km  corresponding  to 
heater  beam  width  of  7.5°.  — 


3.  Si 


itj  of  Current  Observations 


An  analysts  of  the  PEJ  efficiency  has  been  presented  in  Barr  and 
Stubbe ^ .  The  purpose  of  our  comment  is  to  examine  if  the  efficiency  n  can 
be  increased  over  the  one  determined  in  the  above  reference  by  matching  the 
operation  characteristics  of  the  HF  heater  to  the  nonlinear  response  of  the 
ionosphere.  We  thus  write  Ea.  (6)  as 


n 


K 


Ao2L4 

PHF 


(7) 


where  K  is  a  constant  that  incorporates  factors  beyond  our  controL  (i.e.  EQ, 

_z,  ;  etc).  Furthermore  in  order  to  be  specific,  we  consider  the  particular 

case  of  200  Hz,  (3  =*  1.2  x  iO2  sec”1). 

The  particular  observations  were  made  using  an  HF  frequency  f  »  2.76  MHz, 

effective  radiated  power  ERP  =  150  MW,  and  antenna  beam  width  2aQ  =»  15°;  the 

"•8 

estimated  ELF  power  was  6-8  mW,  giving  an  efficiency  of  n  *»  10  .  The 

—8 

observations  were  consistent  with  a  modified  Hall  current  density  of  10 
A/m",  which  corresponds  to  a  conductivity  modification  of  uo0  =  A  x  10  mhos 
at  an  altitude  hQ  =  78-80  km,  and  a  value  of  Eq  »  25  nV/m.  The  value  of 
corresponding  to  the  modified  conductivity  was  saturated  at  T^  =  2000  °K. 

The  oower  density  at  78-80  km  was  of  the  order  of  2  mW/m“.  For  the  ELF/ULF 
;reauencv  range  of  interest,  the  observed  ELF  power  was  consistent  with  a 
horizontal  dipole  in  the  ionosphere  with  magnetic  moment 


'.{  =  3.5  x  10  A-ra. 


i  /.  75 

As  is  well  known  from  the  work  of  Galejs"  ’  horizontal  dipoles  are 
inefficient  exciters  of  the  earth  ionosphere  waveguide  at  100-200  Hz,  having 


excitation  efficiency  z  of  the  order  of  10  “  compared  to  that  of  comparable 

ground  baaed  electric  dipole.  Notice  that  ground  based  ELF  sources  are  also 

—  -  -3 

very  inefficient  with  tvpical  efficiencies  e  of  the  order  of  10  or  less  for 

horizontal  magnetic  dipoles.  rrom  Eq.  (7),  the  measured  efficiency  and  the 

values  discussed  above,  we  can  rewrite  Eq.  (7)  as 


op  a 

,„-8  ,La  '  ,rHF(0,  L  , 
n  -  10  i— ;  p  *  ;  I—; 


■  Ac 


(8) 


HF 


where  -o3  ■  4  x  I0-7  mhos,  LQ  =  hQ  tan  =»  10  km,  and  P^p  Q 


values  of  the  reference  case.  The  on-off  heater  time 


‘Ho 


1  MW  are  the 
is  half  of  the  wave 


oeriod  of  200  Hz.  Furthermore ,  we  should  remember  that  the  _c  tonification 

2 

over  an  L~  area  was  accomplished  with  an  incident  power  flux  So  =*  2  nW/m  at 


V 


4.  Modulated  HF  Heating  and  Associated  Conductivity  Modulation 

The  nodulated  heating  of  the  ionosphere  by  an  X-raode  HF  field  of 
frequency  j  and  peak  amplitude  EQ  is  given  by 


AT 

ill!,  nr.  l  ,,  _  -  v 

2  dt  "en  .  2,2  ;er'‘e  o 

1  +  v  /  w 
en  o 


v  exD  ,  f 
ev  •  L 


T  -  2000 

a 

2000  T 


(9) 


where  is  the  nuiver  energy  defined  as 


2P2 
e  E 

_ o 

2  2 
m  Lu 

o 


(10) 


and  ^  '-3  the  effective  frenuency 


jj 

o 


(11) 


Also  -,la  Is  the  electron  cyclotron  frequency,  .■  is  the  electron  neutral 

collision  frequency  for  momentum  transfer,  v  the  rate  of  inelastic  electron 

neutral  collisions  resulting  in  excitation  of  the  rotational  levels  of  N0,  the 

terra  with  v  .  is  the  vibrational  loss  and  T  =  100°X  Ls  the  ambient 
ev  o 

temperature  of  the  neutrals.  For  the  range  of  electron  energies  of  interest 
here 


,-9 


ven  2,3  X  10  ^  1'100UK;  =  Vq  '100 


(12a) 


-n  ?  “1/2  T  -1/2 

, ,  r  3  ")  v  10  'J  .  =  V  - - -  , 

er  '  '100 ;  "  F.  '100^' 


(12b) 


u  *  3.5  x  10-3  M  =  1068  v  * 
ev  X 


where  N  is  the  number  density  of  the  neutrals.  Defining 


v  =  T  /T  ,  T  =  v_t 
e  o  R 


V  /  uj 
0  0 


we  can  cast  Eg.  (9)  in  a  dimensionless  form  as 


(13) 


2 

2  dr 


,  2  2 
l  +  a  y 


v  -  1 

*7 


10. 68  exp  [f 


(14) 


where  f  ■  5.3  +  3.755  tanh(0,ll  (y  -  18)].  The  modification  of  the 
temperature  results  in  changes  in  the  conductivity  of  the  medium.  For  the 
Lonosphere  the  conductivity  tensor  o  transverse  to  the  magnetic  field  is 
given  by 
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Ju , 
V 


(15) 


were  Op  is  the  Pedersen  conductivity,  i.e.  the  conductivity  along  the  electric 
field  and  is  the  Hall  conductivity.  In  our  dimensionless  units  they 
are  given  by 


(16a) 

(16b) 

(17) 

(18) 


From  the  above  equations  we  note  that  the  values  of  the  dimensionless 
quantities  a  and  3,  7  and  t  are  altitude  dependent .  while  the  vaiue  of  K 
depends  only  on  the  incident  HF  power  density.  Equations  (9-18)  can  and  have 
been  solved  numerically  under  several  conditions.  Important  insight  is 
provided  by  considering  some  limiting  cases. 

(i)  The  regime  of  interest  corresponds  to  3  «  1,  i.e.  neutral  densities 
N  <<  10*^  ?/Cn',  so  chat  for  y  =  l ,  <jp/a^  3  3  «  1.  Furthermore 
heating  always  decreases  forcing  the  Hall  current  to  flow  outside 
the  heated  volume.  The  Pederson  conductivity  has  a  different 
deoendence  on  y.  For  values  of  dy  <  l ,  Op  increases  with  v,  producing 


currents  in  antiphase  with  the  Hall  currents,  while  for  _•/  >  l  it 
decreases  with  y  and  the  currents  are  in  phase  with  the  Hall 
currents.  Notice  chat  since  the  ambient  electrojet  current  flows  at  an 
angle 


5  *  tan  *  1/3 

with  respect  to  the  ambient  it  is  a  Hall  current.  For  iy  =  1  it 
flows  at  a  45°  angle  with  respect  to  This  implies  that  the  induced 

oolarization  electric  field  has  the  same  value  as  E<).  Finally  for 
£y  >>  1  the  Pedersen  current  becomes  dominant. 

(ii)  Inspection  of  the  heating  equation  (Eq.  14)  shows  the  presence  of  a 
strong  barrier  for  values  of  y  a  20-30  due  to  the  third  tern  on  the 
r.h.s.  of  Eq.  (14).  This  term  is  due  to  the  high  value  of  the 
inelastic  electron-^  cross  section  for  vibrational  excitation.  It 
implies  a  heavy  energy  penalty  in  achieving  temperatures  past  y  »  20. 

(iii)  Equation  (14)  shows  that  for  values  of  ay  <  1  and  v  <  20  the  heating 
rate  is  exponential  with  t  while  the  cooling  rate  which  in  this  regime 
is  due  to  rotational  energy  transfer  has  a  weaker  (~~)  dependence. 

These  observations  will  be  utilized  in  the  next  section  in  our  attempt  to 

optimize  the  efficiency  of  the  PEJ  antenna. 

5.  Efficient  Operation  of  the  PEJ  Antenna 

In  order  to  be  concrete  we  consider  the  benchmark  case  discussed  in 

Section  3.  Following  Barr  and  Stubbe^  ^  we  assume  that  the  source  height  is 

14  3 

at  73-80  km,  which  corresponds  to  a  neutral  density  N  «  2  x  10  '/cm  . 

For  this  value  of  density  a  ■  .05  and  £  *  .05.  The  value  of  the  quiver  energy 
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(19) 


in  °K  is  given  by 

T  m  vi  !L ^2_  Ov 

.  >  t  j 

—  mW/m~  to 

For  f  *  2.76  MHz,  ;2g/2i:  =  1.35  MHz,  ERP  of  150  MW  which  corresponds  to 

2  mW/m"  and  T  »  100  °K,  z/T  °  .22.  Using  Eq .  (13),  we  find  that 
o  o 

K  ■  30.  Finally  for  ELF  frequency  of  200  Hz  the  dimensionless  ELF  time  is 

"CLF  =  Figure  3a  shows  the  waveform  of  y  for  a  square  wave  heater 

operation  with  period  half  the  period  of  the  ELF  wave.  It  clearly 
demonstrates  the  inefficient  fashion  in  which  the  heater  power  is  utilized. 
First,  more  than  95%  of  the  power  utilized  to  maintain  the  temperature  at  its 
saturated  value,  while  the  absorbed  energy  is  transferred  to  the  excitation  of 
vibrations  which,  of  course,  do  not  contribute  to  the  conductivity 
modification.  Second,  we  are  overheating  the  plasma.  There  is  little,  if 
anything,  to  gain  by  heating  past  the  value  of  y  *  1/3,  as  is  obvious  from 
Eqs.  (16).  Figure  3a,  however,  also  suggests  the  operating  procedure  which 
can  lead  to  improved  efficiency,  which  is  discussed  next. 

Let  us  consider  a  facility  with  the  same  average  power  as  the  Tronso 
"Heater"  facility  but  equiped  with  the  capability  to  scan  fast  over  an  area  up 
to  a  cone  of  37.5°  (Fig.  4).  The  ratio  of  this  area  to  that  modified  by  the 
7.5°  half -beamwidth  HF  antenna  is 


tan  ^5 Q  2 
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If  using  the  same  average  power  as  the  "Heater"  facility  we  could  modify  the 
conductivity  of  the  accessible  area  A  at  the  desired  ELF  frequency  the 
efficiency  according  to  Eq.  (8),  will  increase  by  a  factor  i  which 
corresponds  to  3.3  x  10^.  It  will,  however^  be  reduced  by  the  relative 


•nodi f icacton  of  conductivity.  The  quadratic  dependence  c,~  on  the  area  Is  a 
consequence  of  the  fact  that  a  larger  area  produces  a  simultaneous  increase  In 
total  current  as  well  as  flux.  Techniques  by  which  this  can  be  accomplished 
are  demonstrated  below. 

We  next  compute  the  power  enhancement  over  the  standard  case  that  can  be 
accomplished  by  dividing  the  area  A  into  a  number  of  smaller  areas  Aq  and 
sweeping  the  transmitter  so  that  during  the  half  cycle  ^£Lp/2  “  ^  it  stays  on 
each  spot  only  A0/A  of  the  time.  The  expected  enhancement  R  of  the  ELF  power 
over  that  of  the  standard  case  is  given  by  the  ratio  of  the  efficiencies 
obtained  from  Eq.  (7)  for  the  two  cases.  The  antenna  with  beamwidth  of  15° 
Irradiates  each  area  A  for  the  duration  Tqjj  with  a  revisit  time  tgpp  and  so 
the  total  heated  area  is  A  =*  ( +  'oFF^  Ao*  Further  the  antenna  gain  K 
is  reduced  from  the  value  of  30,  the  heated  area  is  increased  by  the  factor 
30/K,  Thus  for  the  same  power  PHp,  we  get 
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where  -a  is  the  final  conductivity  modification  averaged  over  all  the  spots 

and  U  =  .72  is  the  conductivity  modification  of  the  standard  case  (Fig. 

3a).  The  heated  area  can  be  increased  by  choosing  small  tgjj  and  Tgpp,  whose 

limits  are  set  by  the  heating  and  cooling  timescales.  The  heating  timescale 

is  given  by  Eq.  (14),  which  for  y  >  1  yields  y  3  y( x  3  0)  exp(t/i0)  with  tq  = 

3/2K.  To  achieve  a  20Z  rise  in  the  temperature  in  one  heating  cycle  we  need 

=<  .  2  t  ,  which  is  0.01  for  K  =  30.  A  value  1 0„  smaller  than  this  will 
‘ON  o’  ON 

lead  to  a  weak  heating.  Since  the  idea  of  the  proposed  scheme  is  to  avoid  the 
strong  vibrational  losses  by  keeping  Tg  <  2000°K,  the  relevant  cooling  time  is 
that  due  to  the  rotational  losses.  In  the  normalized  units  this  has  a  value 
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~  l  and  so  igpp  should  be  <  1  to  prevent  excessive  cooling.  With  these 

considerations  different  cases  are  examined  as  follows. 

We  first  take  the  case  where  the  full  power  (K  =  30)  of  the  antenna 

irradiates  each  area  Aq  over  a  time  =■  .02  with  a  revisit  time  T0Fp.  The 

procedure  was  repeated  three  times  till  t  -  17.  The  enhancements  for 

different  values  of  TqFF  are  given  in  Table  la,  which  shows  that  increasing 

the  total  area,  i.e. ,  the  ratio  tqff/  T0N  resu^-ts  an  increase  of  R.  This  is 

however  accompanied  by  a  decrease  in  Ao  because  the  distribution  of  the 

power  into  a  much  larger  area  limits  the  heating.  This  latter  effect 

overcomes  the  advantage  of  the  increase  in  area  and  beyond  7gFF/ =  40  the 

value  of  R  decreases.  There  is  thus  a  balance  between  these  two  effects  that 

gives  an  optimum  value  of  tgFF  for  a  given  Tg^.  The  value  of  _o  is  the 

average  obtained  in  the  last  (  +  tqff)  before  the  power  was  turned  off  at 

r  *  17.  The  peak  enhancement  for  this  case  was  455  corresponding  to 

t0FF'/T3N  3  40  and  C^e  ^eacin?  of  the  temperature  of  a  typical  spot 

is  shown  in  Fig.  3b.  With  the  same  K  value  and  different  values  of  tgjj,  the 

enhancement  R  will  be  ootinum  at  different  values  of  *_/*,.  The  reason 

OFF  ON 

for  the  relatively  low  enhancement  compared  to  the  theoretical  maximum  of  3  x 
3 

10  ,  is  due  to  the  reduction  in  the  conductivity  modulation,  giving  a  low 
_  2 

value  for  the  (Ac/Aa)"  factor. 

Figures  3c  and  d  indicate  the  results  of  a  trade  off  study  in  which  the 
antenna  gain  was  reduced  bv  a  factors  of  2  and  3  so  that  the  area  of  each  spot 
was  2Aq  and  3Aq  respectively.  This,  of  course,  results  in  a  reduction  of  the 
corresoonding  value  of  K.  In  the  case  presented  in  Fig.  3c,  the  antenna  gain 
Ls . reduced  bv  a  factor  of  2  so  that  K  ■  15  and  the  area  of  each  spot  is  2A0 . 
This  requires  a  longer  tg*j  co  achieve  significant  heating  of  a  soot  and  also 
the  total  number  of  spots  that  can  be  covered-is  reduced  (Table  lb).  Unlike 
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Table  la,  which  shows  the  enrtancenent  for  a  wide  range  of  *Qpp/’  r^,  Table  lb 

shows  che  sensitive  dependence  of  R  on  "oFF,/ T0N  arount*  tne  case  of  maxitaum 

enhancement.  Figure  3d  shows  the  case  with  K  *  10  where  the  maximum 

enhancement  524  in  efficiency  was  achieved.  In  this  case  only  15  spots  each 

with  area  3A^  were  heated.  In  all  the  three  cases  K  =  10,  15  and  30  shown  in 
o 

Figs.  3b, c  and  d  enhancements  ~  500  are  easily  achieved  and  the  waveforms  of  y 
in  these  cases  are  similar  to  each  other. 

6.  Summary  and  Conclusions 

In  this  comment  we  nave  demonstrated  that  operating  the  PEJ  antenna  in  a 
node  consistent  with  the  nonlinear  properties  of  the  ionospheric  plasma  can 
increase  che  HF  to  ELF  conversion  efficiency  by  a  factor  of  500  efficiently 
than  the  current  facility,  without  any  increase  in  the  ground  HF  power  or  the 
antenna  gain.  Furthermore,  this  is  close  to  one  order  of  magnitude  more 
efficient  than  the  Navy  ELF  facility.  Under  strong  auroral  conditions  the 
efficiency  can  increase  dramatically.  Comparable  or  even  higher  efficiencies 
can  be  achieved  for  other  ELF  frequencies.  To  summarize,  we  believe  that  the 
current  experi-mencal  evidence  indicates  that  the  PEJ  antenna  constitutes  the 
first  practical  device  based  on  tapping  the  free  energy  of  the  ionospheric 
olasoa.  Adjusting  its  operation  to  the  high  efficiency  mode,  requires  only 
ohase  shifting  or  other  electronic  means  for  scanning  the  HF  antenna  at  a 
relatively  high  soeed.  Vith  this  complement  PEJ  can  be  a  unique  eeophysical 
exoloration  tool. 
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Figure  Captions 


Figure  1 


Figure  2 


Figure  3 


Figure  4 


A  tvpical  magnetic  field  amplitude  vs.  frequency  measured  on 
ground  at  the  Max  Planck  Troraso  "Heater"  facility  (Stubbed. 


Equivalent  circuit  system  to  describe  the  modulation  current. 
Current  I ^  is  flowing  through  the  modified  volume,  I2  around  the 
-  modified  volume,  along  the  magnetic  field  lines  through  the 
conjugate  ionosphere,  and  I  is  the  current  flowing  in  the  wave 
( Stubbe  st  al.‘^). 

The  ionospheric  electron  temperature  (in  units  of  100  °K)  under 
various  scenarios  of  heating  the  PEJ  and  producing  ELF.  a:  All 
powar  delivered  to  an  area  corresponding  to  a  beam  half-width  of 
7.5°.  b-d:  The  power  is  delivered  over  a  wider  area  consisting 
of  many  spots.  R  is  the  enhancement  of  the  ELF  power  over  that 
of  the  case  a. 

Scheme  for  the  enhancement  of  ELF  power  by  sweeping  over  an  area 
within  the  beaming  angle  of  37.5°. 
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Figure  Caption 


'able  l 


Enhancement  of  the  ELF  generation  for  different  antenna  gains  (K) 
and  ’oFF'' t0N  rati°s*  a*  The  enhancement  R  for  K  =  30  and  a  wide 
range  of  values  of  x gpp/ Tq^.  The  gain  from  increasing  tgpp/ 
is  offset  by  the  accompanying  decrease  in  La,  giving  the  peak, 
value  of  R  =  455.  b.  For  half  the  full  antenna  gain  (K  =*  15)  the 
variation  of  the  enhancement  around  the  peak  value  (R  ■  458)  is 
shown. 
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